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1.0 


OVERVIEW 


"  This  report  constitutes  the  first  phase  of  research  conducted  by 
Axiomatix  wherein  our  aim  was  to  derive,  evaluate  and  compare  various  detec¬ 
tor  structures  whose  purpose  Is  to  Intercept  a  spread -spectrum  communication 
transmitter.  The  transmitter  under  surveillance  employs  a  variety  of  modu¬ 
lation/spreading/transmission  techniques  which  are  Invariably  assumed  to 
emit  the  message- bearing  signal  In  deep  background  noise  or  Interference. 

Thus,  the  challenging  task  of  the  Intercepting  detector  Is  to  reveal  as  best 
he  can  any  transmission  of  the  unfriendly  spread  waveform  In  the  presence 
of  a  strongly  obscuring  noisy  environment. 

Naturally,  the  degree  of  success  achieved  by  the  Interceptor  when  de¬ 
tecting  the  presence  of  the  communicator's  spread  signal  depends  on  the  amount 
of  information  available  to  him  regarding  the  structure  of  that  signal.  At  one 
extreme,  the  interceptor's  most  fortunate  situation  would  be  to  acquire  the 
spreading  code  Itself  ("crack"  the  code).  This  being  toe  demanding,  he  must 
settle  for  less,  such  as  approximate  knowledge  of  the  signal’s  carrier  centet 
frequency,  code  rate,  code  epoch,  spreading  bandwidth,  etc.,  or  a  subset  thereof. 
At  the  other  extreme,  he  might  know  almost  nothing*;  In  which  case,  he  could 
resort  to  a  simple  energy  discriminator  device  ("radiometer")  since  It  Is,  on 
many  occasions,  a  low-SNR  asymptotically  optimal  detector.  In  this  report,  our 
main  interest  is  to  investigate  the  possibility  that  the  Interceptor  could  do 
better  by  optimally  processing  whatever  Information  Is  available  to  him  between 
the  two  extremes  cited  above. 

The  detectability  of  the  primary  candidates  for  low-probabll Ity-of- 
intercept  (LPI)  waveforms,  such  as  direct-sequence  (DS),  frequency-hopping  (FH), 
time-hopping  (TH)  and  their  hybrids,  is  customarily  related  to  two  factors: 

(1)  performance  level  of  the  interceptor's  detector  whose  only  function  Is  to 
monitort  the  communicated  messages  ajid,  (2)  amount  of  signal  processing  ( 1 - e. , 
complexity)  associated  with  such  performance.  Thus,  the  exploitablllty  of  a 
spread- spectrum  waveform  Is  measured  In  terms  of  the  complexity  required  of 
the  interceptor's  receiver,  so  that  it  can  perform  such  surveillance,  as  well 
as  in  terms  of  the  level  of  its  power  spectral  density  (PSD)  and  the  resulting 
probability  of  Intercept. 

* 

Approximate  spectral  band  location  Is  a  minimum  prerequisite  for  any 
interceptor. 

^Intercepting/ jamming  combinations  would  also  be  of  interest  at  a  more 
advanced  stage  of  study. 


In  subsequent  sections,  it  will  be  shown  that  the  detectability 
(exploitability)  of  a  spread  signal  is  greatly  dependent  on  the  interceptor's 
knowledge  regarding  the  form  of  both  the  waveform  and  the  interference.  In 
other  words,  one  should  know  not  only  the  spreading  format  used  (plus,  possib¬ 
ly,  some  other  parameters*),  but  it  is  also  very  helpful  to  be  able  to  identi¬ 
fy  the  kind  of  background  noise  or  interference  involved.  Although  most  of 
this  report  deals  with  white  Gaussian  noise  as  the  dominant  form  of  (omni¬ 
present)  interference,  a  case  is  made  in  the  final  section  about  the  impact 
of  nonwhite  (random  tones,  in  particular)  interference.  It  will  be  shown  in 
Sections  3  through  5.1  that,  when  white  noise  is  the  only  deterrent,  the  gains 
achieved  by  intelligent  processing  of  the  received  waveform  can  vary  from 
modest  to  significant  (always  as  compared  to  the  performance  of  a  radiometer), 
depending  on  the  scenario  at  hand.  However,  when  random  interfering  tones 
are  present  and  a  DS  or  a  hybrid  FH/DS  is  detected,  those  gains  become  impres¬ 
sive  (many  tens  of  decibels,  for  example).  For  details,  see  section  5.2. 

In  searching  for  those  structures  which  perform  the  aforementioned 
tasks,  we  start  from  the  optimal  solutions,  as  derived  from  applying  the  opti¬ 
mal  1 ikelihood- ratio  rule.  The  results  assume  a  good  deal  of  knowledge  about 
the  signal  parameters  (although  never  the  spreading  code  Itself)  and  could  thus 
become  rather  academic  in  most  practical  situations.  Furthermore,  they  typically 
suffer  from  the  common  symptom  of  prohibitive  implementational  complexity.  Thus, 
it  is  important  to  consider  suboptlmal  structures  and  evaluate  their  relative 
loss  with  respect  to  the  (unattainable)  optimal  performance.  This  has  Indeed 
been  the  spirit  which  permeated  the  present  work. 

The  report  is  organized  as  follows:  In  section  2,  we  classify  the 
various  receivers  and  briefly  lay  the  theoretical  groundwork  for  evaluating  the 
forthcoming  detector  structures.  In  Sections  3-  4  and  5,  we  develop  the  theory 
for  DS,  TH  and  FH  signals,  respectively^.  This  Is  also  the  order  by  which  the 
relative  gains  from  sophisticated  receivers  seem  to  increase  when  faced  by  white 
noise.  In  particular,  a  two-reception,  synchronous,  coherent  DS  detector  can 
result  in  a  gain  of  up  to  4.5  dB  in  signal-to-noise  ratio  (SNR)  above  the  radi¬ 
ometer.  Each  of  the  three  aforementioned  features  of  this  detector  contributes 
approximately  1.5  dB;  however,  each  also  imposes  strenuous  conditions  on  the 

See  Section  2  for  a  detailed  discussion. 

^Let  us  note  that  the  hybrid  FH/DS  of  section  5.2  can  also  be  thought 
of  as  a  DS  signal  in  colored  interference. 


system  implementation.  This  kind  of  modest  relative  benefit  must  be 
attributed  to  the  "noise-like"  appearance  of  a  high-rate  DS  code,  which  makes 
it  hard  to  distinguish  from  wideband  thermal  noise.  On  the  other  hand,  TH  or 
FH  can  look  spectrally  different  from  OS,  so  that  higher  gains  can  be  expected; 
this  is  documented  in  Sections  4  and  5.  We  note,  however,  that  those  gains  are 
typically  associated  with  excessive  complexity;  In  which  case,  implementation 
becomes  the  cardinal  issue.  This  prompted  us  to  Investigate  suboptimal ,  but 
still  very  efficient,  schemes  which  are  suggested  by  today's  technology  (e.g., 
real-time  autocorrelation  devices).  The  application  of  such  novel  Ideas  in 
Section  5  has  established  conditions  (both  by  analysis  and  simulation),  under 
which  one  can  anticipate  significant-to- impressive  gains  over  the  radiometer, 
even  in  the  presence  of  white  noise.  It  is  shown  in  Section  5  that  a  proper 

measure  of  performance  improvement  (defined  in  section  5.1)  is  proportional  to 

2 

yh  ,  where  yh  Is  the  hop  SNR  of  an  FH  system.  The  proportionality  constant  de¬ 
pends  on  algorithmic  parameters,  but  is  independent  of  the  assumed  tactical 
scenario.  Thus,  by  properly  translating  yh  into  physical  parameters  (i.e., 
transmitter/intercepter  distance,  signal  power,  antenna  gains,  etc.),  one  can 
arrive  at  certain  threshold  values  in  the  proximity  of  which  the  signal  becomes 
highly  detectable.  Although  this  would  be  the  bottom  line  of  a  global  inter¬ 
ception  analysis,  this  task  will  not  be  undertaken  here,  as  it  is  highly  sce¬ 
nario  dependent. 

In  conclusion,  the  findings  of  this  report  can  serve  as  broad  analyt¬ 
ical  guidelines  for  the  determination  and  evaluation  of  a  variety  of  detectors 
in  the  context  of  LPI  Interception.  They  also  serve  to  motivate  a  further  pur¬ 
suit  of  techniques  which,  in  light  of  today's  (and  tomorrow's)  technology, 
emerge  as  very  appealing,  sophisticated  and  promising  alternatives  to  the  sim¬ 
plistic  solutions  of  the  past. 


2.0 


DETECTOR  CLASSIFICATION  AND  PERFORMANCE  MEASURES 


In  this  section,  we  first  classify  a  variety  of  receivers  which  an 
interceptor  could  employ  to  detect  the  presence  of  a  wideband  (spread)  signal 
imbedded  in  white  and/or  nonwhite  Interference.  Furthermore,  we  discuss  some 
approximate  measures  of  performance  upon  which  the  assessment  for  any  detec¬ 
tor's  merit  is  based.  We  should  note  here  that  such  measures  give  a  quanti¬ 
tative  feeling  of  what  can  be  expected  from  certain  detectors  but  do  not  pro¬ 
vide  any  indication  of  the  structure's  implementational  complexity,  which  is 
definitely  an  important  parameter.  Therefore,  a  total  evaluation  of  the  sys¬ 
tem  effectiveness  should  incorporate  that  as  a  separate,  but  indispensable, 
consideration. 

The  classification  attempted  here  is  sumnarized  in  Figure  2.1.  Wide¬ 
band  detectors  can  first  be  classified  as  optimal  or  suboptimal .  We  term  "op¬ 
timal"  those  structures  which  result  from  a  straightforward  application  of  the 
generalized  likelihood  ratio  theory  to  this  problem  and  Invoke  a  number  of  as- 
suptions  regarding  the  signal  structure.  For  instance,  optimal  solutions  typi¬ 
cally  assume  knowledge  of  certain  parameters  (code  rate  or  hopping  rate,  SNR, 
carrier  frequencies,  etc.)  and  average  over  the  unknown  ones  (timing  epochs, 
true  hopping  or  code  sequences,  etc.).  Since,  however,  optimal  structures  are 
often  hard  to  implement,  one  resorts  to  "suboptimal"  solutions  whereby  either 
one  (or  more)  assumptions  are  removed  or  certain  "atypical"  nonlinearities,  such 
as  the  £n  cosh(*)  function,  are  substituted  by  simpler  ones  as,  for  instance,  a 
square-law  device. 

A  second  classification  results  from  the  nature  of  the  spread- 
spectrum  signal  for  which  the  intercepting  device  is  designed,  so  it  can  be 
a  DS,  FH,  TH  or  hybrid  detector.  Also,  DS  and  TH  detectors  can  be  coherent 
or  noncoherent,  depending  on  whether  or  not  the  carrier  phase  is  known.  FH 
detectors  are  assumed  to  be  noncoherent  from  hop  to  hop  due  to  the  nature 
of  the  communicating  channel  and  the  transmitter  wideband  FH  synthesizers. 

The  presence  or  absence  of  timing  (epoch)  information  about  the  spread¬ 
ing  code  distinguishes  detectors  as  synchronous  or  asynchronous.  It  is  clear 
that  timing  is  initially  unavailable  since  the  detector  does  not  even  know  if 
the  signal  is  there;  hence,  asynchronous  structures  naturally  attract  more  prac¬ 
tical  interest.  However,  synchronous  and/or  coherent  detectors  will  also  be 
considered--not  only  because  they  provide  useful  upper  bounds  in  performance--but 
also  because  of  the  conceptual  possibility  of  improved  detector  structures. 


Wideband  Detectors  for  LPI  Waveforms 


Figure  2.1.  Classification  of  LPI  Detectors 
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aided  by  the  recursive  estimation  of  pertinent  parameters  [8].  Furthermore,  by 
quantizing  the  continuous-epoch  uncertainty  region  into  an  adequate  finite  num¬ 
ber  of  points,  matched-detector  structures  can  be  devised  that  are  synchronized 
to  each  of  those  levels  and  whose  outputs  are  combined  to  yield  the  final  deci¬ 
sion  variable.  Thus,  since  the  maximum  code-epoch  uncertainty  equals  a  chip 
time*,  one  concludes  that  good  nonsynch ronous  detectors  can  be  constructed  from 
combinations  of  a  few  synchronous  ones.  More  discussion  on  this  topic  is  pro¬ 
vided  in  later  sections. 

Finally,  we  classify  the  detectors  according  to  whether  they  expect  to 
operate  in  the  presence  of  white  (and  typically  Gaussian,  such  as  thermal  noise) 
or  nonwhite  interference.  The  maximum-likelihood  approach  renders  itself  to  the 
first  case  most  easily,  while  its  formulation  is  cumbersome  for  the  second.  Ac¬ 
cordingly,  the  theory  is  fairly  mature  for  AWGN  interference,  while  it  is  still 
at  the  exploratory  level  for  scenarios  which  call  for  more  complicated  interfer¬ 
ence  forms. 

Let  us  now  briefly  discuss  analyt'-al  ways  for  measuring  performance 
in  the  present  context.  It  is  well  known  [3]  that,  under  the  common  detection- 
performance  criteria  (most  notably,  the  Neyman-Pearson  philosophy),  generalized 
likelihood  ratios  yield  optimal  hypothesis-testing  solutions,  and  performance  is 
measured  by  the  resulting  pair  of  detection  and  false-alarm  probabilities  (Pp, 

PpA) .  Furthermore,  it  can  be  shown  that,  in  many  cases  of  practical  interest, 
the  overall  likelihood  ratio  (LR)  can  be  expressed  as  the  product  of  individual 
LR's  due  to  the  statistical  independence  of  successive  code  chips  (random  DS 
code)  or  hopping  slots  (TH  or  FH).  When  the  true  model  for  the  problem  does  not 
allow  for  such  assumptions  (as  in  the  case  of  a  random,  but  fixed,  carrier  phase), 
the  resulting  LR  is  too  complicated  for  either  analysis  or  implementation^.  The 
performance  of  such  decision  schemes  is  assessed  by  bounding  arguments  (see  the 
following  sections)  so  that  the  product  LR  can  always  be  thought  of  as  the  ana¬ 
lytical  cornerstone.  Since  the  log-likelihood  ratio  (LLR)  is  not  only  theoret¬ 
ically  equivalent,  but  also  more  convenient  to  implement  than  an  LR  (it  turns 
products  into  sums),  it  will  serve  as  our  point  of  departure. 

The  number  of  terms  entering  the  summation  of  the  LLR  is  typically 
large--it  equals  the  total  number  of  chips  (DS)  or  time  slots  (TH)  or  frequency 

slots  (FH)  observed.  Therefore,  it  is  reasonble  to  conclude,  via  a  central-limit 
* 

This  clearly  assumes  a  purely  random  code  sequence  for  which  a  full 
chip  time  shift  corresponds  to  another  realization  of  the  same  stochastic  process, 
t 

In  addition,  it  can  be  argued  that  the  superposition  of  random  data 
on  DS  or  TH,  such  as  PSK  modulation,  invalidates  the  notion  of  a  constant  carrier 
phase. 


A  A  A  A.V  A  A  A  A  .A  .A  _  A  A  A  A 
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type  of  argument,  that  the  distribution  of  the  resulting  LLR  is  approximately 
Gaussian  under  either  hypothesis.  This  constitutes  the  basic  assumption  upon 
which  performance  is  derived  in  this  report.  In  other  words,  the  pair  (P_,P_.) 
is  expressed  in  terms  of  the  Gaussian  integral  function 


Q(x) 


(2-1) 


and  its  inverse  ()"*(•);  furthermore,  using  the  approximation  that  the  LLR  pos¬ 
sesses  the  same  variance  under  either  hypothesis,  it  is  easily  shown  that  per¬ 
formance  is  quantified  by  the  distance  d,  defined  as 


d 


,-l 


(2-2) 


which  is  a  function  of  the  decision  rule  and  the  system  parameter  values. 

Peterson  et  al  [3,  section  4.9]  have  invoked  the  Gaussian  assumption 

2 

about  the  LLR  in  order  to  show  that  d  is  equal  to  the  difference  of  the  means 
of  the  LLR  under  the  two  hypothesis  (hence,  the  notion  of  distance)  and  Is  also 
equal  to  its  variance.  Most  important,  d  can  be  related  to  the  original  product 
LR  A  via* 

d2  =  £n^A|2  HQ})  =  £n(l  +  var{A|H0}|  (2-3) 

where  the  second  equality  results  from  the  fact  that  <£  {AIH0I  *  1  for  every  LR. 
Therefore,  (2-3)  measures  the  performance  of  any  LLR*,  which  could  Include  arbi¬ 
trary  memoryless  nonlinearities  (indeed,  the  optimal  ones).  Since  many  of  them 
are  hard  to  implement,  however,  it  is  of  interest  to  examine  the  penalty  In  per¬ 
formance  incurred  when  going  from  optimal  to  suboptimal  decision  rules  just  by 
simplifying  (or  approximating)  those  nonlinearities.  Clearly,  the  suboptimal 
rules  are  also  associated  with  some  distance  d,  as  per  (2-2),  which  we  shall 
denote  by  dappr  (approximate)  as  opposed  to  the  1  t-xact"  distance  dfix  derived 
from  (2-3).  The  other  approximations  nonwithstanding,  the  difference  between 
^appr  anc*  ^ex  s*10u^cl  signify  the  loss  paid  for  implementational  simplicity.  In 
addition,  dappr  will  serve  as  the  yardstick  when  comparing  the  approximate  rules 
resulting  from  various  other  assumptions  (or  the  l*ck  thereof). 

_  ★ 

Here,  <?{•}  and  var{*}  stand  for  the  statistical  mean  and  variance, 
respectively. 

^Always  within  the  Gaussian  assumption. 
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3.0  DETECTION  OF  DIRECT-SEQUENCE  SIGNALS  IN  THE  PRESENCE  OF  AWGN 

In  this  section,  we  consider  and  analyze  wideband  detectors  for  DS 
waveforms.  It  will  be  assumed  throughout  that  the  receiver  has  knowledge  of 
the  carrier  frequency  and  code  rate  and  that  the  code  is  biphase  modulating 
the  ccrrier.  For  the  carrier  phase,  the  two  possibilities  of  it  being  known 
(coherent  detection)  and  unknown  (noncoherent  detection)  are  examined  separately 
in  sections  3.1  and  3.2,  respectively.  A  common  assumption  in  both  of  these 
sections  is  that  the  detectors  are  synchronous,  i.e. ,  the  code  chip-timing 
epoch  is  known.  In  most  practical  situations,  this  Is  rather  unrealistic;  the 
synchronous  results  can  then  be  thought  of  as  upper  bounds  on  the  performance 
of  any  asynchronous  detector.  Furthermore,  by  quantizing  the  continuous-epoch 
uncertainty  region  into  an  adequate  finite  number  of  points,  matched-detector 
structures  can  be  devised  which  are  synchronized  to  each  of  these  levels  and 
whose  outputs  are  combined  to  yield  the  final  decision  variable.  More  discus¬ 
sion  on  the  cost  of  asynchronism  Is  given  in  section  3.3.  Section  3.4  discusses 
some  gains  derived  from  combining  two  independent  receptions,  ranks  the  perfor¬ 
mance  of  the  aforementioned  detectors  and  compares  them  to  that  of  the  radiometer. 

In  order  to  introduce  some  notation,  let  the  high-rate,  ±l-valued 
random  spreading  code  c(t)  be  represented  by 


c(t) 


1  cn  p(t  -  nTc  -  eTc) 


(3-1) 


n=-<*> 


where  p(t)  is  a  unit  pulse  of  duration  Tc  seconds  and^n}^  is  a  sequence  of 
independent,  identically  distributed)  (iid)  random  variables  (rv's),  with 
Pr|cn  =  1^  =  Prjcn  ■  -l|  =  0.5.  Furthermore,  the  chip  epoch  is  modeled  by  the 

The  waveform  observed  by  the  detector 


rv  e,  uniformly  distributed  over  (0,1) 
is  therefore  given  by 


r(t)  = 


.  j 


/Z5  c(t)  cos(u0t  +  +  n(t) 

n(t) 


TO 

TO 


(0  <  t  <  T)  (3-2) 


-\J 


V'j 


.  -f, 


& 


[5 
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* 


i 


C 


e 


c 


where  S,  wg  and  $  are  the  average  signal  power,  carrier  radian  frequency  and 
carrier  phase,  respectively,  under  H^,  and  n(t)  is  bandpass  AWGN  with  one-sided 
power  spectral  density  (PSD)  of  Ng  W/Hz.  The  observation  time  is  T  seconds, 
which  we  assume  to  be  an  integer  multiple  of  the  chip  time,  i.e.,  T  «  NTC;  N  a 
positive  integer.  Such  an  assumption  is  the  least  restrictive  since,  in  prac¬ 
tice,  N  is  generally  a  large  number  (of  the  order  of  hundreds  or  higher).  How¬ 
ever,  a  somewhat  stronger  restriction  is  embedded  in  (2),  namely,  the  assumption 
that,  under  hypothesis  Hj,  the  signal  is  present  during  the  whole  observation 
interval.  That  excludes  the  possibility  of  the  signal  either  starting  or  ending 
at  any  random  time  in  (0,T).  Still,  such  a  formulation  Is  important  because  It 
provides  meaningful  and  fairly  simple  comparative  conclusions  which  (a)  would 
otherwise  be  obscured  by  mathematical  complexities  and,  (b)  can  be  argued  to 
extend  to  more  general  models. 

In  terms  of  the  above  notation,  a  synchronous  detector  implies  that  e 
is  identically  zero,  while  a  coherent  one  means  that  *  is  known.  The  detector 
for  which  ♦  is  unknown  but  constant  over  the  observation  Interval  (0,T)  will  be 
called  carrier  noncoherent*.  For  computational  purposes,  we  also  consider  the 
fictitious  chip-noncoherent  detector,  where  phase  is  assumed  to  change  randomly 
from  chip  to  chip.  Although  this  is  totally  unrealistic.  It  is  discussed  here 
because  its  readily  derived  performance  serves  as  a  very  useful  lower  bound  to 
the  performance  of  the  optimal,  carrier-noncoherent  receivers,  which  are  diffi¬ 
cult  to  analyze. 


3.1  Synchronous  Coherent  Detectors 

When  both  the  random  phase  <p  and  random  chip  epoch  e  are  assumed  known, 
the  detector  is  asked  to  perform  the  following  composite  hypothesis  testing  prob¬ 
lem:  Decide  between  the  alternatives  Hg  and  Hj,  where 


r(t) 


/S  c(t)  +  nj(t) 

Hj(t) 


TO 

TO 


QO  *  t  <  NT^ 


(3-3) 


Hypothesis  H1  is  composite  because  it  contains  all  possible  patterns  that  the 
code  can  assume  in  NTC  seconds.  Although  we  consider  only  random  codes  [in 

'  _ ' 1  jf 

This  clearly  assumes  that  the  observation  duration  is  well  within 
the  coherence  time  of  the  spread- spectrum  channel . 


10 


which  case,  there  are  exactly  M  3  2^  such  patterns),  some  of  the  steps  below 
would  also  be  valid  for  deterministic  codes.  The  subscript  I  in  nj(t)  indi¬ 
cates  that  only  the  inphase  component  of  the  noise  contributes,  with  a  flat 
two-sided  PSD  of  Nq/2  W/Hz.  At  this  point,  we  should  emphasize  that  the  prob¬ 
lem  formulation  is  reminiscent  of,  but  not  identical  to,  the  M-ary  communica¬ 
tion  problem  of  detecting  the  presence  of  one  of  M  =  2^  possible  waveforms  em¬ 
bedded  in  noise.  In  the  latter  case,  we  know  a  priori  that  one  of  the  signals 
is  present,  but  we  don't  know  which  one  it  is.  In  the  current  case,  our  only 
interest  is  detecting  the  presence  of  an^  one  of  the  M  equilikely  candidate 
waveforms  without  identifying  it. 

Starting  with  classical  well-known  results  about  detectors  designed 
to  minimize  the  probability  of  error  [3],  it  follows  that  the  optimal  detector 
performs  a  generalized  likelihood  ratio  threshold  comparison.  Specifically, 
the  generalized  likelihood  ratio  is 


2?  f(r(t)|H1#C<(t);0  <  t  <  T)  Hj 
aOU))  -  2- -  .7-7—.  - 1  i 


f(r(t)|H0) 


(3-4) 


where  ■F(^r( t)  ] H indicates  the  conditional  likelihood  functional  of  the  hypo¬ 
thesis  H.j ,  given  the  observation  r(t);0  <  t  <  T,  and  c^t)  is  the  jth  pattern 
out  of  the  2n  total  possible.  For  the  case  of  AWGN,  the  generalized  likelihood 
ratio  in  (3-4)  reduces  to  the  following  form: 


exp{-NYJ  2N  f  NT 

ACr(t))  =  - ^  exPyli~  /  C  r(^)  c*(t)  dt 

2"  i=l  N0  g  1 


)i-Nyr)  2N  N  H, 

V*  2  «p  TT  rj  c1j  i  Ao 

2n  i=l  N0  j=l  J  lf0  u 


(3-5) 


where 


fC  Nr 


(3-6) 


V5 


O'*  * 
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is  the  predetection  (or  chip)  SNR,  c^  is  the  jth  chip  of  the  jth  pattern,  and 


r.  Is  given  by 

J 


A  f  JTC 

rj  *  J  C  r(t)  dt 

(M)TC 


(3-7) 


The  expression  in  (3-5)  can  be  further  simplified  using  certain  symmetries  and 
identities.  Indeed,  it  is  shown  in  Appendix  A  that  (3-5)  is  expressible  as 


aOUD  -  ^“PK}cosh^rjj  | 


(3-8a) 


or,  equivalently. 


ii 

£nACr(t))  ■  -Ny  +  ][  £n  cosh 


{\ri) 


H<  £n  Aq 


(3-8b) 


We  note  that  (3-8)  could  have  been  derived  directly  based  on  the  statistical  in¬ 
dependence  of  successive  chips  and  the  AWGN  (see  [6]}.  The  somewhat  alfferent 
approach  preferred  here,  however,  is  more  easily  generalized  to  models  where 
successive  chips  are  not  independent  (c.f.  section  3.2). 

A  structure  implementing  the  rule  (3-8)  is  shown  In  Figure  3.1  with 
f(x)  =  £n  cosh(x).  Its  operation  Is  as  follows:  A  time- synchronized,  chlp-by- 
chlp,  coherent  integration  of  the  received  waveform  r(t)  results  In  a  sequence 
of  rv's  r^ ,  each  of  which  is  separately  processed  by  the  memoryless  nonlinear¬ 
ity  f(-)»  then  summed  over  the  desired  length  of  N  chips  to  produce  the  decision 
statistic.  Let  us  notice  that  the  optimal  f(*)  is  an  even  function,  which  can  be 
be  interpreted  as  the  receiver's  attempt  to  compensate  for  his  lack  of  knowledge 
regarding  the  sign  of  each  integrated  code  chip. 

The  mean  under  of  the  random  argument  of  the  £n  cosh(  •)  function 
in  (3-8b)  Is  2yc,  which  Is  also  Its  variance  under  Hq.  Therefore,  for  typical 
predetection  SNR  values  below  -10  dB  or  so,  the  approximation  £n  cosh(x)  s  x2/2 
is  applicable;  in  which  case,  (3-8)  reduces  to  the  approximate  (suboptimal)  rule 
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i 


r 


a 


S' 

< 

c 

i. 


N  ?  H, 

I  ri2  £  V 

J-l  J  ° 


(3-9) 


Since  the  aforementioned  approximation  of  the  £n  cosh(x)  function  becomes  In¬ 
creasingly  tight  as  x  +  0,  it  follows  that  (3-9)  Is  asymptotically  optimal  as 

yc  -*•  °- 

Within  the  theoretical  framework  of  Section  2,  the  performance  of 
(3-9)  is  easily  derived.  The  mean  and  variance  of  A'  are  given  by  (for  proof, 
see  Appendix  B) 


*  "(VcKj  +  5kl) 


.arf|Hk}  *  "(Vc)2  (I  +  2^c  skl) 


respectively,  where  6t1  Is  the  Kronecker  delta 


(3-10a) 


(3- 10b) 


■I 
(0  , 


1  ,  k*l 


(3-10c) 


Therefore,  Pp  and  Pp^  are  related  by 


P0  *  Q 


Q'  M-  a  ^  yc 

V1  +  bYc 


(3-11) 


a  ■  ;  b  ■  4 


(Synchronous  coherent  detector)  (3-12a) 


Since  byc  Is  typically  much  less  than  1,  (3-11)  can  be  combined  with  (2-2)  to 
yield 


d  5  a  vTT  y 
appr  'c 


( 3- 12b) 


O  O.  c*  O 


We  note  that  performance  is  dictated  by  the  product  *4T  yc»  which  Is  typical 
of  schemes  employing  post-detection  integration.  In  contrast,  detection  of  a 
known  waveform  would  involve  the  factor  J Nyc  (see  section  3.5);  the  losses  in 
performance  due  to  the  generalized  test  (since  the  signal  pattern  Is  unknown) 
for  low  values  of  SNR  are  evident  by  comparison. 

Equation  (3-11)  has  been  plotted  in  Figure  3.2  for  various  values  of 
N  and  PpA.  In  terms  of  the  overall  observation  time  T  and  the  null-to-null 
spread- spectrum  bandwidth  W  *  2T  Hz,  (3-12)  can  be  rewritten  to  establish 

)  v» 

the  required  CS/NoDreq  order  to  achieve  the  performance  level  dappr  as 

&)  s  (7)Vrdappr  (3‘13) 

'  u/req 


On  the  other  hand,  the  distance  d  for  the  exact  rule  (3-8)  can 

A 

also  be  derived.  Let  y^  =  2SS  r^./Ng.  Under  HQ,  yj  Is  a  zero-mean  Gaussian 
rv  with  variance  oy2  =  2yc-  Using  the  fact  that  * jcosh(yj)}  =  expjoy2/2^  in 
conjunction  with  (3-8a)  and  (2-3)  results  In 

dfix  =  yJH  £n  cosh^Yj  (3-14) 

2 

For  small  x,  £n  cosh(x)  *  x  /2,  which,  upon  substitution  in  (3-14),  verifies 

that  d  -»■  d  of  (3-12)  as  Yr  +  0.  In  conclusion,  approximate  rule  (3-9) 
cX  appr  c 

is  well  suited  for  the  low-predetection  SNR  range  of  interest. 


Synchronous  Noncoherent  Detectors 


Let  us  now  relax  the  unrealistic  assumption  that  the  carrier  phase  of 
the  DS  waveform  is  known  by  the  interceptor,  but  still  retain  the  assumption  of 
a  known  code  epoch.  When  the  observation  interval  is  smaller  than  the  coherence 
time  of  the  channel,  the  carrier  phase  $  can  be  modeled  as  an  rv,  uniformly  dis¬ 
tributed  in  (0,2ir),  resulting  in  the  aforementioned  carrier-noncoherent  system. 

In  addition,  a  fictitious  chip-noncoherent  system  will  be  considered  whereby  the 
phase  is  assumed  to  change  randomly  from  chip  to  chip.  The  reason  for  discussing 
such  a  totally  hypothetical  signal  Is  twofold:  (1)  it  can  be  shown  that  the 
optimal  decision  rule  for  a  carrier-noncoherent  system  is  asymptotically  equiva¬ 
lent  to  the  decision  rule  of  the  chip-noncoherent  system  (hence,  their  corre¬ 
sponding  performances  also  coincide  asymptotically)  as  the  SNR  yc  decreases  to 
zero;  (2)  more  important  is  the  fact  that  the  performance  of  the  chip-noncoherent 
structure  is  easily  obtained  (within  the  Gaussian  approximation),  in  sharp  con¬ 
trast  to  the  unwieldy  analysis  of  the  optimal  carrier- noncoherent  system.  Thus, 
the  (unobtainable)  performance  of  the  latter  can  be  nicely  bracketed  between  the 
(obtainable)  performance  of  the  former,  which  serves  as  an  asymptotically  tight 
lower  bound*,  and  the  upper  bounding  performance  of  the  synchronous  coherent  de¬ 
tector  in  the  previous  section  (3.1). 

Let  us  first  consider  a  carrier-noncoherent  system.  Under  the  usual 
narrowband  assumption  “0»  2ir/T,  the  resulting  decision  rule  is 


where  Iq(*)  is  the  zeroth-order  modified  Bessel  function  and  is  the  Uh  corre¬ 
lation  envelope 

Ri  =  V®!!2  +  \  ’  1=1,2 . ^  (3"16) 


Since  a  chip-noncoherent  detector  utilizes  less  statistical  informa¬ 
tion  than  a  carrier-noncoherent  detector,  it  is  inferior  to  the  latter  for  all 
SNR  values. 
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In  (3-16), 


'i  =  jli  S  C'J  4  ^ ;  a=,-‘ 


(3-17) 


J  Jl 


r  fcos  unt 
c  0 


(j-l)T, 


sin  wgt 


dt  ;  j=l,.. . ,N 


(3-18) 


As  mentioned,  the  optimal  rule  (3-15)-(3-18)  is  hard  to  mechanize. 
Instead,  suppose  that  the  received  waveform  r(t)  is  represented  under  by 


oo 

'(t)  =  /2s  £  Cj  p(t -jTc)cosjw0t  +  )  +  n(t) 


(3-19) 


where  ^ |  is  a  sequence  of  statistically  independent  uniform  phases,  thus  model¬ 
ing  the  aforementioned  chip- noncoherent  system.  It  is  then  straightforward  to 
show  (see  Appendix  C)  that  the  optimal  decision  rule  is 


■"**  *  X  £"  ‘0  (nf 


Hq  A0 


(3-20) 


where  r.  is  the  envelope  of  the  jth  chip 

J 


+ s 


(3-21) 


with  rj.,rn.  as  per  (3-18).  An  implementation  of  (3-20)  is  shown  in  Figure  3.3. 

J  ? 

Again,  under  the  small-argument  approximation  £nlg(x)  =  x  /4,  it  follows  that 

(3-20)  reduces  to  the  suboptimal  rule 


:  ?  Hi 

y  r  2  Z  A  * 

j-i  i  n<i  ° 


(3-22) 


Figure  3.3.  Optimal  Synchronous  Chip-Noncoherent  Detector  for  Direct-Sequence  Waveforms 
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We  note  that  the  rule  (3-22)  is  formally  identical  to  (3-9),  the  difference 
being  that  here  r^  corresponds  to  a  noncoherent  chip  integration. 

Although  the  distance  dfiX  for  the  exact  rule  (3-20)  can  be  derived 
(see  Section  4  for  a  closely  related  case),  we  shall  be  content  with  d,„„„  of 
(3-22)  in  view  of  the  low-SNR  conclusions  of  the  previous  section.  Using  the 
Rayleigh  (under  H(0  or  Rician  (under  nature  of  r.,  it  is  easily  shown  that 
(see  Appendix  D) 


4/Kl  ■  (Vc)[*  * 


(3-23a) 


var<r .  H. 


(Ve)2E+**e  sk] 


(3-23b) 


Combining  (3-22)  and  (3-23)  and  employing  the  Gaussian  assumption  results  in 
(3-11)  (or  (3- 12b)},  where  now 


a  =  1  ;  b  =  2 


(Synchronous  chip-ncncoherent  detector)  (3-24) 


Equivalently, 


(4)  .  * 


(3-25) 


which  reveals  a  loss  factor  of  SZ,  or  1.5  dB  with  respect  to  the  coherent  detec¬ 
tor.  This  is  the  penalty  paid  for  the  lack  of  coherency.  Furthermore,  since 
the  performance  of  the  carrier-noncoherent  system  is  bracketed  between  the  co¬ 
herent  (upper  bound)  and  the  chip-noncoherent  (lower  bound)  ones  that  differ  by 
only  1.5  dB,  it  is  concluded  that  any  implementational  complexity  beyond  that 
of  rule  (3-22)  is  not  justified  from  a  practical  standpoint. 


3.3 


We  shall  now  explore  the  possibility  of  removing  the  rather 
stringent  requirement  of  epoch  synchronization.  In  other  words,  we  shall 
derive  optimal  and  suboptimal  structures  that  do  not  assume  knowledge  of 
the  chip  epoch  e  (see  (3-1));  furthermore,  we  shall  estimate  the  (approxi¬ 
mate)  penalty  in  performance  compared  to  the  previously  discussed  synchron¬ 
ous  detectors.  For  reasons  which  were  explained  before,  attention  will  be 
given  to  totally  coherent  and  totally  (chip)  noncoherent  detectors.  These 
are  the  chip-by-chip  processors  which  readily  render  themselves  to  asyn¬ 
chronous  alternatives. 

In  subsection  3.3.1,  Continuous-Epoch  Uncertainty,  we  discuss  struc¬ 
tures  which  utilize  analog  filters  in  order  to  perform  averaging  over  the  con¬ 
tinuous  rv  e,  uniformly  distributed  in  £0,TC).  Thus,  the  final  system  output, 
which  is  a  continuous  function  of  time,  provides,  after  proper  sampling,  the 
theoretically  optimal  answer  to  the  generalized  (i.e.,  average  over  e)  likeli¬ 
hood  ratio  test  for  discriminating  between  Hq  and  H^.  The  derivation  presented 
here  follows  some  basic  steps  from  Krasner's  work  [4]  in  a  somewhat  more  simpli¬ 
fied  and  expanded  manner.  In  subsection  3.3.2,  Quantized-Epoch  Uncertainty,  we 
explore  suboptimal  structures  that  result  from  a  quantization  of  the  epoch  un¬ 
certainty  region  and  are  much  easier  to  implement  than  the  optimal  detector. 

In  particular,  the  two-point  quantization  scheme  is  discussed  and  analyzed  in 
detail.  It  will  be  argued  that  different  interpretations  of  that  analysis 
could  serve  as  upper  and  lower  bounds  on  the  performance  of  the  optimal  contin¬ 
uous  uncertainty  detector  in  subsection  3.3.1. 


3.3.1  Continuous-Epoch  Uncertainty 

For  the  sake  of  clarity,  let  us  consider  the  coherent  case  first. 

Most  of  the  following  steps  carry  over  without  change  to  the  chip-noncoherent 
decision  rule  (3-22)  as  well. 

The  starting  point  is  tne  coherent  likelihood  ratio  (3-8a),  addition¬ 
ally  conditioned  on  the  random  offset  t.  Following  a  procedure  identical  to 
the  development  of  section  3.1  yields  the  expression 


N  r—  <  T1  (  .  N 

ACr(t) |  e)  =  TT  [cosh(rj(e))  expj-Yc|J  =  exp|-NYcJ  Jf  cosh(>jU))  (3-26) 
J  1  3  1 


for  the  aforementioned  conditional  likelihood  ratio.  In  equation  (3-26),  r ^ ( e ) 
indicates  the  coherent  integration  of  the  received  waveform  r(t)  in  the  jth  chip 
interval,  as  per  (3-7),  where  the  integration  limits  have  now  been  adjusted  to 
reflect  the  conditioning  on  an  assumed  value  of  the  offset  e: 

r.(e)  =  /  c  r(t)  dt  (3-27) 

(j-l)T« 

In  the  following,  the  constant  factor  exp|-Nyc|  will  be  denoted  by  KA ,  for  sim¬ 
plicity.  Let  us  note  that  the  final  decision  rule  (3-8b)  for  the  synchronous 
case  results  immediately  from  (3-26)  for  e  =  0  once  we  rewrite  the  latter  in 
the  equivalent  form 


A(r(t)|0  =  Ka  exp^n^TT  cosh  (rj(e))^j  =  <A  exp|  £  £n  cosh^rj(e)j 


(3-28) 


and  use  the  monotonicity  of  the  function  exp(*).  It  will  soon  be  clear  that, 
from  an  implementational  viewpoint,  (3-28)  is  more  convenient  to  work  with 
than  (3-26). 

The  overall  conditional  likelihood  ratio  A(V(t))  is  formed  by  aver¬ 
aging  the  conditional  one,  i.e.,  ACr(t)|e},  with  respect  to  e,  which  is  assumed 
uniformly  distributed  in  £0,T 

V 

AO(t))  =  <fe  A(r(t)|0  =  Ka  r  c  exp |  £  £n  cosh^rj(e)j|  de 


0 

=  Ka  A  exp  )  |  £n  cosh  / 

0  )Jsl  \(  j~l)T 


r(-r)dx  | >de 


(3-29) 


Although  the  above  decision  rule  is  optimal  from  a  theoretical  stand¬ 
point,  it  is  not  immediately  clear  exactly  how  it  can  be  implemented  in  prac¬ 
tice  without  resorting  to  various  approximations,  such  as  quantizing  the  range 
of  e,  etc.  One  possible  way  will  now  be  illustrated  which  employs  appropriate 
analog  filters  and  devices. 
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First,  let  us  define  the  waveform  y(t)  as 


f  r( t)  d-r 


(3-30) 


The  above  function  is  easily  mechanized  as  the  output  corresponding  to  an  input 
r(t)  of  a  linear  filter  matched  to  the  chip  square  pulse 


(see  (3-1)3  since,  indeed, 

f  *>  d,  -  / 


r(x)  p(t-x)  dx 


It  is  clear  from  (3-29)  that  the  argument  of  the  £n  cosh(«)  function  is  just 
the  value  of  y(t)  at  time  t  =  jT  +e.  Thus,  we  can  rewrite  (3-29)  as 

v 


A(r(t)]  =  Ka  J  c  £n  cosh|y(jTc+e|Jde 


(3-32) 


The  next  step  is  to  somehow  create  the  summation  in  the  exponent  of  (3-32)  as 
the  outcome  of  a  linear-filtering  operation  whose  input  is  £n  cosh(y(t)}.  Such 
a  sum  of  equl spaced  (by  Tc}  samples  of  the  Input  function  arises  at  the  output 
f(t),  with  proper  sampling,  of  a  filter  whose  impulse  response  is  a  sequence 
of  equl spaced  delta  functions: 


N-l  /  x 
h(t)  =  £  5(t  “  mTc) 


(3-33) 


1 

• 

1  ,  0  <  t  <  T 

J 

p(t)  3  j 

c 

^0  ,  otherwise 

(3-31)  W 
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Indeed,  when  the  output  waveform 


f(t)  *  h(t)  ®  £n  cosh(y(t)}  =  ft  <5Gr-mTc)^  *n  cosh(y(t-x))dT 


\  in  cosh(y(t  -  mTc}) 


(3-34) 


of  such  a  filter  is  sampled  at  t  =  kTc+  e,  it  equals 


f(kTc+e)  ■  1  £n  cosh^y((k-m)Tc+e}j 


(3-35) 


or,  changing  variables  j-k-m. 


f/kT  +e)  »  S  in  cosh(y(jT +e)\ 
'  c  '  j®k-N +1  \  V  c  V 


(3-36) 


Clearly,  then,  the  exponent  of  (3-32)  arises  exactly  when  k»N,  i.e.,  when  f(t) 
is  sampled  at  t  »  NTc+e. 

Before  proceeding,  let  us  note  that  an  impulse  response  as  per  (3-33) 
corresponds  to  a  tapped-delay-1 ine  filter  whose  (N-l)  delay  elements  are  Tc  sec¬ 
onds  each  (see  Figure  3.4(a)}.  Such  an  implementation,  however,  bears  the  re¬ 
striction  that  it  is  not  recursive,  I.e.,  higher  values  of  N  require  more  delay 
elements.  An  alternative  structure  is  shown  In  Figure  3.4(b)  which  employs  a 
single  delay  element  in  the  feedback  path.  Ideally,  such  a  filter  would  provide 
an  infinite  series  of  delta  functions,  thus  alleviating  the  need  to  specify  N. 
Typically,  however,  stability  considerations  dictate  the  insertion  of  a  gain 
a  <  1  in  the  feedback  loop;  the  available  maximum  N  is  then  practically  restric¬ 
ted  y  the  infinte-impul se-response  (HR)  filter  memory. 

Returning  to  (3-36)  and  combining  it  with  (3-32),  one  recognizes  that 
the  desired  likelihood  ratio  statistic  is,  in  fact,  the  time  average  of  the 
function  exp{f(t)}  over  the  interval  ^NTC ,(N+l)Tcj .  Again,  such  an  average  is 
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created  by  sampling  at  t  s  (N+1)T  the  output  of  a  filter  whose  impulse 

v 

response  is  a  unit  pulse,  T  seconds  long.  Figure  3.5  provides  a  possible  1m- 

v 

plementation  of  such  a  "moving -average”  filter.  Finally,  the  detector  compares 
the  output  sample  at  t  «  (N+1)TC  to  a  threshold  and  produces  a  decision.  The 
overall  scheme  is  depicted  in  Figure  3.6. 

Several  observations  can  be  made  here.  First,  a  practical  scheme 
would  avoid  the  £n  cosh(*)  and  exp(*)  non! inearlties  since  they  are  hard  to 
construct  and  replace  them  by  a  single  square-law  device  operating  on  y(t), 
i.e.,  the  following  approximation  of  (3-32)  could  be  used: 

ilOtOD  =  K  Jc  f  y2  (JT  *  0  dc  (3-37) 

A  0  U-i  J 

The  basis  for  such  a  simplification  is  the  previously  used  fact  that,  for  low 
SNR,  y(t)  is  a  small  quantity;  in  which  case,  the  approximations  £n  cosh  y(t)  : 

p 

1  +  y  (t)  and  exp(f(t)3  :  1  +  f(t)  can  be  employed.  Furthermore,  let  us  note 
that  the  chip  pulse  shape  need  not  be  rectangular.  An  arbitrary  (but  chip-time 
limited)  shape  can  be  accommodated  by  designing  the  very  first  filter  of  Fig¬ 
ure  3.6  to  be  matched  to  that  particular  shape.  All  other  elements  remain  the 
same. 

Finally,  it  can  be  shown  that  the  arguments  and  detector  structures 
of  this  subsection  apply  directly  to  the  chip-noncoherent  case,  provided  that 
the  first  nonlinearity  £n  cosh(*)  of  Figure  3.6  Is  replaced  by  £n  IQ(*)  and 
the  Initial  chip-matched  filtering  at  the  Input  of  that  figure  Is  done  nonco- 
herently.  This  can  be  accomplished  by  the  noncoherent  chip-envelope  detector 
comprising  the  predetection  integration  of  Figure  3.3,  once  the  Inphase  and 
quadrature  integrate-and-sample  type  filters  are  replaced  by  analog  matched 
filters  with  a  continuous-time  output. 

3.3.2  Quanti zed-Epoch  Uncertainty 

A  step  towards  reducing  the  complexity  of  the  previous  detector  is 
to  quantize  the  epoch  uncertainty  region  of  Tc  seconds  into  a  (small)  finite 
number  of  alternatives,  i.e.,  points.  This  Is  equivalent  to  assuming  that  the 
epoch  rv  can  take  on  only  those  values  and  accordingly  develop  the  optimal  or 
near-optimal  detector  for  the  resulting  finite  hypothesis  problem.  The  quality 
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(as  well  as  the  complexity)  of  such  an  approximation  will  obviously  Increase 
^  with  the  number  of  points  considered.  It  will  soon  be  clear,  however,  that  the 

simplest  case  of  only  a  two-point  quantization,  which  will  be  the  one  we  analyze 
herein,  provides  quite  an  adequate  performance  with  respect  to  the  Ideal  (syn¬ 
chronous)  detector;  thus,  higher  level  quantization  and  complexity  might  even  be 

*  undesirable  from  a  practical  viewpoint.  In  addition,  a  suitable  Interpretation 
of  this  analysis  can  simultaneously  serve  as  a  bounding  technique  (upper  and 
lower)  for  the  optimal ,  continuous-epoch  uncertainty  detector  of  the  previous 
subsection,  whose  exact  analysis  Is  Intractable. 

•  Let  us  assume  that  the  two  points  comprising  the  epoch  uncertainty* 

are  eg  and  corresponding  to  the  offsets  e*0  and  e=  1/2  in  (3-1).  That 

is  to  say,  the  incoming  code  chips  are  likely  to  arrive  either  In  perfect  syn¬ 
chronization  (e  *  0)  or  half  a  chip  off  (e*  1/2)  with  respect  to  the  local  chip 
timing  reference,  and  they  are  both  equally  probable  to  occur.  A  variation  of 
the  optimal -likelihood  ratio  (3-28)  for  the  coherent  case,  adapted  for  the  above 
scenario,  would  Imply  that  the  decision  rule  should  be 


/JT  PTc+Tc/2 

rj(eo)  3  J  r(r)dx  ;  r  Ce  1  / 2^)  *  J  r(x)dx  (3-39) 

(M>TC  (MW? 

r 

Clearly,  the  chip-noncoherent  detector  would  form  the  envelopes  r j  In  place  of 
the  coherent  Integrations  (3-39);  otherwise,  the  development  Is  analogous. 

We  should  note  here  that  the  two  NTc-second  observation  Intervals  for 
the  two  corresponding  exponential  terms  of  (3-38)  are  slightly  off  by  half  a 
chip.  For  large  N,  such  "edge  effects"  are  insignificant  and  will  be  neglected 
in  the  following^. 

1  So,  at  this  point,  the  epoch  uncertainty  Is  assumed  to  be  discrete 

by  nature,  not  by  quantization. 

^For  that  purpose,  small  adjustments  will  be  made  for  convenience 
without  explicit  acknowledgement. 


The  familiar  small  SNR  approximation  applied  to  (3-38)  yields  the 
subopt imal  rule: 


jlj  rJ2(£0rj?1rj2(Cl/2)Hi  V 


(3-40) 


Let  r^  and  r^  Indicate  the  coherent  integrations  during  the  first  and  second 
halves  of  the  jth  chip  interval,  respectively,  i.e.. 


(j-l)T 


( j-l/2)T, 


r(t)dT  ;  rj2 


/jTc  r(t)dr 
(M/2)TC 


(3-41 


Then  (3-40)  can  be  rewritten  in  terms  of  (3-41)  as 


A. 1  * 


jii  (rjl  +  rJ2)2  +  jlj  (rj2  +  rj+l,l)2  £  A0 


(3-42 


or,  expanding  terms,  neglecting  edge  effects  (i.e.,  substituting  r^  1  with 
r^j],  and  dividing  by  1/2,  we  get 


I,  rjl2  *  rj22  *  rjZ  (rJl  *  i  Ao’ 


Figure  3.7  shows  a  possible  implementation  of  rule  (3-42). 

An  examination  of  the  decision  rule  (3-42b)  reveals  the  source  of 
inferiority  of  the  asynchronous  coherent  detector  in  comparison  with  its  syn¬ 
chronous  counterpart,  namely,  rule  (3-9).  The  latter  rule  can  be  obtained 
from  (3-42b)  if  the  factor  r.  .  ,  in  the  last  term  is  substituted  by  r..  be- 
cause  the  resulting  expression  is  then  the  perfect  square  (Viz)  ‘  V  of 
(3-9).  Since  r^  and  r^  ^  are  independent  rv's  (they  belong  to  different 
random  chips)  regardless  of  the  true  epoch  value,  it  follows  that  such  a  sub¬ 
stitution  would  Increase  the  mean  of  the  decision  statistic  A'  under  Hj  and 
slightly  reduce  the  variance— the  effect  of  which  is  a  net  improvement  in  per¬ 
formance  for  the  synchronous  detector.  We  will  now  take  a  closer  look  at  the 
performance  of  (3-42b). 


Without  any  loss  in  generality,  we  can  assume  that  the  true  epoch  is 
e  a  0.  The  symmetry  of  the  problem  suggests  that  the  results  to  be  derived  under 
such  conditioning  also  represent  the  unconditional  (average)  performance  in  this 
two-possible-value  case.  In  order  to  pursue  the  familiar  Gaussian  type  of  anal¬ 
ysis,  one  needs  the  means  and  variances  of  the  rv's  r^  and  rj2  as  per  (3-41). 
Those  are  readily  calculated  to  be 


(3-43) 


Furthermore,  it  is  easily  shown  that  (See  Appendix  E) 


#{rJ2  rj+l,l|Hk;£=S°}  “  0 


*{rjl  r j 2 ( Hk *e~°  } 


(3-44) 


for  k=0,l .  Combining  (3-42b),  (3-43)  and  (3-44),  it  follows  that 


,{ATkie"°}  ■  "(nr)  t1  s»] 


(3-45) 


A  complete  evaluation  of  system  performance  would  require  an  addi¬ 
tional  calculation  of  the  variance  var|A'|Hk|,  as  was  done  for  the  decision 
statistic  (3-9).  However,  matters  can  be  simplified  considerably  if  we  observe 
that  the  result  for  the  asynchronous  variance  should  be  higher  than,  but  close 
to,  the  synchronous  variance  of  section  3.1,  particularly  for  the  low  SNR  range 
of  interest.  The  main  rationale  supporting  such  an  assessment  is  the  compari¬ 
son  between  the  asynchronous  rule  (3-42  b)  and  the  expanded  synchronous  rule 
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which  shows  them  differing  in  only  one  of  four  contributing  terms.  Thus,  a 
variance  equal  to  N(NoTc)2/2  (see  (3-10b))  is  a  good  approximation  for  all 
cases. 

It  follows  from  the  above  discussion  and  the  underlying  Gaussian  as¬ 
sumption  that  the  primary  measure  of  performance  degradation  due  to  asynchronism 
is  the  reduction  in  the  difference  (distance)  between  the  means  of  the  two  dis¬ 
tributions  corresponding  to  the  alternative  hypotheses.  So,  if  we  rewrite  the 
synchronous  expression  of  (3-1  Oa)  as  A*  | ^  =  NCN0TcD[i  +  2  yc  5k{]/2  and 
compare  it  to  (3-45),  we  induce  a  loss  factor  of  4/3,  or  1.25  dB,  as  the  SNR 
penalty  associated  with  the  aforementioned  asynchronous  detector. 

The  above  estimated  loss  can  also  be  viewed  as  a  lower  bound  on  the 
performance  losses  of  the  continous-epoch  uncertainty  detector.  This  is  be¬ 
cause  the  two-point-uncertainty  coherent  detector  discussed  so  far  assumes  a 
random  epoch  that  can  take  on  only  two  values;  thus,  it  faces  less  uncertainty 
than  any  other  multipoint  or  continuous  type  of  epoch  randomness.  Therefore, 
it  is  bound  to  outperform  any  other  asynchronous  detect or,  including  the 
continuous-uncertainty  one. 

It  is  also  of  interest  here  to  derive  an  upper  bound  on  the  perfor¬ 
mance  losses  of  the  continuous-uncertainty  detector.  This  can  be  done  by  assum¬ 
ing  a  truly  continuous-epoch  uncertainty  and  viewing  the  two-point  detector  as 
a  suboptimal  quantized  implementation  of  the  optimal  continuous-uncertainty  one. 
The  performance  losses  of  the  suboptimal  detector  under  a  worst-case  choice  of 
the  unquantized  epoch  will  then  clearly  upperbound  the  losses  of  any  other, 
higher  complexity,  multipoint  or  continuous  detector. 

Some  reflection  could  persuade  us  that  the  worst  epoch  for  the  two- 
point  quantized  detector  is  the  one  straddling  between  e  *  0  and  e  ■  1/2  ,  namely, 
e  *  1/4  (or,  equivalently,  e  *  3/4,  whose  performance,  however,  is  identical). 
This  is  because,  under  Hj ,  rji  will  then  have  the  lowest  possible  mean  (namely, 
zero)  half  of  the  time  while  rj2  remains  unaffected.  Any  other  offset  will  yield 
a  higher  overall  average  for  rji,  thus  improving  performance. 

It  is  shown  in  Appendix  F  that,  under  the  worst  epoch  e  *  1/4,  the 
joint  statistics  of  rji  and  rj2  become 
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4j2rjl|Hk;eal/4} 

=  <f{rj2rj+l,l|Hkie’1/4 

which,  when  substituted  in 

(3-42b) ,  yields 

(3-46) 


(3-47) 


A  comparison  of  (3-47)  with  (3- 10a)  reveals  a  loss  factor  of  8/5,  or  2  dB,  with 
respect  to  the  synchronous  structure. 

In  conclusion,  the  lack  of  synchronism  in  the  coherent  detector  costs 
anywhere  from  1.25  dB  (lower  bound)  to  2  dB  (upper  bound).  Tighter  bounds  can 
be  obtained  by  reiterating  the  above  arguments  for  epoch-quantized  detectors 
with  more  than  two  quantization  points.  As  the  number  of  points  increases,  one 
asymptotically  converges  to  the  actual  performance  of  the  continuous-epoch  de¬ 
tector.  In  any  case,  a  rough  figure  of  an  average  1.5-dB  loss  due  to  asynchron- 
ism  should  not  be  far  from  exact.  Similar  conclusions  can  be  drawn  for  the 
noncoherent  detectors  of  the  previous  sections. 
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3.4  Energy  Detector  (Radiometer) 


* 


The  most  commonly  used  energy  measuring  device,  the  "energy  detector11 
or  "radiometer",  will  now  be  analyzed  as  a  potential  candidate  LPI  detector. 

In  the  battle  for  implementational  savings  and  hardware  reduction,  the  radiom¬ 
eter  is  an  easy  winner:  it  is  simply  the  cascade  of  a  wide  bandpass  filter, 
followed  by  a  square-law  device,  a  zonal  LPF  (possibly  excluded),  an  integrator* 
and,  finally,  a  threshold  comparator  which  produces  the  decision,  as  shown  in 
Figure  3.8.  The  bandwidth  of  the  BP  filter  should  approximately  cover  the  null- 
to-null  bandwidth  of  the  DS  spread  signal  to  be  detected;  thus,  approximate 
knowledge  of  the  chip  rate  is  an  implicit  necessity. 

Clearly,  the  hypothesis  discriminating  power  of  the  radiometer  is 
based  on  simple  energy  considerations:  signal  plus  noise  funder  should 
possess  more  energy  than  noise  alone  (under  Hq}  •  Thus,  by  properly  adjusting 
the  decision  threshold,  a  correct  decision  will,  hopefully,  be  made.  It  goes 
without  saying  that  performance,  as  measured  by  the  pair  (po»pp^)»  increases 
monotonically  with  increasing  SNR  y  and  observation  time  T=NT  .  On  the  other 

w  C 

hand,  selection  of  the  RF  bandwidth  Wgp  in  the  filter  is  subject  to  optimiza¬ 
tion  and  is  discussed  more  thoroughly  in  the  following. 

It  has  long  been  known  that  the  precise  analysis  of  systems  contain¬ 
ing  linear  plus  nonlinear  elements,  as  in  Figure  3.8  for  arbitrary  filter  shapes 
and  bandwidths,  is  generally  unwieldy  [11].  If  the  envelope  detector  of  Figure 
3.8  is  assimilated  to  (or  approximated  by)  an  inphase-and-quadrature  (I-Q)  type 
of  envelope  detector,  it  can  then  be  claimed  that  the  output  statistic  A  has  a 
chi-squared  distribution  (central  or  noncentral  ,  depending  on  whether  a  signal 
is  absent  or  present,  respectively) ,  with  2NWgpTc  degrees  of  freedom.  The  prob¬ 
abilities  Pp^  and  Pg  can  then  be  calculated  by  a  numerical  integration  of  the 
corresponding  densities,  from  the  threshold  Aq  to  infinity,  since  those  inte¬ 
grals  cannot  be  found  in  closed  form.  This  procedure  can  be  very  laborious. 
Fortunately,  when  the  time-bandwidth  product  Wgp(NTc)  is  large  (which  is  the 
case  here,  since  N  is  large),  one  can  very  accurately  approximate  the  proba¬ 
bility  densities  of  A  under  both  hypotheses  by  Gaussian  ones,  by  virtue  of  a 
central -1 imit-type  argument.  Then  calculation  of  the  overall  performance  be¬ 
comes  fairly  straightforward. 
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Since  the  integrator  itself  acts  as  a  lowpass  filter,  the  zonal  LPF 
preceding  it  is  typically  eliminated;  it  is  included  here  as  a  conceptual  aid. 
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The  received  waveform  r(t)  is  described  by  (3-2).  The  BPF  is  assumed 
to  be  rectangular  with  bandwidth  WBp.  The  output  of  the  BPF  in  Figure  3.8  is 

rBp(t)  *  /2Tc(t)  cosQiigt  +  <►)  +  nBP ( t )  (3-48) 

where  nBp(t)  is  filtered  AGN  with  the  typical  BP  representation 

nBp(t)  »  /^(t)  cos  wQt  +  nQ(t)  sin  uQtJ  (3-49) 

In  (3-49),  nj ( t)  and  r>q( t)  are  approximately  Independent,  zero-mean,  jointly 
Gaussian,  lowpass  noise  processes,  each  with  a  PSD  of  Ng  W/Hz  (one-sided)  and 

bandwidth  of  WBp/2  Hz  (one-sided).  The  code  signal,  filtered  by  the  BPF,  is 

designated  c(t)  in  (3-48).  The  zonal  LPF  (shown  for  analytical  purposes  only) 
rejects  the  2fQ  components  of  the  square-law  output  z(t)  and  has  an  output 
given  by 

zLp(t)  =  s(c(t))2  +  nj2(t)  +  np2(t)  +  2vTc(t)  n^t)  (3-50) 

In  order  to  maintain  analytical  tractabil Ity,  we  approximate  and  model 
c(t)  as  an  attenuated,  but  undlstorted,  version  of  c(t),  namely, 

c(t)  -  ac(t)  (3-51) 

2 

where  the  "attenuation  factor11  a  measures  the  (normalized)  power  reduction  due 
to  filtering: 


a2  » 


(3.52) 


2 

since,  for  a  ±l-valued  code,  the  denominator  of  (3-52)  equals  c  (t)  *  1. 

In  the  previous  expression,  Sc(f)  is  the  PSD  of  the  random  code  signal  c(t) 

and  H  (f)  is  the  lowpass  equivalent  of  the  BPF  HDO(t),  i.e., 

L  or 


C 


HL(f) 


i  m<wBP/2 

0  ,  otherwise 


By  suitably  adjusting  H^(f),  (3-52)  can  be  generalized  to  filters  other  than  the 
perfectly  rectangular  one  considered  above. 

From  the  above  definitions.  It  follows  that  CS(t)}2  *  a2.  Thus, 


zLp(t)  *  a2S  +  nj2(t)  +  nQ(t)  +  2^ac(t)  nj(t) 


(3-53) 


under  Hj.  Under  the  alternative  HQ,  zLp(t)  Is  given  by  the  above  expression 
with  S  =  0.  It  follows  that 

*|zLP(t)lHki  *  ^  6kl  +  N0  WBP;  k*°»1  (3-54) 

with  skl  as  per  (3-10c).  Thus,  the  expected  value  A | |  of  test  statistic  A 
at  the  output  of  the  Integrator  in  Figure  3.8  can  be  calculated  under  either 
hypothesis  Hk*,k»0,l,  as 

'HM  *  J*  '{ZLP(t)|  Mdt  ■  NTc[n0Wbp  ♦  a2S«  J;  k-0,1  (3-S5) 

0 

The  next  step  Is  to  obtain  the  second-order  statistics  of  z^p(t)  and 
A.  The  typical  approach,  based  on  the  assumption  WBpNTc  »  1,  Is  to  model 
zLp(t)  as  a  very  wideband  ("delta-correlated")  process  with  respect  to  the  ap¬ 
proximate  bandwidth  ([NT^"1  °f  the  Integrator.  Thus,  one  needs  to  evaluate  only 
the  PSD  of  z|_p(t)  at  the  origin  f  ■  0,  distinguishing  between  the  average  (con¬ 
stant)  contribution  and  the  random  contribution.  Using  the  fact  that  n j ( t )  and 
ng(t)  are  Gaussian,  the  autocorrelation  function  r *  z^p(t)  can  be  obtained  di¬ 
rectly;  from  that,  its  PSD  under  Hj  Is  given  by 

2 

SzLP(f)  -  (a2S  +  N0WBp)  «(f)  +  4Sni(f)  ®  Snp(f)  +  4a2  SSc(f)  ®  Snj(f)  (3-56) 

where  ©means  "convolved  with."  The  coefficient  of  the  Dirac  delta  function 
6(f)  agrees  with  the  first-order  statistics  of  zLp(t)  obtained  directly  In 
(3-54).  The  two-sided  PSD  of  z^p(t)  At  f  *  0  measures  its  random  contribution 


c 


i 

! 


<4 


a 


and  is  obtained  by  evaluating  the  last  two  terms  of  Sz^pff)  in  (3-56)  with  f*0, 
resulting  in 


S,  (01 


Nq2  WBp  +  2a4  SNq  ^  ;  k-0,1 


(3-57) 


The  delta-correlated  random  part  of  zLp(t)  contributes  to  the  variance  of  a  as 


•f|Hk}  ■  KK'0’  ■  (ntc)(mo2“bp  *  2»4  SN0  *kl)  <3-5e> 


The  mean  and  variance  of  A  in  (3-55)  and  (3-58)  suffice  to  characterize  perfor¬ 
mance  under  the  Gaussian  assumption  which  hinges  on  the  fact  that  wNpNTc"  N»  1. 
The  result  is  (3-11)  with 


;  b  *  2a  ‘ 


BP  'c 


(3-59) 


Clearly,  the  choice  of  Wgp  has  an  Impact  on  a  through  (3-59)  and  a  in  (3-52). 
The  quantity 


(I)  (¥)*<* 


(3-60) 


BP  c 


BP  'c 


has  been  plotted  in  Figure  3.9  as  a  function  of  the  product  wbpV  The  max¬ 
imum  a  •  0.77  is  attained  at  W00T  *  1.  This  corresponds  to  a  minimum  SNR 
max  or  c 

loss  of  1.1  dB  with  respect  to  the  chip  noncoherent  detector  (3-24).  Since 
the  above  conclusions  are  based  strictly  on  energy  considerations  and  the  Inter¬ 
chip  interference  effects  due  to  filtering  have  not  been  taken  Into  account, 
the  actual  losses  are  somewhat  higher.  In  the  following  comparisons,  we  select 
a  nominal  loss  figure  of  -1.5  dB.  We  note  that  the  above  losses  pertain  to  a 


wS 


rectangular  BPF;  other  filter  shapes  will  result  In  different  loss  figures, 
typically  of  the  same  order.  For  example,  [10]  can  be  consulted  for  the  effect 
of  filtering  on  a  baseband  rectangular  pulse  In  terms  of  output  SNR  losses, 
compared  to  a  matched  filter. 


3.5 


Comparisons  and  Discussion 

All  the  previous  theory  was  developed  under  the  cardinal  assumption 
that  the  received  waveform  r(t),  on  which  the  decision  is  based.  Is  of  the  form 
(3-2).  Excluding  the  fact  that  the  spreading  Is  done  via  DS  methods  (alterna¬ 
tives  such  as  frequency  and  time  hopping  will  be  studied  In  the  sequel),  three 
key  restrictions  In  the  system  modeling  were  applied:  (a)  under  Hj,  the  whole 
observation  Interval  (0,T)  is  occupied  by  the  spread  waveform,  (b)  there  Is  on¬ 
ly  one  version  of  signal  plus  noise  available,  and  (c)  the  additive  disturbance 
is  white  Gaussian  noise.  Let  us  now  comment  on  each  of  these  restrictions. 

Assumption  (a)  Is  not  as  restrictive  as  it  might  seem  at  first.  Ad¬ 
mittedly,  if  the  starting  signal  transmission  time  Is  unknown,  the  intercepting 
receiver  will  be  forced  to  look  at  T-second  observation  Intervals  that  are  ran¬ 
domly  placed  with  respect  to  the  actual  transmission  time.  One  could  develop 
a  theory  along  lines  similar  to  those  presented  herein,  whereby  the  performance 
losses  due  to  such  realistic  randomness  are  taken  Into  account.  Although  the 
performance  of  each  individual  suggested  system  will  certainly  be  altered,  their 
relative  ranking  in  terms  of  effectiveness  Is  expected  to  remain  the  same.  Thus, 
the  conclusions  and  comparisons  which  we  shall  soon  undertake  should  carry  over 
more  or  less  Intact.  In  order  to  be  definite,  however,  one  should  take  a  closer 
look  In  the  future*. 

Assumption  (b)  is  more  substantial  in  nature.  In  fact.  It  will  now 
be  shown  that  performance  Improvement  can  be  gained  by  utilizing  two  receptions 
which  cover  the  same  geographical  area  (hence,  contain  the  same  useful  signal), 
but  possess  independent  RF  circuits  (hence.  Independent  thermal  noise  compon¬ 
ents).  In  order  to  briefly  Illustrate  the  Idea,  let  us  consider  a  synchronous 
coherent  system. 

Let 


rj(t)  a  ^  c(t)  6kl  +  nj(t)  /o  <  t  <  T\  (3_61) 

r2(t)  =  ^  c(t)  5kl  +  n2(t)  '  k*0,1  ' 


represent  the  two  receptions.  Here  the  two  noise  processes  nm(t);m*l,2  are 
Gaussian  and  Independent,  while  the  signal  component  c(t)  Is  the  same  in  both 


One  could  create  a  multitude  of  possible  solutions,  such  as  over¬ 
lapping  versus  nonoverlapping  observation  intervals,  etc. 


waveforms.  We  should  point  out  that,  when  Implementing  the  cross-correlation 
which  we  are  about  to  discuss,  maximum  care  must  be  exercised  to  ensure  that 
the  two  code  signal  components  are  indeed  in  phase  since  timing  offsets  of  more 
than  a  code  chip  will  result  in  involuntary  signal  spreading  and  total  loss  of 
detection  capability. 

The  decision  rule  adopted  here  is  a  simple  correlation  (see  Figure  3.10) 


/NT  N  /, \  m  "1 

C  r.(t)  r2(t)dt  =  1  r.(1)  r.(2)  l  *' 

1  1  j-1  J  J  hq 


(3-62a) 


where 


•  r 


rjt)  dt 


(j-D 


(3-62b) 


We  note  that  (3-62)  is  analogous  to  (3-9)  once  r?  is  substituted  by 
Clearly,  in  the  absence  of  thermal  noise,  both  (3-9)  and  (3-62)  produce  the  same 
quantity,  NSTC  ;  however,  performance  is  different  In  noise.  Indeed,  from  (3-26), 
it  easily  follows  that  (see  Appendix  G): 


rjA'|Hk[  *  n(n0tc)  yc  5kl 


(3-63a) 


HATki  -  N(¥c)2  H +  \  Q 


(3-63b) 


The  resulting  performance  is  again  described  by  (3-11),  where  now 


a  =  2  ;  b  =  4 


(Two  Independent  receptions) 


(3-64) 


Thus,  a  comparison  with  the  corresponding  performance  parameters  (3-12a)  for 
the  (one- reception)  synchronous  coherent  radiometer  reveals  a  gain  factor  of 
or  1.5  dB  In  SNR  for  the  present  system.  However,  we  would  again  like  to 
emphasize  that  this  gain  is  attained  at  the  cost  of  higher  complexity. 


Before  commenting  on  the  third  assumption  (c),  let  us  attempt  an 

evaluation  of  the  relative  merits  of  all  the  previous  detector  structures. 

The  conclusions  have  been  summarized  in  Figure  3.11  for  the  exemplary  case  of 
_2 

PpA  3  10  and  N  =  1000  chips.  Clearly,  the  ultimate  upper  bound  In  performance 
corresponds  to  the  ideal  (but  unrealistic,  from  the  LPI  viewpoint)  case  of  a 
totally  known  signal  ([under  HQ,  l.e.,  a  code-matched  receiver.  Since  It  Is  an 
unachievable  upper  bound,  it  has  been  Included  merely  to  serve  as  a  comparative 
yardstick.  It  is  represented  by  curve  (D  ,  whose  analytical  description  is  via 
the  distance  dUB  (as  per  (2-2)3: 

dUB  =  0-65) 

where  a  coherent  system  has  been  assumed.  Comparing  (3-65)  with  (3-12b)  implies 
that  there  is  a  loss  factor  of  associated  with  the  lack  of  code  knowledge, 
and  this  factor  is  significant  at  low  SNR's.  On  the  other  hand,  all  performances 
are  lower  bounded  by  that  of  the  radiometer,  which  was  derived  in  the  previous 
subsection  (see  (3-59)3  and  is  represented  by  curve  ©  in  Figure  3.11.  The  re¬ 
maining  curves  0  ,  (3)  and®  pertain  to  alternate  receivers  which  were  consid¬ 
ered  herein:  ©  corresponds  to  a  detector  utilizing  two  independent,  synchron¬ 
ous,  coherent  receptions,  0  to  one  synchronous,  coherent  detector,  and  ®  to 
a  synchronous,  chip-noncoherent  detector.  The  difference  between  0  and  0  is 
roughly  1.5  dB  which  is  also  the  difference  between  0  and  ®  ,  as  well  as  ® 
and  ©.  It  follows  that,  under  fairly  ideal  conditions,  i.e.,  two  Independent 
synchronous  coherent  receptions,  one  could  utilize  up  to  4.5  or  5  dB  gain  above 
the  radiometer;  each  of  the  added  features  (i.e.,  two  receptions,  synchronism 
and  coherency)  can  be  thought  of  as  contributing  1.5  dB  to  the  gain.  We  note 
that,  in  arriving  at  these  gain  figures,  system  Imperfections  and  noncal 1b rated 
parameters  such  as  doppler  frequency  offset  have  not  been  accounted  for. 

In  conclusion,  it  is  seen  that  certain  gains  are  plausible  with  re¬ 
spect  to  the  radiometer  if  careful  designs  are  employed,  but  they  are  certainly 
not  overwhelming— at  least  in  the  SNR  range  of  interest.  This  is  because  a 
high-rate  DS  waveform  is  not  very  distinguishable  from  the  background  thermal 
AWG  noise  In  which  it  is  detected  and,  as  Is  well  known,  the  radiometer  then 
becomes  an  asymptotically  optimal  detector.  This  is  also  in  accordance  with 
Krasner’s  approximate  result  regarding  biorthogonal  waveforms  ([6,  eq.  (29)3- 


1. 

2. 

3. 
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Completely  known  waveform 

Two  independent  receptions,  synchronous  coherent  detector 
Synchronous  coherent  energy  detector 
Synchronous  chip-noncoherent  energy  detector 
Radiometer 
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Figure  3.11.  Performance  Comparisons  for  the  Detection  of  DS  Waveforms 
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However,  if  one  removes  the  biorthogonal Ity  of  the  signal  set  and/or 
the  AWG  nature  of  the  noise,  the  results  can  differ  substantially.  One  illus¬ 
tration  of  the  former  case  can  be  found  in  the  following  section  regarding  time 
hopping.  Furthermore,  if  one  removes  assumption  (c)  about  the  white  Gaussian 
nature  of  the  noise,  i.e.,  if  some  colored  interference  is  superimposed  on  the 
thermal  background  noise,  a  whole  new  class  of  (mostly)  open  problems  is  gener¬ 
ated.  Tackling  such  problems  from  a  theoretically  optimal  standpoint  would  be 
a  much  more  complicated  task  in  general  than  that  which  has  been  performed  up 
to  now  for  AWGN;  nonetheless,  the  benefits  derived  from  applying  certain  sub- 
optimal  solutions  (which  are,  however,  more  sophisticated  than  the  radiom¬ 
eter)  could  be  signficant.  One  particular  case  of  detecting  DS  waveforms  in 
random-tone  i nterference— and  in  the  absence  of  thermal  noise— is  treated  in 
section  5.2,  wherein*  it  is  shown  that  the  application  of  real-time  auto¬ 
correlation  techniques  can  bring  about  impressive  gains  with  respect  to  the 
radiometer.  Although  this  is  a  rather  ad-hoc  technique  and  the  hypothesized 
scenario  is  not  the  most  general  one,  it  still  serves  as  a  good  motivating 
example  of  how  benefical  an  extra  degree  of  sophistication  can  be  when  the 
problem  formulation  and  conditions  vary. 


Although  the  development  there  formally  pertains  to  a  FH/DS  hybrid 
system,  it  essentially  addresses  DS  waveforms  since  it  is  concerned  with  the 
signal  processing  on  a  per-band  (single  hopping  slot)  basis. 
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4.0  DETECTION  OF  TIME-HOPPING  WAVEFORMS 

We  now  focus  on  the  wideband  detection  of  TH  waveforms,  where  the  code 
dictating  the  hopping  is  again  assumed  random.  We  shall  consider  only  synchronous 
detectors  (i.e.,  the  timing  or  epoch  of  the  hopping  slot  will  be  assumed  known); 
asynchronism  can  be  treated  with  methods  similar  to  those  used  in  section  3.3. 
Furthermore,  both  carrier-coherent  and  carrier-noncoherent  systems  will  be  exam¬ 
ined.  In  particular,  the  former  will  assume  a  known  carrier  phase  combined  with 
a  pulse-position-modulation  (PPM)  format  while,  for  the  latter,  any  modulation 
that  randomizes  the  phase  from  frame  to  frame  (e.g.,  BPSK)  is  well  suited.  We 
term  such  cases  "pulse  noncoherent."  Clearly,  one  could  hypothesize  a  nonco¬ 
herent  PPM  system  with  an  unknown  constant  carrier  phase  and  derive  the  optimal 
detection  rule.  However,  the  resulting  receiver  and  pertinent  analysis  are  much 
too  complicated  (for  a  similar  situation,  refer  to  section  3.2  for  DS);  there¬ 
fore,  only  pulse- noncoherent  systems  will  be  discussed  here.  Again,  the  lower- 
bounding  performance  of  the  radiometer  will  indicate  that  very  little  is  lost 
by  using  this  simplification. 

For  our  purposes,  the  received  waveforms  (TH  signal  plus  noise)  can 
be  written  as 

00 

r(t)  =  SZ5  I  p(t  -  kTF  -  pkTH)cos(u>0t  +  ek)  +  n(t)  (4-1) 
k — 00 

where  Tp  is  the  frame  length,  T^  is  the  hop  length  (width  of  each  time  slot),  and 
pk’0k  are  rand°m  variables  that  are  Independent  from  frame  to  frame*  and  of  each 
other,  denoting  the  slot  location  and  carrier  phase,  respectively,  during  the  Jcth 
frame.  Here,  pk  can  take  on  any  one  of  the  equiprobable  values  pk=0,.. .,Np-l, 
where  Np  -  Tp/T^  is  the  total  number  of  slots  per  frame  and  ek  summarizes  both 
the  unmodulated  carrier  phase  $  plus  any  superimposed  PSK  modulation.  For  coher¬ 
ent  systems,  ek  is  assumed  known  for  every  k.  Finally,  p(t)  is  a  unit  pulse  of 
duration  T^  seconds,  while  n(t)  is  the  usual  bandpass  AWGN.  A  noiseless  sample 
waveform  (realization)  of  the  T^  signal  is  shown  in  Figure  4.1.  We  shall  assume 
that  the  total  observation  interval  consists  of  Q  frames,  i.e.,  T  »  QTp  =  QNpTH. 


The  independence  of  is  within  the  aforementioned  spirit  of  a 
"pul se-noncoherent"  system. 
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4.1  Synchronous  Coherent  Detectors 

Under  the  coherent  assumption,  "1"  is  substituted  for  /l  cos(a>Qt+8|<) 
in  (4-1).  Then  a  straightforward  application  of  the  generalized  likelihoood 
ratio  yields 


A0*(t)3  -  {  exp( 

m=l  )  nF  k=l  ( 


(4-2) 


where 


4  STH 
'  ^ 


(4-3) 


is  the  SNR  per  time  hop  (slot)  and 


/(m-l)Tp+kTH 

M  r(t)  p(t  -  kTH}dt 

(m-l)TF+(k-l)TH 


(4-4) 


is  the  integral  of  r(t)  in  the  Jcth  slot  of  the  mth  frame.  If  yh  is  a  small  nun- 
ber,  so  that  exp{x}  *  1+x,  (4-2)  can  be  simplified  to  the  approximate 


A*  Cr(t)3 


f 


H>i 

r(t)<lt  „<  A'0 


(4-5) 


i.e.,  just  the  coherent  integral  of  r(t).  However,  since  y^  is  typically  above 
0  dB  (much  larger  than  y^,  such  approximations  are  not  as  successful  as  in  the 
DS  case— a  fact  to  which  the  comparisons  in  section  4.3  will  attest. 

The  distance  dc0^  for  the  coherent  rule  (4-2)  can  be  derived  based  on 
(2-3)  and  the  following  steps,  as  in  [3]: 


coh 


£n  <?{a2|hq}  =  Q  £n  <*{af2|h0} 


(4-6) 


where  Ap  is  the  likelihood  ratio  per  frame.  But, 
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(4-7) 


since  rkm  is  a  zero-mean  Gaussian  random  variable  (rv)  under  Hq.  Combining 
(4-6)  and  (4-7)  results  in 


dcoh  =  J  +  ^  (*?"  -  l)_ 


(4-8) 


The  distance  d  . 

coh ,appr 

since 


for  the  approximate  rule  (4-5)  is  derived  much  easier 


(4-9a) 


(4-9b) 


(4-10) 


It  is  clear  from  (4-8)  and  (4-10)  that  d^  asymptotically  approaches  dCQ^  ap^r 
as  +  0,  an  expected  result.  Furthermore,  we  can  compare  the  approximate 
coherent  rules  (3-9)  for  DS  versus  (4-5)  for  TH  and  note  that  the  first  Involves 
a  nonlinear  operation  on  the  data  prior  to  integration,  while  the  second  does 
not.  This  has  a-reflection  on  their  corresponding  distances  as  (3-12b)  is  di¬ 
rectly  proportional  to  yc»  while  (4-10)  Is  proportional  to  Thus,  for  the 

very  low  SNR  case,  the  latter  would  outperform  the  former;  therefore,  TH  is  a 
more  detectable  waveform  than  DS. 


4.2  Synchronous-Noncoherent  Detectors 

Starting  from  (4-1)  and  Invoking  the  independent  phase  assumption 
results  in  the  decision  rule 


(4-11) 


where  is  the  Jcth  slot,  mth  frame  envelope 


and 


(4-1 2a) 


eI,km 

®Q,km 


/(m-l)T.+kTu 

F  H  r(t) 

(m-l)TF+(k-l)TH 
k-l,...,Np  ;  m=l,... 


(4-1 2b) 


The  distance  dnonco^  associated  with  (4-11)  is  derived  by  steps  identical  to 
(4-7)  as 


dnoncoh  *  f  *»[‘  »  J^'oC^  '  l)]  <4'13> 

Since  IQ(x)  <  exp{x}  for  every  x  >  0,  it  follows  from  the  comparison  of  (4-8) 
and  (4-13)  that  dcoh  >  dnQncoh  for  every  yh,  as  expected. 

Instead  of  analyzing  the  noncoherent  rule  which  approximates  (4-li) 
for  low  SNR,  let  us  consider  only  the  performance  of  the  radiometer  that  oper¬ 
ates  on  the  time-hopped  waveform  (4-1)  for  seconds.  Following  the  steps 
outlined  in  section  3.4,  it  can  be  shown  that 

drad  *  ^Krad(l^)YH  (4'14) 


.  ^  a4/Wgp  Th  again  measures  the  loss  due  to  filtering. 


where  the  constant  K, 


As  in  section  3.4,  Krad  can  be  set  at  its  optimal  value  Kra(j  ■  (0.77 )z  •  0.6. 

We  should  note  here  that  (4-14)  Is  actually  an  optimistic  prediction  (i.e., 
an  upper  bound)  of  the  radiometer's  performance  since  It  is  based  on  the  assump¬ 
tion  that  the  variance  of  the  test  statistic  under  Hj  Is  approximately  the  same 
as  under  Hq.  For  medium-to-hlgh  yh,  this  Is  not  true  (unlike  the  DS  case  where 
it  is  justified  by  the  low  values  of  YcD*  Although  a  more  meaningful  perfor¬ 
mance  description  can  easily  be  derived*,  we  shall  be  content  with  (4-14)  for 
comparison  purposes. 


It  will  consist  of  a  pair  of  equations  similar  to  (3-11)  and  (3-59) 
under  an  appropriate  interpretation  of  the  parameters  involved. 
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4.3  Comparisons  and  Discussion 

First,  we  are  Interested  In  comparing  the  relative  loss  of  the  approx¬ 
imate  coherent  rule  (4-5)  versus  the  exact  (4-2)  since  the  former  can  be  imple¬ 
mented  trivially  (a  simple  Integrator)  In  contrast  to  the  significant  complex¬ 
ity  of  the  latter.  The  comparison  Is  In  terms  of  the  relative  SNR  values  Y^ 

(in  dB)  required  by  each  in  order  to  achieve  the  same  performance  level  (d^) 
for  the  same  fixed  values  of  Q  and  Np.  It  is  based  on  (4-8)  and  (4-10)  and  is 
shown  in  Figure  4.2  (Np  *  10).  As  expected,  the  difference  diminishes  at  low 
SNR,  but  is  rather  pronounced  at  medium- to -high  SNR.  So,Y^oh  »  0  dB  corre¬ 
sponds  toY^oh,appr  =  3.9  dB— a  gap  that  Increases  rapidly  as  Y^oh  Increases. 

Second,  we  look  at  the  gains  brought  about  by  the  coherent  assumption 
versus  the  noncoherent  by  comparing  the  "exact"  distances  In  (4-8)  and  (4-13). 

The  result,  again  In  terms  of  the  required  SNR's  for  the  same  performance  level, 
is  shown  in  Figure  4.3.  The  difference  here  diminishes  as  SNR  Increases,  while 
it  can  be  substantial  at  low  SNR.  The  final  comparison  Is  between  the  noncoher¬ 
ent  "exact"  performance  (4-13)  (or  Its  equivalent  log-likelihood  sum)  versus 
the  radiometer  performance  estimate  (4-14).  It  Is  shown  in  Figure  4.4  for 
Krad  3  0,6  and  NF  “  10-  The  irreduc^ble  distance  of  1.1  dB  as  the  SNR  goes  to 
zero  is  due  to  the  ^Kra(J  factor*.  Again,  we  should  keep  In  mind  that  the  actual 
SNR  losses  of  the  radiometer  are  higher  than  those  shown  In  Figure  4.4  by  an 
amount  that  Increases  with  SNR  due  to  previously  discussed  reasons.  Furthermore, 
a  common  trend  is  evident  from  Figures  4.2  and  4.4,  namely,  that  the  simplify¬ 
ing  deviations  from  the  optimal  decision  rules  Incur  comparative  losses  that 
increase  fairly  rapidly  with  the  available  hop  SNR  for  values  of  Y^  above 
2  -  3  dB.  Since  this  Is  the  dominant  range  of  importance  In  TH  applications, 
optimal  devices  (albeit  complex)  should  attract  due  attention. 


*However,  this  margin  will  probably  disappear  if  the  loss  due  to 
asynchronlsm  Is  accounted  for  In  an  asynchronous,  noncoherent,  optimal  detector. 
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Figure  4.4.  Comparison  of  the  Noncoherent  Rule  (4-11)  versus  the  Radiometer 
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5.0  DETECTION  OF  FREQUENCY-HOPPING  SIGNALS 

USING  AUTOCORRELATION  TECHNIQUES 

We  shall  now  focus  on  the  detection  of  FH  or  FH/DS  waveforms  embedded 
in  different  forms  of  additive  Interference.  The  two  particular  cases  of  inter¬ 
est  are:  (1)  FH  signals  in  AWGN  (section  5.1)  and,  (2)  FH/DS  signals  In  random- 
tone  interference  (section  5.2).  Although  optimal  detectors  can  be  derived— at 
least  for  case  (1)— following  the  guidelines  of  Sections  3  and  4,  those  will  not 
receive  any  in-depth  treatment  here  because  their  exact  implementation  Is  usually 
hampered  by  excessive  complexity.  Instead,  we  shall  deal  with  Improved  suboptlmal 
wideband  detection  procedures,  such  as  real-time  autocorrelation-domain  algorithms 
Our  approach  has  been  motivated  by  recent  advances  in  the  technology  of  real-time 
autocorrelation  and  convolution  devices  with  large  time-bandwidth  products  [12, 
13].  The  rationale  for  introducing  such  techniques  in  the  context  of  wideband 
detection  is  explained  herein. 

Let  s(t)  be  the  frequency-hopped  signal  to  be  detected  [under  Hj}  In 
AWGN  within  a  total  spread  bandwidth  of  Ws  Hz.  If  RH  ^  T^'1  is  the  hopping  rate 
and  the  hopping  frequencies  are  contiguous  and  equispaced  by  the  minimum  nonco¬ 
herent  orthogonal  separation  Af  *  R^,  It  follows  that  the  total  number  M  of  pos¬ 
sible  hopping  locations  is  M  *  W$/RH  *  W$TH— typically,  a  very  large  number.  If 
the  total  observation  time  T  is  for  Ny  hops,  l.e. ,  T  =  NyTH,  **  can  then  be  shown 
under  suitable  assumptions  that  the  log  likelihood  ratio  reduces  to 


where  Rnm  is  the  envelope  output,  at  the  end  of  the  rith  hopping  Interval,  of 
the  matched  filter  centered  around  the  mth  hopping  frequency  (spectral  slot), 
I0(*)  is  the  zeroth-order  Bessel  function,  and  S,Ng  are  the  signal  power  and 
one-sided  noise  power  spectral  density  (PSD),  respectively.  Equation  (5-1) 

Is  the  dual  to  (4-11)  for  the  TH  case;  it  thus  follows  that  performance  Is 
quantified  by  the  "exact"  distance 


(5-2) 
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where  yh  is  the  hop  SNR 


(5-3) 


analogous  to  (4-13). 

The  most  inconvenient  aspect  of  the  optimal  rule  (5-1)  is  perhaps 
its  complexity:  whether  implemented  as  a  bank  of  matched  channels  or  via  an 
FFT  (which  would  simultaneously  provide  all  the  spectral  estimates  Rnm,ma>l, 

. ..,M  for  each  hopping  time]),  the  resulting  complexity,  measured  by  the  number 
of  required  filters,  can  be  enormous.  One  way  of  alleviating  the  problem  is 
shown  in  Figure  5.1:  the  total  bandwidth  W$  is  subdivided  into  large  contig¬ 
uous  segments  of  B  Hz,  each  being  much  larger  than  the  optimal  bandwidth  of 
Rh  Hz.  Appropriate  processing  produces  a  per-band  decision  which  is  then  fed 
into  an  overall  accumulator  (e.g.,  a  majority-logic  combiner)  for  a  final  de¬ 
cision  (Hj  versus  Hq).  Let  us  note  that  alternate  reduced  configurations  have 
also  appeared  in  the  literature  [4,5],  such  as  the  partial-band  filter-bank 
combiner.  There,  only  a  fraction  of  the  total  number  M  of  the  hopping  slots 
is  being  observed,  but  the  observation  bandwidth  per  slot  is  optimal  (B  »  Rh}. 

A  comparison  between  alternatives  has  not  yet  been  performed. 

The  focus  of  this  section  Is  the  waveform-processing  (WP)  aspect  on 
a  per-band  basis.  Clearly,  algorithms  Improving  performance  on  that  level  will 
also  increase  the  overall  system  performance.  We  shall  assume  throughout  this 
development  that  the  time-bandwidth  product 

G  ®  BTH  *  B/Rh  »  1  (5-4) 


is  very  large,  i.e. ,  of  the  order  of  hundreds.  The  factor  G  also  measures  the 
order  of  the  reduction  in  complexity  when  Implementing  the  suboptimal  structure 
of  Figure  5.1  Instead  of  the  optimal.  As  a  consequence,  the  Input  SNR  per  band, 
Y.jn,  defined  as 


(5-5) 


is  typically  very  small  compared  to  unity  (]yfn  «  l}. 


y. 


Two  questions  arise  immediately: 

(1)  Into  how  many  bands  should  the  spectrum  be  divided? 

(2)  What  kinds  of  WP  alternatives  are  available  or  could  be  used? 

With  respect  to  (1),  It  is  intuitively  clear  that  one  should  use  as 
many  bands  as  the  number  of  devices  which  can  be  aff on  id  since,  by  decreasing 
B  and  thus  increasing  y^n  from  (5-5),  more  reliable  decisions  on  a  per-band  ba¬ 
sis  can  be  achieved;  an  elaborate  cost-versus-payoff  study  is  of  significant 
interest  in  this  area,  but  will  not  be  pursued  here.  As  for  the  second  question, 
we  note  that  the  device  which  immediately  comes  to  mind  is  the  familiar  radiom¬ 
eter,  designed  to  measure  energy  in  the  band  of  B  Hz,  as  shown  in  Figure  5.2. 
Under  the  Gaussian  assumption  for  the  output  decision  statistic  Y  (which  is 
well  justified  for  very  large  G  via  a  central -limit-type  argument),  its  per¬ 
formance  is  easily  derived  to  be 


drad  '  <r%]  *  O'^d]  * 


(5-6a) 


where 


'  **2?  * 


r- 1  2 

G  -yh 


(5-6b) 


Note  that  the  second  equality  in  (5-6b)  signifies  the  small  signal  suppression 
effect  of  square-law  detectors  [ll,  page  267].  Again,  it  should  be  emphasized 
that  [P|)>Pp^3  P®rtain  to  the  per-band  decision— not  the  overall  decision  scheme. 

The  autocorrelation-domain  methods  to  be  elaborated  upon  in  the  follow¬ 
ing  sections  can  be  viewed  as  a  step  towards  higher  sophistication  and,  therefore, 
efficiency,  of  the  WP  component.  In  essence,  a  real-time  correlator  performs  a 
transformation  of  the  detection  problem  from  the  time  domain  (with  parameter  t) 
to  a  new  one,  namely,  the  correlation  domain  (with  parameter  t),  in  the  hope 
that  the  new  "data",  i.e.,  the  correlation  function  y(x),  can  provide  enhanced 
feature  separation  between  the  signal  and  interference.  The  merit  of  such  tech¬ 
niques  will  always  be  judged  against  the  radiometer  which,  by  nature,  takes  no 
account  of  any  signaling  details. 

We  note  that  the  two  scenarios  to  be  examined  involve  the  two  most 
different  types  of  interference  from  a  spectral  viewpoint,  namely,  totally  flat 
(AWGN)  and  totally  peaked  (tones).  Furthermore,  they  are  designed  such  that 
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the  structural  differences  between  the  signal  and  interference  can  indeed  be 
exploited  to  the  fullest.  As  a  result,  significant  gains  can  be  extracted 
from  the  correlation  domain  with  respect  to  the  radiometer.  At  the  other  ex¬ 
treme,  no  intelligent  algorithm  could  distinguish,  for  instance,  a  random  hop¬ 
ping  tone  from  unrelated  extraneous  random-tone  interference,  in  the  absence 
of  additional  information.  But  even  then,  correlation  algorithms  would  per¬ 
form  no  worse  than  the  radiometer.  In  conclusion,  it  can  be  stated  that  the 
algorithms  to  be  analyzed  herein  offer  the  potential  of  performance  improve¬ 
ment  over  the  common  energy  detector,  which  can  range  from  none  to  impressive, 
depending  on  the  scenario.  The  structure  and  performance  of  decision  rules 
which  will  bring  about  those  gains  in  the  presence  of  mixed  types  of  interfer¬ 
ence,  will  be  the  topic  of  future  research. 


5.1 


Wideband  Detection  of  FH  in  AWGN 


5.1.1  System  Model  and  Proposed  Algorithm 


Henceforth,  the  term  "wideband'1  will  be  used  to  indicate  detection 
on  a  per-band  basis  under  the  cardinal  restriction  (5-4).  One  mathematical 
representation  of  a  real-time  autocorrelator  whose  actual  implementation  Is 
discussed  in  [12]  is  shown  in  Figure  5.3.  Let  us  note  that  the  offset  t, 
which  is  the  domain  parameter  for  the  output  autocorrelation  function  y(-r), 
varies  linearly  with  the  actual  time  t  of  the  incoming  waveform  as  a  result 
of  the  real-time  nature  of  the  correlating  device.  Therefore,  the  actual  in¬ 
tegration  time  is  from  t  to  TH.  Since  the  radiometer  output  is  only  y(0), 
i.e.,  just  a  data  point  in  the  new  domain,  improved  performance  due  to  the  new 
scheme  is  naturally  expected. 

The  decision  rule  adopted  here  is  as  follows:  The  output  y(r)  is 
power  sampled  at  multiples  of  B”1,  where  B  is  the  Input  BPF  bandwidth  and  the 
resulting  samples  are  linearly  combined  and  compared  to  a  threshold,  whereupon 
a  decision  is  reached  (Hj  versus  HQ  for  the  corresponding  spectral  band).  No- 
tationally. 


6-1  H1 

Y  =  y  at  |  Threshold  (5-7a) 

k=l  Hq 

where  Y  is  the  decision  statistic,  -  k*B~*,ksl,2,...  ,G-1  is  the  kth  sampling 
point  in  the  t  domain,  and 


W 


k 


y2CTk) 


|Lowpass 


(5-7b) 


G"  1  I 

is  the  power  measurement  around  the  point.  The  set  of  coefficients  {a^- 
can  be  chosen  according  to  any  particular  philosophy  and  is  subject  to  optimi-  j 

zation.  It  can  be  argued  that,  under  a  low-SNR  assumption,  setting  all  coeffl-  j 

cients  equal  to  a  constant  a^*l  is  a  near-optimal  choice.  This,  however,  is  not  J 

a  well-founded  conclusion  as  it  involves  a  number  of  approximations.  Fortunate-  j 

1 y,  simulation  has  indicated  (see  section  5.1.3)  that  performance  of  the  summa-  j 

tion  (5-7a)  is  very  insensitive  to  the  exact  value  of  a.  and  a  number  of  reason-  , 

K  I 

able  choices  would  work  as  long  as  the  upper  limit  of  the  suirenation  is  properly  • 


truncated  at  a  level  of  AG,  which  Is  less  than  G  (i.e.,  A  <  1).  This  Is 
because  the  last  samples  possess  such  an  Increased  variance  that  their  In¬ 
clusion  In  the  summation  Is  detrimental.  Analytically,  this  Is  equivalent 
to  letting  ak  =  0  for  AG+1,...,G  In  (5-7a).  Thus,  a  proper  choice  of  A  Is 
rather  crucial,  while  that  of  jakjkGj  Is  not.  For  reasons  that  will  soon 
be  clear,  it  will  be  mathematically  convenient  to  let 

a  *  I  (TH  “  Tk)  2  ;  k-1 . XG 

/  0  *,  k=AG+l,...,G 


in  the  following  calculations. 


5.1.2 


Performance  Analysis 


Here  we  shall  attempt  a  performance  evaluation  of  the  decision  rule 
(5-7)  when  the  Interference  consists  of  bandpass  AWGN  n(t).  Thus,  during  one 
hop  Interval  0  <  t  <TH,  under  Hj,  the  Input  r(t)  to  the  autocorrelator  of  Fig¬ 
ure  5.3  consists  of  a  tone  at  some  radian  frequency  uq  within  the  observed  band 
plus  noise,  1 .e. , 


r(t)  a  /Z5  cos  wgt  +  n(t) 


(5-9a) 


where  n(t)  assumes  the  typical  bandpass  representation* 


n(t) 


cos 


u0t  - 


nQ(t) 


sin 


to, 


i* 


( 5  -  9b ) 


In  (5-9b),  Nj(t)  and  Nq(t)  are  the  Inphase  and  quadrature  lowpass  noise  pro¬ 
cesses,  respectively,  which  are  approximately  independent  and  Identically 
distributed,  with  zero  mean,  flat  two-sided  PSD  of  NQ/2  W/Hz,  and  two-sided 
bandwidth  of  B  Hz  each.  Since  double  frequencies  are  filtered  out  by  the 
Integrator,  the  autocorrelator  output  y(t)  around  the  signal  frequency  Is 


y(x) 


r(t)r(t-x)dt 


ySS^  +  ySN^  +  yNS^  +  yNN^  (5-1°) 


with 

^SS  “  S(TH  “  T)  cos  “0T  (5-lla) 

ySj»|(T)  s  >^(Nj(t)  cos  uqT  +  Nq(t)  sin  u>qx)'  (5-llb) 

yNS(^)  *  cos  u>qt  -  Nq(t)  sin  u>qx)  (5-llc) 

yNN^*^  =  +  NQQ^t))  V  +  (^IQ^  ”  NQI(x))s1n  uigx  (5-lld) 

* 

Since  we  shall  be  concerned  with  noncoherent  detection,  an  arbitrary 
phase  8  can  be  Inserted  In  the  bandpass  description  of  the  noise  without  alter¬ 
ing  the  results. 


where  the  correlation-domain  noise  processes  have  been  Introduced,  l.e., 


"mfT>  • 

IPJ 

T 

npj-|(t-x)dt 

(5-12a) 

Nrn(T)  * 

Q 

/’■ 

njY)(t)dt 

(5-12b) 

Njj(t)  *> 

T 

n j ( t)  nj(t-t)dt 

(5-12c) 

Nqq(t)  * 

/'• 

nQ(t)  nq(t-T)dt 

(5-12d) 

NjqU)  = 

/■ 

T 

nj(t)  nq(t-T)dt 

(5-12e) 

Np,(t)  • 

f>q(t)  nj(t-T)dt 

(5-12f) 

T 


Much  of  the  following  theory  deals  with  the  statistical  characterization  of  the 
above  noise  processes.  For  Instance,  It  Is  clear  that  the  first  four  processes 
In  (5-12a)  and  (5-12b)  are  zero  mean  and  Gaussian,  while  the  remaining  are  not 
Gaussian.  However,  they  will  be  approximately  treated  as  such  In  the  analysis, 
particularly  for  0  <  t  «  T^,  by  virtue  of  a  central-limit  type  argument,  as 
follows:  The  bandwidth  B  of  n j (t)  and  nq(t)  is  much  larger  than  TH"1;  thus, 
each  Integral  in  (5-12c)  through  (5-12f)  can  be  approximated  by  a  large  sum  of 
noise-product  samples,  each  sample  taken  B"*  seconds  apart.  Since  we  shall  be 
Interested  in  values  of  t  that  are  multiples  of  B"1,  It  can  be  shown  that  those 
samples  are  mutually  uncorrelated.  Although  this  does  not  Imply  statistical 
Independence  (which  Is  a  prerequisite  for  a  rigorous  application  of  the  central- 
limit  theorem).  It  nonetheless  reinforces  the  validity  of  the  approximation. 
This  will  permit  us  to  calculate  higher  moments  using  Gaussian  Identities  whose 
exact  evaluation  would  otherwise  be  unwieldy.  We  furthermore  note  that  simula¬ 
tion  has  confirmed  the  accuracy  of  the  approximation  to  a  pleasantly  surprising 
degree  (see  section  5.1.3). 


The  first  point  of  interest  Is  the  mean  value  of  y(x)  In  (5-10).  We 
shall  assume  that  the  input  BP  filter  in  Figure  5.3  has  a  perfectly  rectangular 
transfer  function  ("brick  wall"),  which  implies  that  the  autocorrelation  func- 
tion  Rnn(x)  for  both  nj(t)  and  ng( t)  is  given  by 

R„„(t)  4  ^(nj(t)nj(t-x)}  »  *{nQ(t)ng(t-,)}  *  ^)SaUB'r)  (5‘13) 

where  Sa(x)  =  (sin  x)/x.  Using  the  fact  that 

'{"i}  ’  *{Nq}  -  *{«§  ■  *{NI<)}  *  «{*Ql)  -  0  (5- 14a) 

and 

*  4t{Nqq(,)}  ■  (T„  -  t)  R„n(T)  ( 5-  14b) 

we  conclude  from  (5-10)  through  (5-14)  that 

tf|y(t)j  a  [s  +  (N0B)  Sa(irBt)j(TH  -  t)  cos  u>qt  (5-15) 

Certain  interesting  observations  can  be  made  on  (5-15).  First,  although  the 
quadrature  component  (sin  oigx)  has  been  eliminated,  the  mean  *{y(t)|  still 
Includes  an  oscillation  of  unknown  frequency*.  That  is  why  noncoherent  (power) 
sampling  must  be  performed  at  the  signal -processing  unit.  Second,  we  note 
that  the  mean  value  at  zero  offset  «f{y(0)}  *  £S  +  NqB)Th  Is  just  the  average 
energy  measured  by  the  radiometer  at  the  end  of  T^  seconds,  as  It  should  be. 

Since  we  have  assumed  that  «  1,  Implying  that  S  «  NqB,  It  follows  that 
the  radiometer  output  is  dominated  by  the  noise  contribution— definitely  not 
a  helpful  situation.  Third,  and  most  important,  it  follows  from  (5-15)  that. 

If  the  output  y(r)  is  power  sampled  at  multiples  of  B”*,  l.e. ,  *  kB_1;k  an 

integer,  the  contribution  of  the  mean  noise  will  be  eliminated  due  to  the 
nulls  of  the  Sa(x)  function.  This  fact  will  be  exploited  in  the  following. 


The  only  thing  known  about  fg  *  ug/2ir  Is  that  it  belongs  in  the 
frequency  band  observed,  but  it  is  otherwise  unspecified. 
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Let  us  return  to  (5-10)  for  a  moment.  Since  we  shall  be  Interested 
in  the  low- SNR  case,  we  make  the  simplifying  assumption  that  the  signal  times 
noise  terms  y sn(t)  an<*  are  ne9^9^e  compared  to  the  noise  times  noise 

term  yNN(x).  For  Y^n  below  -15  dB  or  so,  such  a  simplification  is  quite  justi¬ 
fied,  as  also  Indicated  by  the  simulation.  When  is  higher,  the  neglected 
terms  Increase  the  variance  of  the  output  statistic  under  Hj  ( they  clearly  are 
not  present  under  Ho},  with  a  resulting  deterioration  In  overall  performance. 

At  this  stage,  however,  we  shall  derive  a  performance  estimate  for  low  SNR, 
leaving  more  refined  (and  much  more  cumbersome)  calculations  for  later.  Within 
this  formulation,  y(x)  can  be  rewritten  as 


y(x)  =  C(x)  cos  wqt  +  N^(t)  cos  oiqt  +  Nq^(t)  sin  u>qX 


(5-16) 


where 


C(t)  = 


sCTh-0 


(5-17a) 


and  the  equivalent  noises 


N®q(x)  4  Nn(x)  +  NgpU) 


Ngq(x)  4  Njq(x)  -  Nqj(x) 


(5-17b) 


(5-17c) 


As  mentioned,  y(x)  will  be  power  sampled  at  x^  =  kB  ,  at  which  points, 


r{NfOO}  *  '{Nq’M)  ■  0 


(5-ia) 


According  to  our  previous  discussion,  N®q  and  Ngq  will  be  treated  as  Gaussian 
noises.  Furthermore,  it  is  shown  in  Appendix  H  that  Nj^xiQ  and  N^(xj}  are 
uncorrelated  for  every  k,j=l,2,3. . . ,G;  similarly,  the  pairs  (Nf\Tki}t  Nl\Tk2}) 
and  (^o^ki}*  ^(^2))  are  uncorrelated  as  long  as  xiq  /  x^.  Coupling  this 
result  with  the  approximate  Gaussian  assumption,  it  follows  that  the  sequences 


of  noise  samples  and  {nq<,CtO^(cssi  contain  2G  Gaussian  rv's  which 

are  zero  mean  and  approximately  independent.  Thus,  to  complete  their  statisti¬ 
cal  description,  the  second  moments  (or  variances,  for  this  case)  are  required. 
Those  are  obtained  in  Appendix  H,  wherein  it  is  shown  that 


where 


'{("foo)2}  ■  (V)2(VTk)2  (Fi<k> +  F2<k>) 
*{("!TOk:))2j  ■  (w)2(W2  (Fi<k>  -  F2<k>) 

F1(k)  =  f1  (l-p')  Sa2[rG£l  -  Ck>]dP’ 


( 5- 1 9a ) 


(5- 19b) 


(5-20a) 


>(k)  =  (l-p')Sa^G(Cl  -  ck>*  +  Ck)]  Sa[^G((l-ck>'-ckj]dp'  (5-20b) 


and  the  normalized  parameter  $k  is  defined  as 


_  -  T|t  -  k  _  k  .  ksl  r  , 

Ck  ’  "  STjj  ‘  E  ’  k*1--M 


(5- 20c) 


so  that  0  <  ck  <  1.  Some  meditation  on  the  functional  form  of  F2(k)  reveals 
that  It  has  a  negligible  contribution  compared  to  F ^ ( k) ,  so  it  will  be  dropped 
henceforth.  From  (5-19),  we  then  have  that 


f(k>  •  '{(ntVo)2}  •  '{  kW} 


CotfUn-’k)2  Fi<k>  ■  W^O-tk^iW 


(5-21 


which  is  a  fairly  simple  expression,  but  also  fundamental  for  the  following. 


.*  >  “  >  .* 


Let  us  now  return  to  (5-16).  As  mentioned,  the  lack  of  knowledge 
regarding  wq  forces  us  to  noncoherent  sampling.  Let  Wk  -  ^(TOIlowmss  denote 
such  a  sample,  which  is  also  defined  in  (5- 7b).  We  note  that  this  kind  of  sam¬ 
pling  can  be  implemented  via  a  square-law  device,  followed  by  a  wideband  lowpass 
integrating  filter;  its  bandwidth  should  be  several  times  larger  than  the  Input 
bandwidth  B,  but  still  much  narrower  than  twice  the  carrier  frequency.  Then 
(5-16)  and  (5-17a)  imply  that,  under  H^, 


2  2 

(STH)2  (l  •  <k)  *  (NftTO)  +  (N<f  00)  i  k-1 . G-l  (5-22) 


By  virtue  of  the  approximate  statistical  Independence  of  the  second  and  third 
noise  terms  in  (5-22),  {W|()jj"j  is  a  sequence  of  approximately  Independent  noisy 
samples,  upon  which  the  decision  is  based.  As  mentioned,  this  valuable  property 


of  independence  gradually  diminishes  as  k  -*•  G-l  because  the  fundamental  assump¬ 
tion  in  assessing  it  (i.e.,  Gaussian  equivalent  noise)  weakens  toward  the  end. 
This  Is  also  why  those  latest  samples  should  be  Ignored.  Nonetheless,  such  a 
property  is  the  key  ingredient  in  the  anticipated  superiority  of  the  correlator 
versus  the  radiometer— namely,  the  fact  that  additional  information  can  be  ex¬ 
ploited  by  incorporating  all  those  new  samples  in  the  decision.  Their  approx¬ 
imate  Independence  keeps  these  samples  from  being  statistically  useless. 

In  order  to  proceed  with  the  evaluation  of  rule  (5-7),  we  need  th« 
mean  and  variance  of  the  Wk's.  From  (5-19)  and  (5-22),  we  have  that  the  dif¬ 
ference  of  the  means  of  under  the  two  hypotheses,  A*{W|<},  is  given  by 


a'W  4*{*k|Hl}  -  *|Wk|H0}  *  (STH)2  (l  -tk)! 


(5-23) 


Furthermore,  using  the  independence  between  and  NqPO'O  with  the 

Gaussian  assumption  yields  the  variance  as 


-M  -  ^{Woo)2}  *  *2{(«fW2}] 

=  2[3oN4(k)  -  aN4(k)]  =  (2aN2(k)) 


(5-24) 


with  ouc  as  per  (5-21). 
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As  seen  in  (5-23),  (5-24)  and  (5-21),  the  quantities  of  interest  for 
Wk  include  a  multiplicative  factor  (l  -  a  reflection  of  the  fact  that  the 

integration  time  reduces  proportionally  to  xk  as  xk  -*■  T^.  The  insertion  of  the 
coefficients  ak  in  (5-8)  purported  to  remove  this  factor  without  essentially 
altering  the  results.  Thus,  if  we  define 


u  ■  4 

wk 


fTH  -  \)‘ 


V(‘  ■ 


(5-25) 


it  follows  from  (5-23)  and  (5-24)  that 


'M  ■ 


(5-26a) 


Hv}  *  4(v) 


4  F,z(k) 


(5-26b) 


while  the  decision  rule  (5-7a)  reads: 


XG 

*  I 


(5-27) 


We  shall  now  define  performance  for  the  correlator  in  a  manner  that 
Is  analogous  to  the  radiometer  as  manifested  by  (5-6).  In  other  words,  within 
the  Gaussian  assumption  for  Y,  performance  is  quantified  by  (5-6a),  where  the 
output  SNR  for  the  correlator  is  now 


-|>{YlHo))2 

HvlHo} 


(5-28) 


We  wish  to  emphasize  again  that  (5-28)  Is  meaningful  only  when  Y  possesses  a 
distribution  that  Is  sufficiently  close  to  Gaussian  and  the  variance  under 
either  hypothesis  is  the  same.  The  latter  assumption  Is  justified  for  very 
low-input  SNR,  while  the  former  can  be  validated  only  though  simulation.  For 
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somewhat  higher  SNR,  the  proper  measure  of  performance  and  comparison  Is  the 
pair  (T’D'Pfa};  In  which  case,  var|Y|Hj|  Is  larger  than  var|Y|Ho|  by  a  signifi¬ 
cant  factor.  More  discussion  on  this  topic  can  be  found  In  the  following 
section. 

Based  on  (5-26),  (5-27)  and  the  uncorrelatedness  of  W^'s,  It  follows 
from  (5-28)  that 


i 


V 

I 


Fk 


SNR 


corr 

out 


CxGS2]' 


x2g2 


XG 


75" 


4CW  I  Fi  (k)  4  l  F/(k) 

k*l  1  k=l  1 

or,  using  the  fact  that  yh  =  G-Y.jn, 


»Crr  ‘  *(>.8)  S^-Th4 


where  the  coefficient  A(X,G)  Is  defined  as 

.2 


A(X,G)  = 


x‘ 

75" 


4G  I  F  ‘(k) 
k=l  1 


(5-29) 


(5-30) 
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and  will  be  called  the  "correlator  processing  gain"  or,  for  brevity,  CPG.  The 
terminology  is  easily  understood  if  we  combine  (5-6b)  and  (5-29)  into  a  ratio: 


SNR 


corr 


SNR 


out  _  2 

raff"  -  A(x ,G)  *yh 

out 


(5-31) 


The  interpretation  of  (5-31)  Is  that,  within  the  framework  of  performance  estab¬ 
lished  herein,  the  correlator  will  outperform  the  radiometer  if  the  product 

O 

A(x,G)»yh  is  above  unity  and  vice  versa.  We  note  that  this  product  signifies 
some  kind  of  utility  factor  and  Is  separated  Into  two  components:  (1)  A(x,G), 
depends  on  features  of  the  device  used  (G)  and  the  signal-processing  algorithm 
employed  (X),  and  (2)  yh  a  ST^/Nq,  depends  on  scenario  parameters  (signal  and 
noise  strength,  hopping  duration).  Thus,  the  acronym  CPG  Is  naturally  attached 


to  the  component  (gain),  depending  on  the  signal -processing  aspects  of  the 
detector.  Let  us  briefly  focus  on  it  in  the  following. 

Exact  evaluation  of  A(x,G)  unfortunately  requires  numerical  integra¬ 
tion  by  computer,  which  somewhat  destroys  any  anticipated  insight.  However,  a 
bounding  argument  can  be  developed  which,  in  essence,  warrantees  a  minimum  value 
for  A(x,G) ,  namely,  Amin(X,G).  This  minimum  value  (or  lower  bound)  can  be  used 
in  conjunction  with  (5-31)  to  warrantee  a  minimum  gain  for  the  correlator  versus 
the  radiometer  and  is,  in  fact,  very  tight  to  the  exact  result  for  low  values  of 
X.  The  argument  is  based  on  the  observation  that  the  quantity  F^(k)  in  (5-20a) 
is  a  monotonically  increasing  function  of  the  argument  ?k;  thus, 

F^k)  <Fx(xG)  =  y1  ( 1-p ' )Sa^ [jrG( l-x)p^j dp*  ;  k*l,...,XG  (5-32 

0 


since  cx6  *  x  from  (5-20c).  The  Inequality  in  (5-32)  can  be  strengthened  by  the 
following  steps: 


jfl  (l-p')Sa2[^(l-x)pJ]dp' 

<  ^  Sa2^G(l-x)p^dp'  =  G— 1  ( 1— X ) ~ 1  / 


G(l-X)  2 

Sa£[irx]dx 


<  G”1 ( 1 


/" 


Sa^[irx]dx 


=  G’1(l-x)"1(l/2) 


(5-33 


Thus,  from  (5-32)  and  (5-33), 


XG 


I  Fj  (k)  <  XG  Fj  (XG)  <  - ? 

ksl  1  1  4G(l-xr 


which,  when  combined  with  (5-30),  yields 


*  Am.  (x)  =  x(l-x)‘ 


A(X,G)  >  A„4_(X,G) 


(5-34 
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Vie  note  that  the  lower  bound  A  -n  depends  on  X  only,  while  the  exact  gain  A(x,G) 
also  depends  on  G.  This  is  just  a  fortunate  coincidence,  and  indicates  that 
Amin  can  be  used  for  every  G»  as  l0n9  as  G  »  1.  Thus,  we  arrive  at  the  follow¬ 
ing  bound: 


cNOcorr 

iNKout  ^  .x2  2 

>  X(l-X)  yH  (5-35) 

SNRout 


which,  as  mentioned,  should  be  tight  for  small  X. 

In  addition  to  this  general  bound,  we  can  create  a  good  approximate 
formula  for  the  above  ratio  when  G(l-x)  »  1,  as  follows:  since  the  steps 
used  to  derive  (5-33)  are,  in  fact,  tight  successive  approximations,  we  con¬ 
clude  that  Fj(k)  =  G“1(l-kG*1)'1(l/2),  which  leads  to 


G  . 

I  4 

m=G(l-x)  m 


The  last  summation  can  be  well  approximated  by  the  integral 


G  1 
m=G(l-x)  m 


•f 

G(l-X) 


X 


provided  that  G(l-X)»l.  Substitution  into  (5-30)  leads  to 

> 


(5-36) 


which  differs  from  the  bound  (5-35)  in  the  exponent  of  (1-x).  For  low  x, 
(5-35)  and  (5-36)  agree  closely.  In  any  case,  (5-35)  and  (5-36)  indicate 
that  the  relative  merit  of  the  correlator  increases  proportionally  to  the 
square  of  the  hop  SNR  and  could  therefore  reach  significant  levels,  depend¬ 
ing  on  the  application  (coded  or  uncoded  systems,  slow  or  fast  hopping,  etc.). 


I 
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5.1.3 


Simulation  Results  and  Discussion 


The  proposed  correlation  algorithm  plus  the  radiometer  were  simulated 
by  computer;  a  detailed  description  of  this  simulation  can  be  found  in  Appendix 
I.  Here,  we  shall  summarize  the  findings  and  discuss  them. 

First,  the  tightness  of  the  predicted  bound  (5-35)  was  examined.  The 
following  parameters  were  chosen  so  that  the  theoretical  assumptions  involved 
in  the  analysis  would  be  well  justified:  (1)  large  time-bandwidth  product  (6* 
1000  or  30  dB),  (2)  small  input  SNR  ^Y-n  <  -10  dB  j ,  (3)  small  X  (x*0.1)  and, 

(4)  large  number  of  trials  (200),  so  as  to  assume  statistical  confidence.  Since 
those  numbers  were  in  accordance  with  the  assumptions  made,  any  deviations  of 
the  simulation  from  the  theory  would  indicate  an  error  in  the  modeling  process 
(Gaussianness ,  etc.).  Fortunately,  that  was  not  the  case,  as  evidenced  by  the 


closeness  of  the  theoretical  and  simulation  results,  as  shown  in  Table  5.1, 
which  lists  the  output  SNR  difference  &SNR  (in  dB)  between  (sNRout*")^  and 
(snr5S?)„b.  both  from  theory  (5-35)  and  the  simulation.  In  other  words,  the 
theoretical  lower  bound  from  (5-35)  is 


For  X=0.1 ,  the  first  term  on  the  right-hand  side  of  (5-37)  is  equal  to  -11  dB 
while,  for  G*10^,-  (YH)dB  =  30  ^B+(Y1n)dB*  note  t*iat  simu^ation  result 
in  Table  5.1  is  slightly  below  its  theoretical  counterpart  (lower  bound)  for 
the  first  two  columns,  while  it  exceeds  the  latter  in  the  last  column.  This 
can  be  attributed  to  both  the  many  analytical  approximations  involved  and  the 
possibility  thatYfn,  in  the  first  two  columns,  is  insufficiently  small.  For 
Yin<-20  dB,  the  simulation  seems  to  agree  perfectly  with  the  theory.  In  any 
case,  the  gains  depicted  in  Table  5.1  for  the  correlator  range  from  9  dB  to 
29  dB— definitely  a  major  improvement  over  the  radiometer. 

Next,  we  look  at  the  sensitivity  of  the  decision  rule  (5-27)  with 
respect  to  the  choice  of  X.  This  is  illustrated  in  Figure  5.4,  where  the  detec¬ 
tion  probability  PD  is  plotted  versus  x(Q  <  >  <  1)  for  various  values  of  Yinand 
^FA*  while  G  = 100  is  constant.  For  each  set  of  three  curves  associated  with 
a  single  value  of  Y^n,  the  corresponding  values  of  PpA  are  10-1,  10"^  and  10-6, 
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from  top  to  bottom,  respectively.  In  all  cases  examined,  the  trend  is  clear: 
performance  is  insensitive  to  changes  of  X  over  a  very  wide  range  as  long  as  X 
is  neither  too  small  nor  too  large.  So,  any  X  in  the  range  (0.1, 0.8)  would  be 
a  suitable  choice.  Although  the  curves  shown  in  Figure  5.4  are  for  a  particu¬ 
lar  choice  of  a^'s  Cthe  ones  in  (5-8)),  an  Identical  kind  of  behavior  was  also 
found  for  other  choices.  This  is  a  rather  convenient  conclusion  since  one  does 
not  have  to  worry  about  optimizing  X  whenever  the  constants  are  changed.  Fur¬ 
thermore,  it  was  found  as  a  more  general  conclusion  that  the  overall  performance 
(Pq.Pfa)  was  insensitive  to  the  particular  choice  of  a^'s.  Thus,  one  can  set 
a^-l,  k=l,...,XG,  which  would  simplify  the  decision  rule  to  a  mere  accumulation 
(summation)  of  the  power  samples  Wk  and  a  threshold  comparison. 

A  set  of  simulated  design  curves  for  both  the  radiometer  and  the  cor¬ 
relator  are  provided  in  Figure  5.5,  parts  (a)  and  (b),  respectively.  There,  Pp 
is  plotted  as  a  function  of  Y.  (dB),  for  various  values  of  PfA  in  the  range  of 
10  to  10"  .  For  comparison  purposes,  the  Pp^-10  and  10-^  curves  for  both 
have  been  redrawn  in  Figure  5.6  so  that  the  relative  gains  of  the  correlator 
versus  the  radiometer  in  input  SNR  can  be  assessed  easier.  We  reiterate  that 
direct  invocation  of  (5-36)  or  (5-31)  would  not  be  meaningful  in  this  case, 
where  G  is  relatively  low  and  Y^  is  relatively  high,  because  the  variance  of 
the  decision  statistic  under  is  many  times  greater  than  under  HQ.  Thus,  any 
gains  in  SNR  should  be  directly  read  off  the  operating  curves.  An  alternative 
comparison  is  offered  in  Figure  5.7,  which  depicts  the  familiar  "receiver  oper¬ 
ating  characteristic"  (ROC)  curves  for  both  systems,  parameterized  by  The 

parameter  G  was  set  at  G*10^.  Clearly,  for  Y-n  >-20  dB,  the  correlator  will 
outperform  the  radiometer  by  an  amount  which  Increases  with  y^n< 

Finally,  let  us  note  that  attention  should  be  given  to  implementing 
the  almost  instantaneous  power  measurement  W^;  this  can  be  accomplished  using 
a  square-law  device,  followed  by  a  fairly  wide  lowpass  filter  whose  bandwidth 
is  many  times  that  of  the  observation  bandwidth  B  (i.e.,  its  response  time  is 
a  small  fraction  of  B"*)-  Still,  that  filter  should  be  narrow  compared  to  the 
IF  frequency  in  order  to  avoid  undesirable  power  measurement  errors. 
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Figure  5.7.  Receiver  Operating  Characteristic  Curves  (P0*Pp^)  For  Both  The 
Correlator  and  Radiometer,  Parameterized  by  y. 
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5.2  Wideband  Detection  of  FH/DS  in  Random  Tone  Interference 

In  section  5.1,  a  scheme  was  proposed  for  improving  the  detection  of 
frequency-hopping  (FH)  waveforms  In  wideband  additive  Gaussian  noise  (AWGN) 
using  samples  from  the  autocorrelation  domain.  It  was  shown  there  that,  under 
fairly  general  operational  assumptions,  an  approximate  gain  of  tt(x,G)yh  in  de¬ 
cision  SNR  can  be  achieved  over  the  energy  discriminator  (radiometer),  where 
Yh  Is  the  hop  SNR.  This  gain— being  defined  on  a  per-band  basis— has  a  direct 
positive  Impact  on  the  overall  system  detection  capability.  The  proposed  algo¬ 
rithm,  albeit  inferior  to  the  optimal  likelihood  ratio  test,  had  the  advantage 
of  greatly  reduced  (hence,  manageable)  complexity. 

The  overall  approach  has  been  motivated  by  the  recent  implementational 
feasibility  of  large  time-bandwidth-product/real -time  correlators  such  as  SAW 
devices.  The  purpose  of  this  section  is  to  show  that,  depending  on  the  specific 
scenario  at  hand,  the  use  of  such  devices  can  bring  about  impressive  gains  in 
detector  performance,  especially  when  compared  to  simpler  alternatives  such  as 
the  radiometer.  In  particular,  we  shall  consider  the  case  of  detecting  a  fre- 
quency-hopping/direct-sequence  (FH/DS)  hybrid  waveform  when  the  dominant  compo¬ 
nent  of  the  observation  "noise"  is  random  tone  interference,  arbitrarily  located 
within  the  hopping  bandwidth.  In  addition,  it  will  be  assumed  that  the  total 
power  of  this  tone  interference  is  much  greater  than  that  of  the  background  ther¬ 
mal  noise,  so  the  latter  can  be  ignored.  It  will  then  be  shown  that,  for  obser¬ 
vation  data  with  a  large  time-bandwidth-product  G  ([defined  by  (5-4)3,  a  very 
simple  algorithm  operating  on  the  output  of  a  real-time  autocorrelation  device 
can  achieve  almost  perfect  performance— in  stark  contrast  with  the  poor  perfor¬ 
mance  of  the  radiometer. 

Hybrid  schemes  have  become  increasingly  popular  spread-spectrum  com¬ 
munication  choices  due  to  the  enhanced  anti  Jam  margin  which  they  offer.  Adding 
DS  modulation  to  FH  also  improves  the  anti -intercept  capability  because  the 
"noiselike"  appearance  of  DS  makes  detection  more  difficult.  On  the  other  hand, 
multiple  tones  constitute  a  common  model  of  nonwhite  interference  and  can  emerge 
in  a  number  of  scenarios,  i.e.,  it  can  be  intentional  (jamming  of  the  band)  or 
unintentional  (multiple  users  in  a  broadcast  environment,  adjacent  rauar  sources, 
etc.).  It  is  also  conceivable  that  the  tone  interference  has  been  deliberately 
inserted  by  the  communicator  in  a  pseudorandom  manner  so  as  to  impede  the  inter¬ 
ceptor's  task,  while  it  can  be  pseudorandomly  avoided  by  the  intended  receiver. 
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In  all  of  the  cases  previously  cited,  the  tone  Interference  could  be 
filtered  out  (notch  filtering)  and  various  techniques  for  doing  so  effectively 
have  recently  been  presented  [14].  The  key  provision  of  these  techniques  is, 
however,  that  each  tone  stays  at  the  notched  frequency  long  enough  for  the  fil¬ 
ter  to  adapt;  in  the  context  of  LPI  detection,  then,  fast  FH  or  1 ow-duty-cycl e 
pulsed  interference  could  present  a  severe  challenge  for  such  systems.  It  Is 
also  clear  that,  as  the  number  of  Interfering  tones  increases,  the  required 
complexity  (number  of  filters)  wou1  !  soon  surpass  acceptable  limits. 

Here  we  address  exactly  those  severe  cases  wherein  the  number  of 
tones  per  observation  Interval  Is  both  large  and  arbitrary  (random),  with  un¬ 
known  frequencies,  durations  and  phases.  It  will  soon  be  seen  that.  In  such 
environments,  the  radiometer  is  doomed  to  fail  as  a  detector.  It  therefore 
comes  as  a  pleasant  surprise  to  conclude  that  the  autocorrelation-domain  algo¬ 
rithm  proposed  here  can  provide  excellent  performance  Independent  of  all  those 
interference  parameters.  The  section  Is  structured  as  follows:  Section  5.2.1 
presents  the  system  model  and  proposed  algorithm,  while  section  5.2.2  proceeds 
with  the  analysis  of  its  performance;  that  of  the  radiometer  is  also  outlined 
here.  Section  5.2.3  provides  the  simulation  results  and  concludes  with  a  brief 
discussion. 

5.2.1  System  Model  and  Proposed  Algorithm 

The  overall  system  structure  is  similar  to  that  described  previously. 
Therefore,  we  shall  illuminate  only  those  aspects  which  are  unique  to  the  pres¬ 
ent  scenario.  Again,  our  concern  is  the  decision  rule  performance  on  a  per- 
band  basis,  as  defined  below.  The  autocorrelation  domain  will  be  employed  to 
produce  a  decision  statistic  since  an  optimal  rule  similar  to  (5.1),  would,  in 
this  case,  be  hard  to  establish  and  too  complicated  to  implement. 

Let  the  FH/DS  hybrid  signal  to  be  detected  be  represented  by  s(t)  * 

m  c(t)  cos  oigt,  where  wg  is  an  unknown  frequency  within  the  observed  spectral 

band,  c(t)  is  the  DS  code  of  rate  R  -T  and  S  is  the  signal  power.  The  un- 

c  c 

known  interference  consists  of  M  tones  CM  is  a  random  variable  in  each  observa¬ 
tion  interval,  which  is  equal  to  the  hop  time  with  1^,  and  ^ 

denoting,  respectively,  the  power,  radian  frequency  and  phase  of  the  kth  tone. 
The  total  received  signal  in  Q3,Th3,  under  hypothesis  (signal  present),  is 
given  by 


3 


ft 


v.y 
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r(t)  55  /2!T  c(t)  cos  u)Qt  +  ^  cos^t  +  4k) 


(5-38) 


Note  the  absence  of  thermal  noise  in  (5-38),  as  previously  discussed.  The  code 
c(t)  can  be  modeled  as  either  a  random  sequence  of  independent,  identically 
distributed  ±l's  with  Pr[c(t)  *  1]  *0.5,  or  a  PN  code  with  a  full  period  equal 
to  Th.  As  per  the  ratio  N  =  T^/Tc  *  Rc/R^  Increases  (N  denotes  the  number  of 
code  chips  per  hop),  the  performance  difference  becomes  Insignificant,  a  fact 
also  verified  by  simulation. 

Let  B  denote  the  input  observation  bandwidth.  The  presence  of  the  DS 
code  implies  that  B  should  be  at  least  equal  to  Rc  or  higher,  but  definitely 
much  larger  than  R^.  Equivalently,  the  time-bandwidth  product  G  ^  BTH  »  1. 
Furthermore,  for  simplicity,  we  shall  assume  that  all  tones  have  equal  power 
Ik = I/M,k=l,. . . ,M,  where  I  is  the  total  interference  power,  and  they  are  equi- 
spaced  within  the  bandwidth  B.  In  other  words,  the  frequency  separation  |f|<  - 
f" lc+i |  between  adjacent  tones  equals  B/(M+1),  which  is  much  greater  than  TH_1, 
i.e.,  |f|c~f|(+i|  »  Th*1.  None  of  the  above  assumptions  is  critical  in  the 
forthcoming  conclusions;  they  simply  ease  the  analytical  burden. 


The  real-time  autocorrelation  operation  produces  the  output 


y(f)  * 


r(t)  r(t  -  T)dt  i  0  <  t  <  T,, 


(5-39) 


Substituting  (5-38)  into  (5-39)  and  rearranging  redefines  y(x)  as 

y(T)  *  Syc(T)cosu.0T  +  (i)CTH-^  ^cos  u^x  +  nj(x)  +  ng(x)j  (5-40) 


where  y  (t)  Is  the  code  partial -correlation  function  (a  random  variable) 

V 


t* 


ycU) 


c(t)  c(t  -  x)dt 


(5-41) 


and  nQ(T),n£;k*l,2,...,M,  are  approximately  Gaussian  (via  a  central -limit- 
theorem-type  argument),  bandpass  processes  (signal  x  interference  terms) 
defined  by 


a,  \ 

VT> 


a’(t)  cos  «ikT  -  oj  <T)  sin  »kt 


(5 -42a) 


and 


nQV)  *  |  2  ^(T))  coswqt  +  eJ(T)j  sin<V  (5-42b) 

where 


“£(■0  =  f  Hc(t)  cos[(uj0 -o>k)t  -  tjldt 

(5-43a) 

«r  A- 

T 

“k(T)  "  y  H  c(t)  sin [(]“o"wkDt ’*k]dt  * 

(5-43b) 

ek(t)  *  ^Hc(t-x)  cos[C.0-Uk3t-.k]dt  t 

(5-43c) 

ek(T)  *  f  Hc(t  -  x)  s1n[Cu0  -  uk}t  -  ♦jldt 

(5-43d) 

J  L  J 

Before  discussing  the  statistical  characterization  of  the  above  noise 
processes  (section  5.2.2),  let  us  examine  the  noiseless  (mean)  part  of  y(x)  in 
(5-4CJ,  as  shown  in  Figure  5.8.  Of  particular  interest  here  are  the  envelopes  of 
the  useful  signal  (small  shaded  triangle)  and  the  Interference  (large  triangle), 
respectively,  since  the  actual  components,  e.g.,  the  first  two  terms  in  (5-40), 
are  modulated  by  the  unknown  frequencies.  The  structural  difference  between 
the  two  correlations  is  evident:  the  OS  code  superimposed  on  each  hop  creates 


a  narrow  mean-autocorrelation  function  for  |x|  >  Tc  since  the  expected  value 
of  the  function  yc(x)  in  (5-41)  is  zero.  Contrary  to  that,  the  interfering 

,T^j .  Clearly,  then,  a  power  (nonco¬ 
herent)  sample  at  ^  *TC  would  measure  interference  only*;  this  sample  could 
be  subtracted  from  the  corresponding  one  at  Tq *  0  so  that,  under  H^,  contains 
the  full  signal  power  plus  interference.  This  subtraction  would  approximately 
cancel  the  interference  contribution  at  tq s 0  so  that,  under  Hj,  only  the  sig¬ 
nal  would  emerge  while,  under  Hq,  the  statistic  would  be  almost  zero.  Thus, 
the  adopted  decision  rule  is  (see  Figure  5.9) 


tones  correlate  for  the  whole  Interval 


y2(  0) 


LP 


*Z(Tc3 


>1 

t  A0 
LP  H0 


(5-44) 


where  aq  is  a  fixed  threshold.  In  the  absence  of  thermal  noise,  aq  can  be  set 
at  a  very  small  (positive)  level  in  order  to  maximize  the  detection  probability. 
The  next  section  shows  that  the  performance  resulting  from  (5-44)  is  excellent. 


5.2.2  Performance  Analysis 

n  ft 

Me  focus  here  on  the  statistics  of  nk(x)  and  ng(x),  defined  in  (5-42) 
and  (5-43).  Let  and  var{*}  indicate  the  mean  and  variance,  respectively. 

It  is  then  easily  shown  that 


rH(x)}  a  *{sj[(T)}  *  o  ;  j*I»Q 


(5-45) 


and  that 


var 


K<t>)  ■  varW(,)l  •  '{(“k(T)) 


(5-46) 


In  the  following,  we  shall  concentrate  on  offsets  x»  mTc;  m=0,l,2,...,N.  In 


fact,  only  x  =  0  and  are  of  interest  ([see  (5-44)3 •  Then  lengthy  manipula¬ 


tions  can  establish  the  following  facts  (for  proof,  see  Appendix  J) 


Any  x,  >  Tc  would  also  do,  especially  in  the  face  of  some  uncer¬ 
tainty  about  T  . 
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(a)  The  variance  in  (5-46)  Is  given  by 


KC"Tc3)‘ 


(N-m)T 


c  _  — 

—  Sa2  »Cfo  -  fk)Tc 


(5-47) 


where  Sa(x)  ^  sin  x/x  Is  the  sampling  function. 

(b)  Each  noise  process  has  mutually  uncorrelated  Inphase  and  quad¬ 
rature  components  (therefore,  by  the  Gaussian  assumption.  Independent). 

(c)  Noises  corresponding  to  different  frequencies  £e.g.,  nj*(t)  and 

Ct  ” 

n£(T)  with  k  f  £)  are  uncorrelated. 

(d)  Processes  n£(t)  and  ng(i),  however,  have  different  properties, 
i.e.,  n£(t)  Is  a  highly  correlated  process  (as  a  function  of  t),  while  n®(t) 
generally  Is  not;  in  fact,  samples  of  ng(t)  taken  Tc  apart  could  be  uncorrelated. 

Of  the  above  conclusions,  (d)  Is  probably  the  most  Interesting  from  a 
performance  viewpoint:  If  the  noise  n^(t)  had  been  as  highly  correlated  as 
n£(x),  the  decision  rule  performance  (5-44)  would  have  been  perfect  since  the 
same*  random-noise  sample  would  be  obtained  at  t*  0  and  t  «  T  ;  thus,  they  would 

t 

cancel  out.  This  not  being  true,  a  slight  degradation  In  performance  Is  expec¬ 
ted,  as  was  also  observed  In  the  simulation.  Nonetheless,  since  the  mean  part 
of  the  Interference  C(second  term  In  (5-40)3  does  cancel  out,  the  performance  of 
this  scheme  is  far  superior  to  that  of  the  radiometer,  which  Is  oblivious  to 
that  term.  In  fact,  the  radiometer  output  is  merely  the  value  y(0),  which  Is 
dominated  by  the  power  of  the  random  Interference  (see  Figure  5.8).  Without 
further  assistance.  It  Is  Impossible  to  determine  If  there  Is  any  signal  In  the 
total  observed  power  since  the  Interference  contribution  Is  random  and,  hence, 
unknown. 

As  mentioned  above,  the  high  correlation  of  a^(t)  Implies  that  aj[(0)*» 
“kOc^’  the  ot*'er  hand,  the  degree  of  correlation  between  0^(0)  and 

sk(Jc3  varies  with  the  frequency  difference  Af^-fg-f^.  For  the  special  case 
wherein  one  of  the  Interfering  frequencies  f^g  coincides  with  fg.  It  Is  easily 
seen  that  B^Q(0)  =  B^g(Tc^;  j»I,Q.  Those  facts  are  used  in  the  subsequent 
analysis. 


Same  within  a  totally  insignificant  change 


We  can  now  return  to  the  decision  rule  (5-44):  upon  squaring,  taking 
the  difference,  lowpass  filtering  (e.g.,  rejecting  double-frequency  terms), 
assuming  that  f^Q  =  fQ  for  one  frequency  and  using  the  above  conclusions,  it 
follows  that 


(5-48) 


(5-49) 


is  the  total  BP  equivalent  noise.  Let  us  note  that,  at  least  for  the  case  of 
a  full  PN  code  period  per  hop,  ycOc3  is  approximately  zero,  which  considerably 
simplifies  (5-48).  Furthermore,  even  for  the  random-code  model,  a  comparison 
between  analysis  and  simulation  for  the  one-tone  case  (see  next  section)  has 
indicated  that  the  impact  of  the  rv  y  fTc}  is  unnotlceable ;  hence,  setting 
ycQTc}  *  0  seems  to  be  a  reasonable  approximation  for  the  general  case.  Still, 
a  full  analysis  of  (5-48)  without  further  simplifying  assumptions  is  extremely 
complicated.  In  order  to  gain  some  Insight  here,  we  shall  focus  on  the  one- 
tone  random  interference  (M*l),  with  the  reasonable  conjecture  that  the  multi¬ 


tone  case  should  provide  analogous  conclusions. 

Since,  in  chls  case,  =  Bj(°)  a  a}(0)  =  a{CTcD  *  ai  and  n2(°) 

n2(Tc},  ^  follows  from  (5-48)  that 

*  -  ^  -H—  1  +  2Yj”^l  +  2  ^Yj-  aj.normjj 


2 


(5-50) 


c 
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where  we  have  defined  the  signal -to-lnterference  ratio  Yj 
normalized  random  variable 


S/I  and  the 


1  .norm 


\  ■  W 


cos  ♦ 


1 


N 

£ 


(5-51) 


In  (5-51) ,  cn;n«l,...,N  are  the  code  chips  and  ♦j  Is  the  random  phase  of  the 

1 


interference.  Clearly,  under  HQ,  A«0  with  probability  1.  Under  H. ,  It  Is 


shown  In  Appendix  K  that 


Mhi> 


Csth3‘ 


and 


var{4|Kj}  •  2Csth54C',iD"1 


(S-S2a) 


(5-52b) 


Based  on  (5-52)  and  the  Gaussian  assumption  about  A,  the  performance  of  this 
decision  scheme  is  predicted  by  the  detection  probability 


,corr 


-  {V?(£ +  v!  -  vj 


(5-53) 


where  Ag*  »  Ao/(STh)2  Is  a  normalized  threshold  and  Q(x)  is  the  Gaussian  inte¬ 
gral  function.  In  the  absence  of  thermal  noise,  Aq*  can  be  set  arbitrarily 
close  to  zero;  thus.  It  always  yields  the  zero  false-alarm  rate  Pp^«0.  In 
practice,  Aq*  would  be  set  according  to  the  thermal-noise  level  and  the  degree 
of  uncertainty  about  the  power  S  of  the  detected  signal^. 

The  performance  predicted  by  (5-53)  Is  Indeed  excellent.  We  note 
that  Increasing  the  Interference  power  (hence.  Increasing  Tj-1)  actually  he! ps 
detection  Instead  of  deterring  It,  while  It  has  no  effect  on  false  alarm.  For 
Instance,  It  can  be  shown  that.  If  N  >  60,  then  PQ0rr  >  99%,  independently  of 


^Alternatively,  this  can  be  expressed  In  terms  of  the  uncertainty 
regarding  the  transmitter/intercept  receiver's  true  distance. 


$ 


Yj,  as  long  as  Ao*  <  0.4,  1  .e.,  one  can  tolerate  40%  uncertainty  about  the 
signal  power  and  still  expect  excellent  detection  capabilities,  regardless  of 
the  Interference  power.  The  above  conclusion  Is  rather  Insensitive  to  N  In 
that,  for  N*10,  the  corresponding  minimum  pjjorr  Is  90%. 

The  radiometer  performance  (shown  In  Figure  5.2)  Is  easier  to  analyze 
for  an  arbitrary  M  and  Is  based  on  the  fact  that  Its  output  Yra<^  «y(0)  can  be 
written  as 


V.d  ■  j,  '5-M> 

Let  M  be  the  maximum  number  of  tones  which  can  be  expected  In  any  hop.  Here, 

RlpX 

for  simplicity,  we  consider  only  the  case  M  »1,  where  Prob[one  Interfering 

max 

tone]  3  Pr[no  interfering  tone]  ■  1/2  for  each  hop.  Assuming  that  the  thresh¬ 
old  Yq  Is  set  at  YqMjT^,  so  that^  PFA*®*  **  can  then  1)6  shown  that  (s*e 
Appendix  L) 


Thus,  the  obtainable  performance  decreases  with  decreasing  signal -to-lnterference 
ratio  Yj  and  In  the  limit 

"•  pDad  ■  i  CUx  ■  \  (5-56) 

Yj  -*-0 

rorr 

which  Is  certainly  poor  compared  to  Pjj  .  It  has  also  been  shown  that  the  ra¬ 
diometer  performance  is  a  decreasing  function  of  M  so  that  the  above  results 

max 

constitute  an  upper  bound  for  the  general  case.  Finally,  we  note  that  the 

^Thls  Is  selected  so  as  to  match  the  zero  Pp^  of  the  correlator; 
another  choice  of  Yrt  would  lead  to  P,-.  3  0.5. 
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question  of  threshold  setting  Is  much  more  crucial  for  the  radiometer  than  the 
correlator  since  the  performance  of  the  latter  is  effectively  independent  of 
the  jamming  power. 

5.2.3  Simulation  Results  and  Discussion 

The  correlator  with  one  tone  interference  was  simulated  by  computer; 
the  results  shown  In  Figure  5.10  are  based  on  10,000  Independent  trials.  (For 
simulation  details,  consult  Appendix  M.)  Also  shown  In  this  figure  are  the  an¬ 
alytical  predictions  (dotted  lines)  whose  agreement  with  the  simulation  Is 
quite  striking,  even  for  such  low  values  of  N  as  N«3  and  N*10.  Since  detec¬ 
tion  probability  Is  monotonically  Increasing  with  N  (see  (5-53)),  those  values 
represent  worst-case  designs  which  nonetheless  yield  excellent  performance  re¬ 
sults.  It  was  somewhat  surprising  to  find  that  the  Gaussian  model  provides 
such  an  accurate  analytical  prediction,  even  for  N*3;  furthermore,  setting 
yc(Tc)  equal  to  zero  proved  to  be  a  well -justified  simplification.  Note  that 
no  false  alarm  was  observed  (Pp^  *  0)  and  that  performance  Is  practically  Insen 
sltive  to  the  amount  of  interference  inserted.  In  contrast,  the  radiometer 
performance  (as  evidence  by  (5-55 J) deteriorates  rapidly  with  decreasing  S/I,  as 
expected.  Finally,  let  us  mention  that  the  performance  shown  In  Figure  5.10  Is 
for  a  nonoptimlzed  (arbitrarily  chosen)  threshold  Aq*»S  TH  /N.  In  the  absence 
of  thermal  noise,  further  improvement  can  be  attained  for  the  correlator  by  de¬ 
creasing  Aq*  to  a  very  small  (but  positive)  value. 

Although  analysis  and  simulation  are  not  yet  available  for  the  multi- 
tone  case.  It  is  anticipated  that  the  gap  between  the  radiometer  and  the  corre¬ 
lator  performances  will  Increase  as  the  Interference-to-nolse  ratio  Increases, 
Independently  of  That  again  Is  due  to  the  relative  Insensitivity  of  the 

correlator  to  the  Interference  nuisance  parameters. 

Comments  pertaining  to  the  power  measurement  of  (5-44)  are  the  same 
as  those  in  section  5.1.3.  Furthermore,  we  note  that.  If  the  carrier  frequency 
fg  of  the  signal  sought  is  known,  narrow  bandpass  filtering  of  the  output  y (t ) 
In  (5-40)  prior  to  the  power  measurement  will  further  enhance  performance. 

As  mentioned,  the  theory  of  this  section  does  not  include  thermal 
noise.  The  determination,  analytical  evaluation  and  simulation  of  algorithms 
which  operate  satisfactorily  In  a  mixed  environment  (i.e.,  both  thermal  noise 
and  random  interference)  will  be  the  topic  of  future  research. 


Truncation  at  0.9999 


Real-Time  Autocorrelator  (Nonoptimized  Threshold) 
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APPENDIX  A 


PROOF  OF  EQUATION  (3-8) 


Let  us  start  from  (3-5),  i.e.. 


A  (r(t))  - 


exp 


-NY. 


2N 

S, 


exp 


H 


1 


r.  c. ft 
j»l  J  H0 


(A-l) 


with  yc  as  per  (3-6),  c..  denoting  the  jth  chip  of  the  1th  code  pattern  and 

J 

r.  as  per  (3-7).  Consider  the  N-dimensional  linear  space  consisting  of  all  bi- 
J  N 

nary  patterns  c- ;i=l,. . . ,2  of  length  N.  Now  observe  that  this  space  can  be 
divided  into  two  disjoint  complementary  subsets,  each  of  cardinality  2  ~  ,  by 
the  rule  that,  for  each  possible  vector  c^  belonging  to  one  of  the  subsets, 

-c..  belongs  to  the  other  subset.  Exactly  which  vectors  c^  are  included  in  which 
subset  is  immaterial  as  long  as  the  above  rule  holds.  We  can  now  perform  the 
summation  in  (A-l)  ove;  one  of  the  two  complementary  subsets  instead  of  the 
whole  linear  space,  with  the  equivalent  rule. 


i3l 

ie» 


H1 

r  Ao 

Hn 


(A-2) 


where  r  =  , . . . ,rNj,»  is  one  of  the  aforementioned  disjoint  subsets  and  any 
resulting  scale  factors  are  absorbed  into  the  threshold  A'^.  in  (A-2),  cosh(x) 
is  the  hyperbolic  cosine  function 


x  -x 
cosh(x)  =  - - 5-^ — 


(A-3) 


Since  cosh(x)  is  an  even  function  of  x,  the  summation  in  (A-2)  can  be  extended 
over  the  whole  linear  space 


,-svv\ 


a 


^  A<> 


(A-4, 


In  order  to  proceed,  we  first  need  to  introduce  the  superscripts  r^  and  c|N^, 
which  will  indicate  the  length  of  the  corresponding  vectors,  and  then  prove 
the  following: 

Lemma  A.  It  is  true  that: 


I  cosh  (r<">.c(N>)  =  2"  TT  cosh/r.) 

i=l  jsj  v 


(A-5) 


Proof:  The  proof  is  by  induction.  First,  let  N»l.  Then, 


2 

]T  cosh  »  cosh|r^)  +  cosh  |-r,  j  =  2cosh|r*j)  (A-6) 


since  cosh(x)  is  even  function.  Thus,  (A-5)  is  satisfied  for  N=l.  Next,  assume 
that  it  is  true  for  N=k,  i.e.. 


^  cosh  (rfk)-c(fk))  =  2k  JT  cosh  (r^) 


(A— 7) 


and  prove  it  for  N=k+1.  Indeed, 


2j  cosh  (r<k+1>.c<k+1>)  .  I  gosh(rfk)-c«k>  *  rktl)  *  cosh(r(k>.c‘ 


(A-8) 


cosh(a+b)  + cosh(a-b)  = 


e  +  e 


+  e 
T“ 


+  e 


*  e'8  -  (<•*  ♦«•*) 


=  2  cosh(a)cosh(b) 


which,  upon  substitution  in  (A-8),  yields 


,k+l 


,2  cosh(r.lk+1).c<k+1))  =  ^  2  cosh(r(k).  c|k))  cosh(rk+1) 

2k 

=  2  cosh(rk+1)^  cosh(r(k).c1(k)) 


k  k  lf+1 

2  cosh(rk+1)2  JT  cosh(rj)  a  2  IT  cosh(rj+1)(A-10) 
j  1  j  1 


where  the  truth  of  (A-8)  has  been  used.  Equation  (A-10)  completes  the  proof  of 
Lemma  A. 

If  Lemma  A  is  applied  to  (A-4) ,  it  renders  it  equivalent  to 


:? 

S 


(A-ll) 


which  is  (3-8). 


v'  l’ 


APPENDIX  B 

PROOF  OF  EQUATION  (3-10) 


Let  us  rewrite  the  decision  rule  (3-9)  as 


n  i 

'  yJ  $  A° 

j*l  J  H0 


(B-l) 


where  the  independent,  identically  distributed  random  variables  (iid  rv's)  y.  4 

2  J 

r.  assume  the  form 
J 


J  FTcV"/ 


under  H, 


under  H. 


(B-2) 


In  (3-14),  n.  represents  the  Gaussian  rv 
3 


"j  '  J  1 


(t)  dt 


(B-3) 


(j-l)T 


whose  statistics  are  independent  of  ,j. 

A  precise  analysis  should  account  for  the  exact  distribution  of  y. 

J 

under  either  hypothesis.  So,  for  instance,  under  Hq,  y^  is  chi -squared  rv 
with  N  degrees  of  freedom.  For  large  N  (of  the  order  of  hundreds  or  more,  that 
is  the  practical  case),  very  satisfactory  approximate  results  can  be  obtained 
which  circumvent  the  difficulties  of  an  exact,  but  enormously  complicated,  anal¬ 
ysis.  This  is  done  by  invoking  a  central -1 imit-type  argument:  The  decision 
statistic  A',  being  the  sum  of  a  large  number  of  lid  rv's  y.,  is  approximately 

J 

Gaussian  distributed,  with  mean 


t'Kf  a  K<rtalM  ;  k’0-1 


(B-4a) 


and  variance 


a&oisfi 


under  either  hypothesis.  Thus,  having  the  second-order  statistics  of  the  in¬ 
dividual  rv's  y.  and,  subsequently,  of  A'  from  (B-l)  enables  us  to  determine 
J 

performance  immediately. 

Let 


m  =  ^c  cj  6kl 


(B-5) 


so  that 


y.  =  (m  +  n.j  *  m1-  +  n.1-  +  2mn . 
'  J/  J  J 


(P-6) 


Furthermore,  the  Gaussian  rv  n.  has  zero  ,e?  and  variance 

3 


4  ■  4/1 


N0Tc 


(B-7) 


Thus,  from  (B-6) , 


•hi 


”2  *  4  ■  STc2 


Finally, 


from  (B-7)  and  (B-9), 


)  (? *  4 S 

1  =  (N0  Tc) 

(|  +  Vu) 

vHn/  +  z”nji 

=  vrjnj2! 

+  4m2  varj 

and  n^  are  uncorrelated  |<f| 

"/!  s  °) 

=  2on2[°n2+  2m[ 

■  (N0Tcf 

\j  *  Vklj 

(B-8) 


Equation  (3-10)  results  from  combining  (B-8)  and  (B-10)  with  (B-4). 


APPENDIX  C 


PROOF  OF  EQUATION  (3-20) 


Let  the  received  waveform  be 


r(t) 


I  c.  pCt  -  jTcO  cosCuQt  +  ♦J  +  n(t) 


n(t) 


ChxD 

CHo) 


0  <  t  <  Th 


(C-l) 


where  is  a  sequence  of  independent  phases  and  the  synchronous  model  (  e  »  0) 
is  used.  Let  c  3  ,c2, . . . .c^J  and  4  3  p^x . *nJ  denote  the  random  vec¬ 

tors  whose  components  are  the  successive  random  chip  and  phase  values,  respec¬ 
tively.  Then,  by  virtue  of  the  assumed  independence,  the  generalized  likeli¬ 
hood  ratio  reduces  to 


(f(r(t)|Hlfc,*  ;  0  <  t  <  T  )} 
Ar(t)  3  *  - L-1-,,--- - - - 


TT  .  -  j*(r(t)|H1,cJ,*;  (j-l)Tc  <t<jTc); 


JT 

J=1  J  J 


(f  (r(t)|H1,cj;<(>j;  (j-l)Tc  <t<JTg^ 
f(r(tJfH0) 


J  rj)|  «i>K! 


Hl 

*  A0 
Ho 


(C-2) 


where  the  jth  chip-interval  envelope  r.  is  given  by 

J 


[I] 


t 

ij 


and  rjj ,  rq^  have  been  defined  i*»  (3-18).  Furthermore,  we  note  that,  although 

rT .  and  rn.  depend  on  c.  =  ±1  through  the  combination  of  (3-18)  and  (C-l) ,  their 

/  J  2  \ 

squares  do  not  (simply  because  Cj  3  lj  .  Thus,  rj  of  (C-3)  is  Independent  of 


c.,  which  means  that 
J 


■  '°(vrj) 


substituting  (C-4)  into  (C-2)  yields 


*(r(t))  »  exp{-NYcJ  TT  l„  i1  A0 


whose  equivalent  logarithmic  version  is  (3-20). 


(C-4) 


(C— 5) 


APPENDIX  D 


PROOF  OF  EQUATION  (3-23) 


A  Gaussian-approximation-based  analysis  of  the  low-SNR,  chip-noncoherent 
detector.  Illustrated  In  Figure  3.3,  can  be  performed  following  guidelines  similar 
to  those  in  Appendix  B,  as  follows:  Let  the  bandpass  AWGN  n(t)  be  represented  by 

n(t)  =  i/Z’-jnjU)  cos  wgt  -  ng(t)  sin  iDgt|  (D-l) 

where  n^(t)  and  ng(t)  are  baseband,  independent,  Gaussian  processes  with  a  flat 

one-sided  PSD  of  NQ  W/Hz  each.  Conditioned  on  the  unknown  phase  and  the 

±l-chip  c.  of  the  interval  (j-l)T  >  t  >  jT  ,  the  Inphase  and  quadrature  variables 
J  c  c 

rj .  and  rq.  of  (3-18)  are  Gaussian  with  means 

J  J 

2 

and  common  variance  o 


2  a  ( 

*  var| 

where  was  defined  In  (3-10c).  Thus,  averaging  the  phase*  ^  ,  it  follows 
that  each  of  the  Independent  envelope  rv's  rj  of  (3-21)  has,  under  ,  a  Rician 
distribution 


Again  it  is  easily  seen  that  the  value  of  c.  becomes  irrelevant  due 
to  the  squaring.  J 


^  Tc  cj  cos 


*  Tc  CJ  \l  s,n 


.  ,N  /n 

k*Ofl  C) 


■j}  ■  ,arK( 


N0Tc 


(D-3) 


c 


Assuming  that  an  instantaneous  power  measurement  can  be  formed  at 


time  t  =  0  and  i  =  T  ,  we  have  the  decision  rule: 

c 


/  N  \2  Hj 

'  (w=t)  PCTc3  Z 


(M— 11) 


with  P(t)  the  power  at  time  t  and  with 


S^T  ^ 

i  iH 


(M-12) 


The  factor 


(A): 


is  used  to  make  a  *  0  under  Hg.  Each  chip  of  c(t)  is  modeled  as  an  independent 
random  variable  with  values  ±1,  each  with  probability  0.5. 

Figure  5.10  from  the  text  presents  the  probability  of  detection  (Pq) 
for  various  values  of  N  at  a  false-alarm  probability  of  zero.  Each  point  used 
to  determine  the  plot  of  PQ  utilized  10,000  trial  runs.  Besides  the  code  se¬ 
quence,  the  phase  <|>j  was  randomized  over  0  to  2ir  radians. 


fD 


Ml 


,N  (D-4) 


where  yc  is  the  predetection  SNR  while,  under  Hq  (l.e.,  signal  absent,  Yc  »  0) , 
(D-4)  reduces  to  a  Rayleigh  distribution 


f 


(D-5) 


Therefore,  the  statistics  of  the  decision  rv  /K'  in  (3-22),  being  the  square 
root  of  the  sum  of  N  iid  rv's,  can  be  precisely  found.  For  example,  under  Hq, 
it  is  a  chi -squared  rv  with  2N  degrees  of  freedom,  resulting  in  a  false-alarm 
probabil ity 


PFA  "  TiTTyr 

while  the  corresponding  detection  probability  is  given  in  terms  of  the  general¬ 
ized  Marcum  Q  function 


jf 

x‘/2 


xN_1  exp[  -x]  dx  =  1 


-  (incomplete  Gamma  function)  (D-6) 


PD  =  qn^’x0^ 


/”x(?)  exp(-  IN-i^x)dx 

x0 


(D-7) 


In  the  above  equations,  xn  is  a  normalized  threshold 


Although  the  above  results  are  exact,  they  are  cumbersome  to  use  from 
a  computational  viewpoint,  expecially  for  large  N  (N>100).  Besides,  the  Gaus¬ 
sian  approximation  about  A'  becomes  sufficiently  tight  in  that  region  and  thus 
emerges  as  an  attractive  and  simple  tool. 

Let  us  start  from  the  fact  that 


rij  ■  Vkic°s*j  +"ij 


(D-9a) 


and 


rQ.  *  *Wkl 

V 


(D-9b) 


In  (D-9),  the  Gaussian  rv's  n^  and  nq^  are  defined  by 

,'3T' 

n. 


r 


f  c 

J  nr(t)  dt 
( j”1 )TC 


(D-lOa) 


and 


/  Mt)  dt 
(j-l)Tc 


(D-lOb) 


and  are,  therefore,  mutually  independent,  each  with  zero  mean  and  variance  o  . 

2  n 
Then,  the  squared  envelope  r^  is 


rj  =  rIj  +  rQj  =  STc  6kl  +  nl.  +  nQ.Z  +  2/5>  Tc  cj  6kl(nIj  ♦j  +  nQj  sin  +j)  (D-H) 


It  immediately  follows  that 


=  ST"2  6l'1  +  2°“2  =  ST‘2  +  NnT' 


‘c  kl  n 


0,  ,  ^  I 

c  kl  0  c 


("oTc>  t1  *  \  ski)  (D-12) 


which  is  ( 3- 23a ) . 


In  order  to  calculate  the  variance  of  r^  in  (D-ll),  we  note  that  the 


first  term  is  a  constant  (thus,  it  can  be  neglected),  while  the  second,  third 

2  2 

and  fourth  terms  are  pairwise  uncorrelated.  This  is  because:  (a)  nT .  and  oq. 


•Il 


are  functions  of  Independent  rv's  and,  (b)  the  fourth  term  Includes  both  c^  and 
which  are  Independent  of  the  noise  and  zero  mean.  As  a  conclusion. 


Hr/}  s  *{(r/  -4/})2}  ■  »rKVvirM 


♦  4ST2  su  cos4j  +  "Qj  stnaj 


V  *  4STc2  Skl  ■  ("oToli*2STc,"o'kl 

(N0Tc ) 2  f1  *  2Vkl) 


(D-13) 


since 


•  jcos  a  var  jsin  3  j 


Equation  (D-13)  is  identical  to  equation  (3-23b). 


APPENDIX  E 


PROOF  OF  EQUATIONS  (3-44)  and  (3-45) 


First,  let  us  note  that,  once  Tc  is  substituted  by  T./2,  (3-43)  is 
derivable  from  equation  (3-10),  i.e.. 


Y  6. . 
c  kl) 


and 


i* 


V 


J 


as  per  (3-43).  We  shall  now  show  that 

‘,{rJ2  '  0 

and  2 

*{rjl  rj2|Hk;es°}  *  5kl 

Consider  the  situation  in  Figure  E.l.  From  (3-41), it  follows  that 

rj2  =  ^  (t)  cj  6kl  +  nj2 

rj+l,l  =  /5('f)  cj+l  6 kl  +  nj+l,l 

where  noises  n^  and  n^  ^  are  mutually  independent,  as  are  and  c 
Thus, 


(E-l) 

(E-2) 


(E-3) 

(E-4) 

j+1* 


#{rj2  rj+l,4  "  S("r)  6kl  *{cjcj+l}  +  <f{nj2  nj+l,l}  =  0 


In  an  analogous  way. 


r 


jl 


(E-5) 


1 


y 


which,  combined  with  (E-3)  yields 


(E-6) 


by  virtue  of  independence  between  and  n^.  Finally,  substituting  (3-43) 
and  (3-44)  into  (3-42b)  yields 


*{A’} 


which  is  (3-45). 
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APPENDIX  F 


PROOF  OF  EQUATION  (3-46) 


We  follow  a  procedure  similar  to  Appendix  E;  however,  the  picture  is 
changed  to  Figure  F.l,  under  the  assumption  of  a  worst  epoch  e  =  1/4.  Then,  we 
can  write  that 


rj2  ■  /T 


5.  •.  + 
kl  j2 


(F-l) 


and 


rji  =  [c j-i +  cj] +  nj 


(F-2) 


Note  that  r^2  is  identical  to  (E-3);  thus,  <£<jrj22  ;  e*l/4 j  is  as  per  (3-43) 

Furthermore, 


4jl2|Hk;e5l/4}  =  (?)  ^{rjl2|Hk;e=1/4;Asame}  +  #{rjl2|Hk;e’1/4;Adiffjf 


(F-3 


where  the  conditioning  events  and  A^.^  refer  to  c^  and  c^  having  the 

same  or  a  different  sign.  Clearly, 


Pr  {Asame}  -  Pr  {Adiff}  ~  7 


(F-4) 


a  fact  used  in  (F-3).  Let  us  now  note  that,  under  Ad^ff,  the  mean  of  r^  is 
zero  while,  under  A  »it  is  equal  to  T _/2.  Thus, 

Sadie  C 


fjl2  lHk-e'l/4}  *  {  kl)  *  (^j  ■  (^)  [  +  7  YC  \l] 


(F- 


r. 


jl  rj2 


Hk;e=l/4 


In  exactly  the  same  way, 

I  ■  '  'lil’H  "'St  "i; 

■ 

which  completes  the  proof  of  (3-46). 


I 


kl 


Definition  of  r^'s,  Assuming 


APPENDIX  G 


PROOF  OF  EQUATION  (3-63) 


Let 


N 

A'  =  I  Vj 
j=l  ] 


H1 

*  AC 


be  the  decision  rule  (3-62a),  where 


y,  -  r.(‘>  r.(2> 

J  J 


and  r.^;i=l,2  is  given  by  (3-62b).  Then, 
J 


(0-1) 


(0-2) 


rJ(i>  *  *'cc3  Skl  ‘"j'0  (M) 

where  n^1)  aid  nj^  are  zero  mean»  mutually  independent  Gaussian  rv's  with 
common  variance  a  2  =  nqTc/2.  From  (G-2)  and  (G-3),  it  follows  that 


■  S  Tc2  6kl  *  ^  Tc  cj  5kl  (nj(1>  *  "j<2))  +  "j(1)  nj(2)  (S-4) 


Thus, 


<? 


{yj}  -  5Tc2ski  ■  ("oTcK*ki 


which  justifies  (3-63a).  Furthermore, 


var 


W  • 


STc  sklvar 


l"j(l)  *  "j(2,i *  var["j(1)  nj<2). 


S  Tc2  ‘kl  2°n  *  %  *  4cn4 


(lY 


ST 


2a_ 


kl 


"  (N0Tc)2|j  +  Yc  6kl] 
as  per  (3-63b). 


(0-5) 


APPENDIX  H 


STATISTICAL  CHARACTERIZATION  OF  NOISE  PROCESSES 
N®q(t)  AND  N®q(x) 

We  define  N®q(x)  and  N®q(x)  as  per  (5-17)  and  (5-12).  It  then  follows 
from  (5-13),  (5-14)  and  (5-15)  that  *Tk  =  *cB"^  are  zer°-mean 

random  variables,  i.e.,  (5-18).  Furthermore,  by  virtue  of  a  central -1  imit- 
type  argument  (summations  of  a  large  number  of  rv's),  they  are  assumed  to  be 
approximately  Gaussian. 

We  shall  first  show  that  and  Ngq(jj}  are  uncorrelated  for 

every  k  and  j;  in  which  case,  they  are  also  approximately  independent  (due  to 
the  Gaussian  assumption).  Indeed,  from  (5-17), 


*{NiqCTkD  nq\t  j)}  *  { N i r CT niqCt j  3}  +  !NQgCTk)  NiqCTp} 

-  *{"nCO  "qtM  -  «{mqqOiO  "giM  t"-D 

Let  us  consider  the  first  term  on  the  right-hand  side  of  (H-l);  a  similar  line 
of  argument  applies  to  the  next  three  terms.  By  definition. 


^niiCtO  niqM  -  J  f H  <?{niCti)niCtrTk3niCt23nQCt2‘Tj3}dtidt2  (H_2) 


Tk  Tk 


However,  ng(t)  is  independent  of  nT(t),  which  implies  that  the  integrand  of 
(H-2)  contains  the  factor  ^{ngCt2"Tp}*  1S  zero-  Thus, 


'{"iiWiqM  *  0 


(H-3) 


Along  with  (H-l)  and  the  above  remarks,  (H-3)  establishes  the  uncorrelatedness 

nf  Ne<V_  "\  anH  NeCV_ 


of  N^Ctk3  and  Nq  (jj}  • 


i  !*,«%• 1 


Let  us  now  consider  NjH(rk£)  and  Ni  (rk2)  for  rkl  j*  t|c2 .  We  have  that 


■  'KiC^iC^dMv^iAqC^:)} 

*  <f{NIlCckiDN()QC'lt2D}+  <f{NQQC,klDNuCTk2)i  (H‘4) 

But, 

<y{NnC'rkiDNn  CTk2D}  ■  /■/■  <f{l'lCtl>lCtl  •  ,ki>lCt2>lCt2  -  Tk2)}dtldt2 

'kj^  ’k!  \2  (H-5) 

We  shall  approximate  the  integral  (H-5)  as  the  double  sum 


where  At  =  ti=B~*.  Approximating  integrals  as  per  (H-5)  by  summation  is  a  com¬ 
mon  and  well-founded  practice.  We  note  that  n j ( k )  and  n^k-k]^)  are  independent 
rv's  since  kj^O,  as  are  nj(m)  and  ni(m-k2).  Furthermore,  due  to  the  indepen¬ 
dence  of  the  rv's  involved,  the  expectation  in  (H-6)  will  be  zero  whenever  all 
four  sampling  times  k,  k-kj,  m,  m-k2  are  pairwise  different,  as  shown  in  Figure 
H.l(a).  If  we  assume  that  k^  <  k2  ,  we  then  observe  that  the  three  remaining 
cases  depicted  in  Figure  H.l  are  (b)  k-k^=m-k2;  in  which  case,  k/m,  or  (c)  k=m; 
in  which  case,  k-k^m-kg,  or  (d)  k=m-k2;  in  which  case,  k-k^.  In  all  those 
cases,  the  common  conclusion  is  that  at  least  two  rv's  exist  which  do  not  coin¬ 
cide  with  the  others.  The  net  result  is  that  the  expectation  in  (H-6)  will 


a. r.  ys  be  zero,  establishing  the  approximate  uncorrelatedness  of  NliCTkO  and 
N I T CTk2^ *  T*ie  same  exact  argument  establishes  that  <?|NQq(jk£)  NQQ(rk2}j  s  °* 
while  it  is  even  easier  to  show  that  *jNuCTkO  NQQCTk2}|  s  NQQCTklDNlICTk2D| 

=  0.  Thus  all  of  the  above  can  be  combined  in  (H-4)  to  illustrate  that  differ¬ 
ent  samples  (tk^  ^  Tk 2]  the  same  process  N®^(t)  are  approximately  uncorrelated. 
It  can  similarly  be  demonstrated  that  the  same  is  true  for  Nn^{T)»  defined  in 


Finally,  we  are  concerned  with  the  variance  of  and  N^tjQ  , 
i.e.,  (5-19).  Clearly,  since  the  respective  means  are  zero,  those  variances 
coincide  with  the  second  moments.  Furthermore, 


*j(NiqCTK))[  3  <f{NnCTk)l  +  ^{nqqCt O}  =  2<?{NnCxk)i 


(H-7) 


-  r  p  ‘?jnlCtDnlCtl-xlOnlCt23niCt2-xk)[dt1dt2 


Tk  Tk 


l+Tk)nl(t‘l)nI (t2+'rk)nI  (^2)}  dtl  dt2 


(H-8) 


where  a  simple  change  of  variables  has  been  performed.  We  shall  now  use  a  fa¬ 
miliar  property  of  four  jointly  Gaussian  rv's  (xf}|=1,  namely,  that  [15] 

<?jlTXij  =  dr{Xlx2}  <*r{x3x4}  ^{x1x3}42x4}  (H-9) 

Applying  (H-9)  into  (H-8)  results  in  (for  any  x): 

Th-t 

dt'ldt2 

Th-t 

+  nj(t^+T)|  njCt23}  dtldt2 

Th-t 

+  yy*^{ni  Cti+‘0niCt2)}<y{niCt,D  niCt2+‘0 }  dtidt2 


tH-x 

Rnn<l>  cW  +  /  f  ^‘1-9  dtidt2 

T  0 
Th-t 
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+  /  /  RnnCtrt2+‘0  "nr^i^'D  dtidV 


(H-10) 


where  Rnn(x)  s  <^nj(t)nj(t+x)£  is  the  correlation  function  of  n j ( t ) »  as  per 
(5-13). 

Equation  (H-10)  can  be  further  simplified.  First,  we  note  that,  for 
t  =  xk  =  kB-1,  Rnn(x|^  =  0,  so  the  first  term  drops  out.  Furthermore,  one  can 
employ  the  even  symmetry  of  Rnn(x),  i.e. ,  the  fact  that  Rnn(x)  *  Rnn(-t)t  in  or¬ 
der  to  reduce  the  two-dimensional  integrals  to  one-dimensional  ones  ([see  [15], 
page  325  for  details).  The  result  is  then 


4i2ic4’  "Wv^y)  Sa2(’B‘,)dp 


0 

+  2 


if1*  T,CCTh  -  xk  -  Sa^p+xo]  Sa^BCp-T^jdp  (H-ll) 


With  the  change  of  variables 


and  defining  the  normalized  parameter  ?k  as 

r  A  Tk  _  k  _  k 
k  BTh  tr 


(H— 12) 


(H-13) 


we  can  rewrite  (H-ll)  as 


Kioo}  ■  ^-cm2  f1  d-p,)sa2[.Bp'CTH-Tk)]dp' 


f1  (1-P')Sa[’B(p'CTH-’k)+Pk)]Sa[*B(p'CTH-|c)-'k)l 


(H-14) 


where 


(k)  =  P  (l-p')Sa2^BTH^l  -  ^p'  dp'  =  (l-p'jSa^GCl-cOp^p'  (H-15) 


F2(k)  =  f1  (IV)  Sa[4'Cl-Ck)  ♦  gJsa^p-Cl-cO  -  ck||*'  (H-16) 

exactly  as  In  (5-20).  Combining  (H-14)  with  (H-7)  yields  (5-19a). 

Quite  similarly,  it  follows  from  (5-17c)  that 


|  =  2#{NIQ(t)}  "  2^{NIQ^t^NQI^t^ 


(H-17) 


Njq(t) 


•n 

ff  ■»{"iCtO"qCti-')"iCt2)nQCt2-T)f<itidt2 

tH  o  Th_t 

JjK nCtr‘£)dtldt2  -  ff  "nnCM-ta  d4d4 


(H-18) 


since  nj(t)  and  ng(t)  are  Independent,  with  the  same  correlation  function.  The 
integral  in  (H-18)  is  Identical  to  the  second  term  in  (H-10),  which  corresponds 
to  the  F1(k)  term  in  (H-14).  Furthermore, 

tH 

*{N,q(,)  Ng,(t)}  ■yy4IC‘i)"qCti-'DnQCt2D  "qC4-*i  dtidt2 


’  ff  Rnn(ti  *  4  +  T)  Rnn(ti  ‘  4  '  *)d4‘ 


(H-19) 


^s-7vV> 


which  is  identical  to  the  third  term  in  (H-10)  corresponding  to  the  FgU) 
term  in  (H-14).  Combining  the  above  with  (H-17)  (note  the  minus  sign  of  the 
second  term)  results  in 


*j^(T)2j  -  Cnob)2(th  '  Tk)2  (Fl(k)  '  F2tk)) 


as  per  (5- 19b). 


APPENDIX  I 


COMPUTER  SIMULATIONS  OF  THE  RADIOMETER  AND  AUTOCORRELATION  TECHNIQUES 


FOR  DETECTING  FH  SIGNALS  IN  GAUSSIAN  NOISE 


1.0  INTRODUCTION 

Theoretical  performances  of  the  autocorrelation  technique  and  the 
conventional  radiometer  approach  in  detecting  frequency-hopped  (FH)  signals 
in  Gaussian  noise,  discussed  in  earlier  sections  of  this  report,  have  been 
verified  through  computer  simulation.  The  simulation  results  are  summarized 
in  this  appendix.  Close  agreement  was  observed  between  theoretical  predic¬ 
tions  and  simulation  results. 

2.0  SYSTEM  MODE! 

The  system  model  studied  here  Is  illustrated  In  Figure  1.1.  The 
total  spread  bandwidth  W  Is  segmented  into  M  B-Hz  bands,  where  M  =  W  /B. 

The  received  signal  from  each  of  these  B-Hz  bands  is  processed  by  a  waveform 
processor  which  determines  whether  or  not  an  FH  signal  is  present  in  that  seg¬ 
ment  of  the  spread  bandwidth  during  the  time  Interval  of  observation.  The 
decisions  made  by  the  waveform  processors  are  accumulated  to  form  the  final 
decision  with  respect  to  whether  or  not  an  FH  signal  is  present. 

Two  forms  of  waveform  processing  are  studied  by  computer  simulation: 
(1)  the  radiometer,  and  (2)  the  correlation  detector.  In  each  case,  the  wave¬ 
form  processor  output  statistics  (mean  and  variance)  are  collected  under  either 
hypothesis  regardless  of  whether  or  not  the  signal  is  present.  The  channel  is 
assumed  to  be  corrupted  by  additive  Gaussian  noise  in  all  cases  under  consider¬ 
ation  herein. 

Based  on  the  waveform  processor-output  statistics,  the  detector  out¬ 
put  signal -to-noise  ratio  (SNR)  defined  here  as 


is  first  computed  as  a  function  of  the  bandwidth-observation  time  product  BT^ 
and  the  input  SNR  yin,  which  is  defined  by 


In  (1-1)  and  (1-2),  the  following  notations  are  used: 


=  mean  value  of  the  waveform  processor  output  when  signal  is  not 
present 

*0MHjD  =  mean  value  of  the  waveform  processor  output  when  signal  is  present 
var(Y|H^)  =  variance  of  the  waveform  processor  output  when  signal  is  absent 
S  =  received  power  of  the  FH  signal  to  be  detected 
B  =  input  bandwidth  of  the  waveform  processor 
Th  =  the  FH  time 

Nq  =  one-sided  noise  power  spectral  density  (PSD)  of  the  receiver. 

The  input  SNR  is  related  to  the  hop  SNR  (Yu)  by  the  BT^  product,  as 

f ol 1 ows : 


rH  s  1J-  '  CBTh)  Yin 


(1-3) 


Also  computed  in  this  appendix  are  the  probabilities  of  detection 
and  false  alarm  CpF/0  eac**  wave^onn  detectors.  They  are  related 

to  the  first  two  moments  of  Y  by: 


1  1 


7 ' 7 


So-  ^0^1 


1 2  var(Y|H03 


1  1 


7 ' 7 


a0  -  tfon  Hi)' 


f2  var(Y|  H^) 


(1-4) 


(1-5) 


M 


In  (1-4)  and  (1-5),  aq  denotes  the  threshold  setting  and  erf  denotes  the  error 
function 


For  each  desired  PpA,  the  probability  of  detection  Pg  can  also  be  determined, 
as  follows: 


o  i  hco 


/2  var (Y  |Hi) 


erf_1(l  -  2PF/0| 


Gaussian  approximations  have  been  made  on  the  waveform  processor  out¬ 
put  Y  in  arriving  at  the  expressions  for  PQ  and  PpA  in  (1-4)  and  (1-5).  This 
is  justified  by  the  fact  that  Y  is  the  sum  of  a  large  number  (-BT^O  of  indepen¬ 
dent  random  variables  (see  next  section)  and  by  invoking  the  central-limit 
theorem. 

3.0  RADIOMETER  SIMULATION 

For  the  case  of  the  radiometer,  the  waveform  processor  performs  the 
following  operation  on  the  received  signal  r(t)  to  obtain  the  detector-output 
variable  Y: 


H  p2(t)  dt 


where  p(t)  is  the  envelope  of  the  received  bandpass  signal  r(t),  given  as 
r(t)  =  cos  wgt  +  n j f t)  cos  u>gt  -  ng(t)  sin  wgt 


In  (1-9),  ojq  is  the  center  frequency  of  the  received  signal  and  nj  and  ng  are 
the  quadrature  components  of  the  bandpass  noise  process  with  a  bandwidth  B. 
Thus,  nj  and  ng  are  independent,  zero-mean  Gaussian  processes  having  one-sided 
noise  PSD  Ng  and  variances  NgB/2.  We  assume  here  that  oig  coincides  with  the 
center  frequency  of  the  BPF  with  bandwidth  B  (see  Figure  1.1)  so  that  there  is 
no  filtering  loss  on  the  received  pulse.  This  assumption  is  not  unduly  restric 
tive  since  most  channels  in  the  bandwidth  B  will  not  experience  excessive  fil¬ 
tering  loss  except  for  the  ones  near  the  band  edges. 


Since  r(t)  is  a  bandpass  process,  it  can  be  expressed  as 


r(t)  =  Im</?  u(t)  e 


I 


(1-10 


where  Im(z)  stands  for  the  imaginary  part  of  z  and  u(t)  is  the  complex  envelope 


u(t)  =  SS  +  nT(t)  -  j  nn(t) 


JVW  J  "QV 

where  j  =  /T.  The  envelope  p(t)  is  given  by 


(1-11 


p(t)  =  ^|u(t)| 


( 1-12; 


Since  n^  and  nq  (and,  thus,  p(t)  also}  are  bandpass  processes  with 
bandwidth  ±B/2,  by  invoking  the  (stochastic)  sampling  theory,  they  can  be  rep¬ 
resented  as  samples  taken  at  time  intervals  separated  by  AT  =  1/B.  In  other 
words,  a  valid  simulation  of  the  integration  process  (1-8),  from  which  the  de¬ 
tector  output  variable  Y  is  derived,  is  the  discrete-time  summation 


Y 


(1-13! 


where  6  =  BTH  is  the  number  of  samples  to  be  accumulated  over  the  hop  time  T^. 
The  discrete-time  simulation  (1-13)  is  equivalent  to  the  actual  process  (1-8) 
in  the  sense  that,  if  the  sample  values  r(k/B)  are  taken  from  the  actual  band- 
limited  analog  signal  r(t),  the  discrete  system  outputs  at  the  sample  time 
instants  (in  this  case,  after  G  =  BT^  samples  have  been  accumulated}  will  be 
identical  to  the  analog  system  output  at  the  same  instants*.  The  noise  samples 
of  nj  and  nq  input  to  the  discrete-time  simulation  model  can  be  taken  to  be 
uncorrelated  and,  thus,  independent.  This  assumes  that  the  impulse  response  of 
the  bandpass  filter  preceding  the  waveform  processor  has  zero  crossings  at  1/B. 
One  BPF  response  which  will  satisfy  this  assumption  is  the  ideal  response 


H(u) 


1  |  U>  —  0)g|  <  B/2 

0  el sewhere 


(1-14 


Reference  A.  Papoulis,  Signal  Analysis,  McGraw-Hill  Book  Company, 
1977,  pp.  25-27. 


This  ideal  response  is  closely  approximated  by  actual  BPF  s  with  sharp  cutoffs 
which  will  probably  be  the  case  in  actual  radiometer  implementations. 

Table  1.1  summarizes  the  theoretical  predictions  on  the  radiometer 
output  statistics  which  are  required  to  compute  SNRQ.  They  are  given  here  as 
functions  of  BT^  and  S/NqB.  The  results  summarized  in  this  table  are  those 
derived  from  the  discrete-time  model,  which  is  of  interest  when  calibrating 
the  accuracy  of  the  simulation. 

Table  1.2  summarizes  the  radiometer  simulation  results.  It  is  ob¬ 
served  that  the  simulation  results  agree  very  well  with  the  computed  results 
from  theory  (with  SNRq  error  <  0.25  dB  in  200  trials). 

Based  on  the  simulated  results  of  the  means  and  variances  of  the 
radiometer  output  under  Hg  or  Hp  the  probability  of  detection  [Pp)  of  FH  sig¬ 
nals  by  the  radiometer  technique  is  computed  by  (1-7),  for  each  desired  proba¬ 
bility  of  false  alarm  [Pp^j ,  as  a  function  of  input  SNR  y^n.  The  results  for 
BT^  =  100  and  BT^  =  1000  are  summarized  in  Figures  1.2  and  1.3,  respectively. 

4.0  AUTOCORRELATION  TECHNIQUE  SIMULATION 

In  the  autocorrelation  technique,  the  waveform  processor  output  Y  is 
characterized  by  the  following  operations: 


G 

■  & 


where  is  the  squared  value  of  the  envelope  of  the  autocorrelation  function 
Rr( x)  of  the  received  signal  r(t)  at  t  =  k/b,  and.a^’s  are  weighting  coeffi¬ 
cients 


Cth  -  W)~  if  k  <  ag 


if  k  >  AG 


where  G  =  BT^  and  A  is  a  positive  fraction  <  1,  which  serves  to  truncate  the 
sum  (1-15)  to  (A  BTh  -  1)  terms.  The  autocorrelation  function  Rr(x)  of  the 
bandpass  signal  r(t)  with  complex  envelope  u(t)  [see  (1-9),  1-10)  and  (1-11)) 
is  given  by 


WW 


Table  1.2.  Comparison  of  Radiometer  Simulation  Results  to  Theory 
(200  Trials) 


Input  SNR 


-20  dB 


■10  dB 


-5  dB 


or  S^muVation  Theory  Simulation  Theory  Simulation  Theory  [simulation 


^ (Y IHi)  4-04 


var  (y|H0)  16.00 
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■20.14  0.00 


3.96  4.00  4.01  |  4.00 
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4.00  4.40  4.41  !  5.27 
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15.96  16.00  15.63  16.00  |  15.96 
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0.25  10.00 
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Theory  Simulation  Theory  Simulation  Theory  Simulation 


4.003  4.00 


4.00 


4.04 


16.00 


4.003  4.00 


4.044  4.126 


16.14  16.00  16.14  !  16.00 


16.20 


-9.73  0.00 


16.69  j 


0.14  ;  10.0 


4.003 


4.407 


16.14 


18.57 


10.06 
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Rr(x)  =  Re|Ru(T)  e  u  j 

where  Ru(x)  is  the  autocorrelation  function  of  u(t): 

ru(t)  =  /  H  u*(t)  u(t-x)dt  (1-18) 

0 

Since  Rr(x)  is  itself  a  bandpass  function  of  x,  the  envelope  of  Rr(x)  Is  | Ru ( T ) | 
Thus,  we  have 

Hk  =  |Rr(l)|2  ■  K(4)|2  •  f"  u*(t)  »(t  -  {j)«  (1-19) 

k/B 

Now,  sincp  u(t),  given  by 

u(t)  =  ^  +  rij(t)  -  J  nQ(t)  (1-20) 

is  bandlimited  to  B,  the  simulation  of  (1-19)  can  again  be  performed  by  a 
discrete-time  model  similar  to  that  used  in  the  radiometer  simulation.  The 
discrete-time  model  is  characterized  by  the  following  equations: 

G  2 

Wk  =  4  I  u*[£/B]u[(£-k)/B]  (1-21) 

°  £=k 

G 

Y  =  1  a.W.  (1-22) 

k=l  K  K 

and 

)  (A)  If  k  ^  AG 

ak  =  <  (1-23) 

/  0  if  k  >  AG 


The  mean  values  of  and,  thus,  Y,  under  either  hypotheses  Hq  or  Hj, 
can  be  computed  relatively  easily  from  the  discrete-time  equations  (1-21),  (1-22) 
and  (1-23).  These  theoretical  results  are  summarized  in  Table  1.3. 

The  variance  of  Wk  and,  thus,  Y,  is  difficult  to  compute  exactly,  how¬ 
ever,  due  to  the  complexity  of  W^.  (Computations  of  expected  values  of  products 
of  eight  Gaussian  random  variables  are  required  here.)  Under  hypothesis  HQ, 
the  sum  in  of  (A-21)  can  be  written  as 

Q 

1  “*(1)  "pry  *  (Nn  +  nqq)  +  J(niq  -  nqi)  (‘-“l 


where  Njj,  Nqq,  Njq  and  Nqj  are  the  sums  defined  by  (with  k  >  1): 


Ntt  »  S  nT(£)nT(£-k) 
£=k  1  1 


Nqq  =  2  nQ(£)nQ(£-k) 


NIQ  =  2  nI(£)nQ(£-k)  ; 


Nqi  *  2  nQ(£)nj(£-k)  (1-25) 

£— k 


The  above  noise  terms  are  zero  mean,  with  variance  (G-k)(NoB/23  .  They  are  un¬ 
correlated.  For  k  «  G,  they  are  sums  of  a  large  number  of  independent  random 
variables.  Thus,  for  k  «  G,  they  are,  to  a  good  approximation,  Gaussian.  We 
can  therefore  write  W^  as 


m 


2  .  2' 


\  +  ^ 


(1-26) 


where  and  ng  are  given  by 


nA  =  Nn  +  Nqq 


nB  =  NIQ  ‘  NQI 


(1-27) 


and  are,  to  a  good  approximation,  independent  Gaussian  random  variables  them¬ 
selves,  with  mean  zero  and  variance  2(G-k)QioB/2)2.  Thus,  Wk  is  a  squared  Ray¬ 
leigh  random  variable,  to  a  good  approximation.  The  mean  square  of  Wk  is  equal 
to 
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4k2)  3  32(G-k)2(-fj  (})  =  2(1)  (G-k)2CN0B)4  (I-2£ 

The  variance  of  Wk  under  HQ  can  thus  be  written  as  (see  Table  1.3): 


varpy  .  *(uk2)  -  £(>y]2  -  ($)*  CG-k)2  ptoB)' 


(I-2S 


It  can  also  be  shown  that  W^'s  are  uncorrelated  (thus,  independent).  Based  on 
this  approximation,  the  variance  of  Y  in  (1-22)  can  be  written  for  Hq  as: 


varQT  fHoD  =  (|*q<)  var(W(0  =  ^  (^)  CN0B)4 


( I-3C 


Combined  with  the  results  in  Table  1.3,  we  finally  obtain  an  approximated  for¬ 
mula  for  the  autocorrelation  detector  output  SNR  (for  the  discrete-time  model): 


[xGS2  +  2SN„B  f  WLil 


Again,  these  results  are  summarized  in  Table  1.4  for  easy  reference.  For  cases 
when  GYin  »  2/x  £nCl/(l-x)),  where  0  <  X  <  1  (i.e.,  large  BT^  products),  the 
improvement  of  the  autocorrelation  technique  over  the  radiometer  approach  is 
given  by  the  ratio 


correlation 


radiometer 


x(l-x)  (Gy J: 
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Table  1.4.  Autocorrelation  Technique  Detector- Output  SNR 
Theory  for  the  Discrete-Time  Model 
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2 

It  is  of  interest  to  note  that,  for  these  cases,  the  optimum  choice  of  x  (the 
summation  truncation  factor)  is  approximately  =  1/2.  For  X  =  1/2,  the  improve¬ 
ment  of  the  autocorrelation  technique  over  the  radiometer  approach  is  given  by 

1  ?(G*1„)2  O'34* 

x=l/2 

(The  exact  result  is  given  in  Table  1.4  for  all  X.) 

Tables  1.5  and  1.6  compare  the  simulation  results  to  the  theoretical 
predictions  with  Gaussian  approximations  on  variance  calculations.  The  follow¬ 
ing  observations  can  be  made: 

(1)  Simulation  results  agree  exceptionally  well  with  theory  for  the 
mean  values  of  the  correlator  detector  outputs  for  all  values  of  G,  Yin  and  x. 
This  is  true  since  the  theory  for  this  part  of  the  calculations  is  exact. 

(2)  Simulation  results  show  that  the  variance  of  the  detector  output 
Y  under  H|  can  be  quite  a  bit  larger  than  that  under  Ho  for  high  Yin  values. 

Also,  simulation  results  on  the  variance  of  Y  is  generally  larger  than  the  ap¬ 
proximated  theoretical  predictions  for  x  >  0.1,  where  Gaussian  approximations 
used  in  computing  the  variance  of  Y  began  to  lose  validity. 

Based  on  the  autocorrelation  detector  mean  and  variance  simulation 
results,  the  probability  of  detection  versus  Y^n  for  various  PFA  settings  are 
computed.  These  results  are  summarized  in  Figures  1.4  and  1.5.  The  Improve¬ 
ment  in  performance  of  the  autocorrelation  technique  as  BTH  increases  from  100 
(Figure  1.4)  to  1000  (Figure  1.5)  is  apparent.  The  truncation  factor  x  in  Fig¬ 
ure  1.4  is  0.5,  while  that  in  Figure  1.5  is  0.1.  The  performance  is  relatively 
insensitive  with  respect  to  X  as  long  as  X  is  >  0.1  and  is  <  0.8.  This  X  depen¬ 
dence  is  illustrated  in  Figure  1.6,  where  the  probability  of  detection  Is  plotted 
against  X  for  various  pfa(  10"2,  10"4  and  10-6)  settings  and  at  various  Yin  (-5, 
-10,  and  -20  dB)  for  the  case  BT^  =  100.  Thus,  we  can  regard  Figures  1.4  and 
1.5  to  be  the  representative  performance  of  the  autocorrelation  technique  for 
BTh  =  100  and  1000,  respectively. 
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Table  1.5.  Comparison  of  Autocorrelation  Detector  Simulation  Results 
to  Approximated  Analysis 

(a)  BTh  =  100;  X  *  0.01 


Yin  (dB) 

-20  dB 

-10  dB 

-5  dB 

Theory 

or  Simulation 

Theory 

Simulation 

Theory 

Simulation 

*C»|h0D 

0.38 

0.378 

0.38 

0.392 

0.38 

0.384 

0.39 

0.390 

0.82 

0.826 

m 

4.24 

«r(V|H„3 

0.016 

0.0251 

0.016 

0.0215 

0.016 

0.021 

- 

0.0274 

- 

0.0235 

- 

4.1 

Autocorrelation 
Detector 
SNRQ  (dB) 

-21.1 

-22.5 

10.74 

9.42 

29.64 

28.57 

(b)  BTh  =  100;  X  =  0.5 


Yin  <dB> 

-20  dB 

-10  dB 

-5  dB 

Theory 

or  Simulation 

Theory 

Simulation 

Theory 

_ 

Simulation 

Theory 

Simulation 

■»CYl"o3 

D 

2.76 

O 

2.69 

m 

2.695 

*on  Hi3 

2.786 

2.84 

5.22 

5.22 

24 

23.75 

varCY|  Hq) 

0.156 

wa 

0.156 

0.633 

0.156 

0.708 

varCYl  Hi3 

- 

0.665 

- 

6.22 

91.4 

Autocorrelation 
Detector 
SNRq  (dB) 

_ 

-14.56 

-19.5 

16.0 

10.04 

34.6 

28 

Table  1.5.  Comparison  of  Autocorrelation  Detector  Simulation  Results 
to  Approximated  Analysis 

(c)  BTU  =  1000;  A  =  0.1 


Yin  <dB> 


-20  dB 


-15  dB 


-10  dB 


Theory 

or  Simulation 

Theory 

Simulation 

Theory 

Simulation 

Theory 

Simulation 

*aiHo3 

0.417 

0.419 

0.417 

0.419 

■ 

0.417 

0.42 

^ctihj 

0.47 

0.472 

0.84 

0.863 

4.46 

4.57 

«ror|H03 

1.76  xlO"3 

2.61xl0"3 

1.76xl0-3 

2.6  x 10"  3 

1.76  xl0“3 

2.61xl0~3 

var(Y|Hj3 

4.97xl0“3 

- 

5.34  xlO'2 

- 

1.41 

Autocorrelation 
Detector 
SNRq  (dB) 

1.16 

0.17 

20.06 

18.76 

39.7 

38.2 

> .  ■ 


3-u  RELATIVE  COMPARISON  OF  RADIOMETER  AND  AUTOCORRELATION  TECHNIQUES 

IN  DETECTING  FH  SIGNALS 

Relative  performances  of  the  radiometer  and  autocorrelation  techniques 
in  detecting  FH  signals  can  be  compared  two  different  ways:  (1)  their  detector 
output  SNR  (as  defined  in  (1-1)),  or  (2)  their  operating  characteristics  Pp  and 
Pp/\  (in  particular,  in  terms  of  the  required  input  SNR  Y^n  for  the  specified  Pp/\ 
and  PD}. 

Table  1.6  compares  the  two  techniques  in  terms  of  their  respective  out¬ 
put  SNR's  as  functions  of  BT^  and  Y^.  As  shown  in  Table  1.6,  the  SNRq  of  the 
autocorrelation  technique  is  significantly  improved  over  that  of  the  radiometer, 
especially  in  large  BTH  cases  and  with  high  Y^n.  In  typical  systems,  BT^  is 
large  while  Y-n  is  small  (BT^  =  1000,  Y^n  =  -20  dB).  In  that  case,  the  autocor¬ 
relator's  output  SNRq  is  still  z  10  dB  over  that  of  the  radiometer.  As  shown  in 
Table  1.6,  the  theoretical  prediction  on  SNRq  improvement  agrees  very  well  with 
simulation  results.  Thus,  this  relationship  can  be  described,  to  a  very  good 
approximation,  by  the  equation  (see  Table  1.4): 

rkin  autocorrelator 


c,,n  radiometer 


Equation  (1-35)  indicates  that  this  improvement,  for  a  fixed  X,  is  essentially 

-C3  YiJ2>  or  V* 

A  more  detailed  evaluation  of  the  relative  performances  of  these  two 
approaches,  which  is  probably  more  useful  in  terms  of  actual  system  design,  is 
to  compare  their  respective  required  input  SNR  (Yin)  for  a  specified  set  of  Pp 
and  Pq.  These  data  can  be  extracted  from  Figures  1.2,  1.3,  1.4  and  1.5.  The 
results  are  plotted  in  Figures  1.7  and  1.8.  The  relative  improvements  of  the 
autocorrelation  approach  over  that  of  the  radiometer,  in  terms  of  required  Y-n 
to  give  specified  Pcfl  and  Pn,  are  summarized  in  Table  1.7. 


vvy* 


Mnltftl 


Table  1.6.  Relative  Comparisons  of  Detectc.  „,r _ 

(a)  BTh  =  100;  X  *  0.1 


Simulation 

Theoretical  Prediction  | 
of  Autocorrelation  SNRq] 
Improvement  Over  ' 

Radiometer 

m 

Radiometer 
SNRq  (dB) 

Autocorrel ation 
SNRq  (dB) 

Improvement  of 
Autocorrelation 
over  Radiometer 
in  SNRo 

.  ..j 

-20 

-20.14 

-22.5 

-2.36  (dB) 

-1.1  (dB) 

-10 

0.25 

9.42 

9.17 

10.74 

-5 

9.98 

28.57 

18.6 

19.64 

(b)  BTh  =  1000;  X  =  0.1 


1 

Simulation 

Theoretical  Prediction 
of  Autcorrel ation  SNRq 
Improvement  Over 
Radiometer 

Radiometer 
S-NR0  (dB) 

Autocorrelation 
SNRq  (dB) 

Improvement  of 
Autocorrelation 
Over  Radiometer 
in  SNRq 

-20 

-9.73 

0.17 

9.9  (dB) 

11.6  (dB) 

-15 

0.14 

18.76 

18.62 

20.06 

-10 


10.06 


38.2 


28.14 


29.7 
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Table  1.7.  Relative  Improvement  of  the  Autocorrelation  Technique  Over  the 
Radiometer  Approach  in  terms  of  Required  Input  SNR 
For  Desired  PQ  and  PpA  Performances 


bth 

PFA 

10'5 

100 

10"3 

10" 1 

10"5 

1000 

10"3 

10"1 

Improvement  of  Autocorrelation  Technique 
over  Radiometer  (dB) 


8.4 

8.5 

7.3 

7.5 

5.4 

5.5 

where 


APPENDIX  J 

STATISTICAL  CHARACTERIZATION  OF  NOISE  PROCESSES 


awn 


nj(t)  AND  ng(x) 


n£(-r)  *  ak(-r)  cos  “kT  ”  S1n  “k1*  . M  ( J~1 ) 


nn<T) 


"  (jj  6k,T>)  cos  “O'  +  (k2j  6k(t)) 


sin  ujqT 


(J-2) 


“k(t> 

=  y  H  c(t)  COS^Aa)kt  -  ♦fcjdt 

T 

(J-3a) 

“k<T) 

3  /"  c(t)  s1n^Au>kt  -  4kjdt 

T 

(J-3b) 

8k(T> 

=  /  "  C(t-r)  COS ^Aukt  -  4k) dt 

T 

(J-3c) 

<8«*> 

*  f  H  c(t-x)  s1n|Au>kt  -4k)dt 

(J-3d) 

with  =  uq  -  u)k.  Equations  (J-l)  through  (J-3)  are  also  given  by  (5-42)  and 
(5-43).  First,  since  c(t)  Is  a  random  code  with  <f{c(t)}  =  0,  it  follows  that  the 
noise  processes  in  (J-l)  and  (J-2)  are  zero  mean,  and  that 


v>rk(T,f =  vari6Jk(T)j  ■  'j(«k<T))2j 


(0-4) 


J2 


We  shall  first  establish  (5-47),  i.e.,  claim  (a)  In  section  5.2.2. 

We  have  from  (J-3a)  that 

T  r\  r  “1  N"1  /*(n+l)T_  r-  -| 

ak(mT0  =  /  cos  OHOt  -  4k  dt  s  I  cn+l  J  cos|^Auk)t  -  4k  dt 

V  L  -1  n=m  <  L  -* 


2  cn+1  1 


n+1  n+1 


( J-5) 


where  cn+1  is  the  (n+l)tn  code  chip  and 


'  cos[cAailOt  '  4k] dt 


( J-6) 


Let  c  =  [c, . cN].  Then, 

2[  ■  \U|J(W 


N-l  N-l 


\j^l*kjn?mn^nCn+1  C"‘+1  *"+1  X",+1 


(J-7) 


r{cn+l  cn'+l}  =  5nn' 


1  ,  If  n*n' 

0  ,  otherwise 


which  reduces  (J-7)  to 


rj(“kOnT0)2}  ■ 


jJL  1,1+1 


N-i  (  9 

^  *4i,  )  *n+l 

n-m  ¥k  « 


( J-8) 


However, 


\{icos  *k 


r(n+l)T 

(X n+1)Tc  \21 

J  cos^Au^dt  +  sin  $k 

J  c  sin(Aw|ct)dt)  > 

"Tc 

"Tc  /  J 

•  (1) 


/  /■(n+1>Tc  \2  If 

I  J  COS£Aii)|<t}dt\  +  I  J 

VTc  /  \"TC 


(n+l)T 


slnfAw^dt 


J 


(J-9) 


But 


J *  cos^Awfct^dt  =  |-^jsinx  =  ^j|lln(Au>)b  -  sin(Aoj)aJ  (JlO-a) 


( Au>)a 


and 


J"  slnfAufcQdt  *  |^j|cos(Aa))a  -  cos(Aw)bj 


(JlO-b) 


Upon  inserting  (J-10)  into  (J-9)  and  combining  with  (J-8),  it  follows  that 

\2 


r|“kO"Tc3)2}  ■ 


2  (A. 


1  N"1  n  \2 

- — jr  I  (sln(n-H)  ^Au)k3Tc  -  sin  n(Au>|0Tj 

iukD  n=Tn  L_  / 

+  |cos  n  (4*^)1^  -  costn+lJ^Aw^^^J 

fii 


N-l 


2(>|0 


" — 5-  i  2-2  cos(Aw./)t1 
wirV  n»m  L  KJ  c] 


2  sir 


CHJ 


(N-m) 


or 


a. 


itCmTc}) 


2)  (N-m)T 


^Sa2 


OO  T 


(J-ll) 


which  is  (5-47). 


J4 


We  shall  now  proceed  to  prove  claim  (b),  namely,  that  each  noise 
process  n£(T);k=l,...,M  and  n^(t)  have  uncorrelated  (and,  by  the  Gaussian  assump¬ 
tion,  independent)  inphase  and  quadrature  components.  From  (J-3a)  and  (J-3b),  we 
can  write 


°k(T)  = 


cos 


+  Ak  sin 


and 


“k<T> 


■  Ak  cos  <j»k  -  Aj[  sin  <j>k 


where  we  define  the  random  variables  Ak  and  Ak  as 

Ak  =  J  H  c(t)  cos  Acukt  dt 


and 


AQ 

Ak 


c(t)  sin  Aa>kt  dt 


Thus,  from  (J-12), 


(J-12a) 


(J-12b) 


(J-13a) 


(J-13b) 


^H(T)ak(T>f  *  'cKicK(t)^(,)I|  ■  'cjAkAk\icos2vsin\l 

*  BAk  )2  -  (Ai!f]-\ls1n*k  cos*kf{ 


“  * 


2*»f 


'ckf  -  K)1  , 

2*kf +  —  2— — ^  \lsin  2*l 


I  sin  2<frkj  *  0  (J-14) 


which  proves  the  uncorrelatedness  of  ak(x)  and  ak(x).  In  a  similar  manner, 

Bk(x)  and  Bk(x)  are  shown  to  be  uncorrelated.  Furthermore,  Bk(x)  and  b^(t)  for 
kjhn  are  also  uncorrelated  since  they  correspond  to  independent  phases  <{>k  and  $m. 
This  implies  that  the  random  variables 

M  t.  .  !!  o.  . 


which  constitute  the  inphase  and  quadrature  components  of  n^(T),  are  uncorrelated. 
This  concludes  the  proof  of  claim  (b). 

The  proof  of  claim  (c)  is  based  on  the  simple  fact  that 


<?jajj(-c)  0^(t)[  -  0  ;  k?*m;j=I,Q 


(J-15) 


since  $k  and  are  independent  ([see  (J-3)[).  Finally,  in  order  to  substantiate 
claim  (d),  we  look  at  n£(t)  at  t  =  0  and  t  *  Tc-  Then, 


“k(0)  =  r  ^  c(t)  COS^Aa^t  -  fyjdt  *  J"  C  c(t)  COS^Ao^t  -  ^jdt  + 


5  “kFcD 


(J-16) 


since  the  first  term  ([integral  from  0  to  T  )  is  negligible  compared  to  the 
(N-l )TC  integration  involved  in  ak(Tc[).  Thus,  <^(0)  is  almost  equal  to  (i.e., 
very  correlated  with)  a£(Tc[),  and  the  same  holds  for  ajj(-t).  It  follows  that 
n£(x)  is  a  highly  correlated  noise  with  respect  to  the  Tc  time  Interval.  On 
the  other  hand. 


jJ(°)  *  y  c  c(t)  cos(Aaikt  -  ♦|c)dt 


(J-17a) 


while 


SkPc)  ■  f  H  c([t  -  Tc3  cos^Aukt  -  «(»k Jdt 
c 

s  J  H  c  c(t')  COS^Aoikt'  +  Ao>kTc  -  $k)dt‘ 
0 

S  ek(0)  cos(AaikTc)  -  pj(0)  sin(AukTc) 


$ 

V 

V, 

V, 


(J-17b) 


v*  v1 


‘  V*  .VV  .  v.vvy.v  '  V.*. 


.V.V.VC.V.V.  K. 1  .  _V7sw 


Thus,  the  degree  of  correlation  between  n^( 0)  and  0k(Tc}  depends  on  the  phase 
shift  AwkTc.  If,  for  instance,  s  0  or  */2  or  then  Bj^0)  an<^  ek  Tc 

are  highly  correlated  while,  if  C^0"^k3^c  =  1/4  or  3/4,  then  &k(Jc3  =  ^(0), 
which  is  uncorrelated  with  3^(0)!  Since  similar  rules  hold  for  0^(t),  we  con¬ 
clude  that  the  correlation  time  of  the  noise  n^-r)  varies  with  (Af|QTc.  For  the 
particular  case  of  one  interfering  tone  on  the  signal  frequency,  i.e.,  f,  =  f0, 
it  is  concluded  that  nj(0)  a  n^T^,  resulting  in  a  highly  correlated  noise 
(which  is  a  favorable  case  here). 


i 


APPENDIX  K 


PROOF  OF  EQUATION  (5.52) 

From  (5-50),  we  have  that 

4  *  E  +  2ti"1(1  *  ^  al.no"3 

where 

al,norm  *  ft)  «  *1  l  S 

Clearly,  ^(aj  norm)  =  0;  which,  upon  substitution  in  (K-l), 
(5-52a).  Furthermore,  from  (K-l)  and  (5-52a), 

var{i|Hi}  ’  var)2CSTl02  l{W  ai,no™|  =  4Pth34  rf1  w)a 

-  4™4  vf 1  ft)  (l)'N 

since  varjcos  =1/2  and  the  cn's  are  independent,  so  that 

Var|n=l  Cn|  ’  N,var  {cn}  =  N 

If  we  recall  that  N  -  TH/TC,  we  conclude  from  (K-3)  that 

1  "  2^Th^ 

varr|Hn  — 


(K-l) 

(K-2) 

imnediately  yields 

.norm} 

(K-3) 


(K-4) 


APPENDIX  L 


RADIOMETER  PERFORMANCE  FOR  THE  FH/DS  CASE  WITH  RANDOM  TONE  INTERFERENCE 
The  radiometer  output  Yra(j  =  y (0)  can  be  written  s 


Y 


rad 


sth  +  miith  +  2>™ 


"M 

k=l 


(L-l) 


where  Ij  is  the  power  per  tone,  so  the  total  power  for  M  tones  is  I  *  MIj.  In 
our  formulation,  M  is  a  random  variable,  taking  values  in  the  set  [1,2,...  ,^ax3 
with  equal  probability  in  each  hop.  Conditioned  on  M,  the  mean  and  variance  of 

vrad  are 

■  STH  +  MI1TH  ■  (S  +  hil)th  <l-2> 

and 


where  we  have  used  the  fact  that  oj(o)  and  a^(0)  are  uncorrelated  for  Mm.  (See 

k  ni  » 

Appendix  J.)  Furthermore,  the  second  moment  of  atj^(0)  is  given  by  equation  (5-47) 
which,  upon  substitution  in  (L-3)  for  m=0  yields 

var{Vrad|M}  ‘  2SIj  NTC2  Jj  sa2[<f0-fOTc]  *  2SI1  TH  Tc  L(H)  C"4' 

where 

KM)  «  jt  Sa2  [<f0-fk]Tc]  (L-5) 


and  the  fact  that  Tu  =  NT,,  has  been  used. 


r 

9 


3 

5 


9 


•2 
L 

L 

E* 


Under  the  Gaussian  assumption  for  Yrad  (conditioned  on  M),  we  can 
easily  derive  the  conditional  detection  probability  under  hypothesis  Hj  as 


p5|H  ■  Pr{Yrad>Yo|M'Hlf  ' 


-  (S+MIl)T| 


(L-6) 


where  YQ  is  the  detection  threshold  and  equations  (L-2)  through  (L-5)  have  been 
used.  Clearly,  under  Hg,  the  first  and  third  terms  in  (L-l)  disappear,  so  that 


mlTH 


under  Hr 


(L-7) 


with  probability  1.  This  implies  that  the  conditional  false-alarm  probability 
is  two  valued: 


1  ,  if  MIjTH  >  Yq 
0  ,  if  MI^  <  Yq 


(L-8) 


In  deriving  (L-8),  we  have  implicitly  assumed  that  the  detector  decides  HQ 
whenever  the  radiometer  output  Yrad  is  exactly  equal  to  the  threshold  Yq. 

If  we  now  insert  the  equidistribution  assumption  about  M  between  1  and 
,  we  can  average  (L-6)  or  (L-8)  with  respect  to  M  in  order  to  obtain  the 
overall  Pp3d  and  pj^*.  Thus, 


PDad  *  ■ 


1  V  Qpo  •  (S*M11>TH 

£  M=1  ^2SI1ThT<;L(M) 


(L-9) 


while 


prad  _  .  ^0 

PA  '  ‘‘t 


(L-lOa) 


where  My0  is  the  largest  integer  not  exceeding  Yq/I^,  i.e.. 


/  * 


(  Mmax  lf  >  %ax 


In  (L-lOb),  1_*J  denotes  "integer  part  of."  Furthermore,  for  a  desired  Pp^a  level, 
Myn  can  be  calculated.  Then,  YQ  should  be  chosen  at  the  lowest  allowable  level 
in  order  to  maximize  .  Clearly,  this  choice  is 


Y0  ■  %'lTH  =  ^naxl1  -  PHd)1lTh 


(L-ll) 


We  note  that.  In  the  case  treated  here,  i.e.,  no  thermal  noise,  P„  can  assume 

*  **  M9/| 

only  discrete  (quantized)  values,  as  dictated  by  (L-lOa).  In  other  words,  PFA 
can  only  be  of  the  form 


>rad  =  k  1 


;  k-0,1,...,  Mmflx 


(L-12) 


We  can  now  combine  (L-9)  with  (L-ll)  in  order  to  get  the  radiometer's 
operating  characteristic  as 


As  an  interesting  side  result,  let  us  consider  the  limiting  case  where 
i.e.,  the  DS  code  rate  is  much  larger  than  the  FH  rate.  Since  this  would  cor¬ 
respond  to  Tc  -*■  0  £for  a  fixed  T^.  we  conclude  that 

M  M 

11m  2  Sa2Rf0-fOTd  "  I  11"  Sa2Rf0-fOTcl  *  M  (L-16) 

N-+~  k-1  L  CJ-  k=l  T  L  c-l 

c 

2 

since  lim  Sa  (x)  *  1.  Therefore,  the  whole  denominator  of  the  Q(>)  functions  in 
x^O 

(L-15)  tends  to  zero, 

2yi*M 

lim  — n —  =  0 


which  implies  that 


lim  Q(  • ) 
N-*-+“ 


1  ,fH>ir»x(l  -^d)  -n 
0 


(L— 17) 


Therefore,  from  (L-16)  and  (L-17), 


lim  P 


Mmax  -  MD  +  1 


(L-18) 


where  Mq  is  the  minimum  M  such  that  M  >  Mmflx(l  -  Pj?^)  -  Yj-  Equation  (L-18) 
can  be  rewritten  in  the  final,  more  compact,  form,  i.e., 


lim  Pi 


max  (l^maxj1  "  FFA  )  *  YlJ  +  *•*  )  "* 


So,  for  instance,  for  M_.w  =  2,  it  is  easily  derived  from  (L-19)  that 

iTlaX 


lim  Pi 


1/2  if  PpA  *  0 

1  if  PFA  -  1/2 


as  expected.  More  generally,  for  an  arbitrary  Mmw, 

max 

lim  -  P”d  +  ' 

».*.  0  FA 


The  case  treated  in  the  text,  f^-1,  is  slightly  different  from 
that  which  is  presented  in  this  appendix.  In  particular,  let  us  assume  that 
in  each  hop,  either  none  or  one  Interfering  tone  can  be  present  with  equal 
probability,  i.e., 

Pr[M*0]  *  Pr[M*l]  -  1/2  (L 


Then, 


"maxKT)-  T,-i 


[2  N* 


1  -  Q 


where 


(L-22) 


and  we  have  assumed  that  f^  * fg  (interference  centered  on  signal  frequency). 
Here,  M  =1.  Furthermore,  if  we  wish  to  have  p£f^  *  0,  we  should  set  the 

max  r  H 


1  -  tQ 


.  1  /prad 
7\PD|H-1 


D I M*0  j 


If  Yj  >  1 


If  Yj  <  1 


APPENDIX  M 

SIMULATION  OF  THE  FH/DS  CASE  IN  RANDOM  TONE  INTERFERENCE 


The  received  signal  is  given  by 

r(t)  =  y'ZS  c(t)  cos  o)Qt  +  JZT  cos^ujt  + 

Recall  from  the  text  that 


N  -  Rc/Rh 


(M-l) 


(M-2) 


with  Rc  representing  the  code  rate  and  R^  being  the  hop  rate.  The  real-time 
autocorrelator  is  specified  by 


y(t)  .  p 


r(t)  r(t-t)dt 


0  <  r  <  Tu 


(M-3) 


Letting  *  uj  with  unknown  phase  difference  ^  and  putting  (M-l)  in  (M-3)  yields 
y(x)  =  S  yc(x)  cos  oigT  +  l(TH-T)cosa»jT  +  >^§T  |n“(x)  +  n^(x)j  (M-4) 

with 


nf(x) 

*  oj(x)  COS 

u»l  -  Oj(t)  sin  aijx 

(M-5) 

n§(t) 

=  ej(x)  cos 

“oT  +  s^n  “o1 

(M-6) 

=  cos  / 

•Tu 

H  c(t)dt 

(M-7) 

«/ 

T 

aJ(T) 

*  -sin  J 

^  c(t)dt 

(M-8) 

X 

e!(x) 

■  cos  ♦,  / 

r  Tu 

n  c(t-r)dt 

(M-9) 

